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Abstract 

The use of data-random graphs in statistical testing of spatial patterns is introduced recently. In this 
approach, a random directed graph is constructed from the data using the relative positions of the points from 
various classes. Different random graphs result from different definitions of the proximity region associated 
with each data point and different graph statistics can be employed for pattern testing. The approach used 
in this article is based on underlying graphs of a family of data-random digraphs which is determined by a 
family of parameterized proximity maps. The relative edge density of the AND- and OR-underlying graphs 
is used as the summary statistic, providing an alternative to the relative arc density and domination number 
of the digraph employed previously. Properly scaled, relative edge density of the underlying graphs is a 
[7-statistic, facilitating analytic study of its asymptotic distribution using standard (7-statistic central limit 
theory. The approach is illustrated with an application to the testing of bivariate spatial clustering patterns 
of segregation and association. Knowledge of the asymptotic distribution allows evaluation of the Pitman 
asymptotic efficiency, hence selection of the proximity map parameter to optimize efHciency. Asymptotic 
efficiency and Monte Carlo simulation analysis indicate that the AND-underlying version is better (in terms 
of power and efficiency) for the segregation alternative, while the OR-underlying version is better for the 
association alternative. The approach presented here is also valid for data in higher dimensions. 
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1 Introduction 



Classification and clustering have received considerable attention in the statistical literature. In this article, a 
graph-based approach for testing bivariate spatial clustering patterns is introduced. The analysis of spatial point 
patterns in natural populations has been extensively studied and have important implications in epidemiology, 
population biology, and ecology. The patterns of points from one class with respect to points from other 
classes, rather than the pattern of points from one-class with respect to the ground, are investigated. The 
spatial relati onships amon g tw o or more classes h ave important implications especially for plant species. See, 
for example, [pieio^ ilMll) and lDixonI (|l994 12002|) . 

The goal of this article is to derive the asymptotic distribution of the relative edge density of underlying 
graphs based on a particular digraph family and use it to test the spatial pattern of complete spatial randomness 
against spatial segregation or association. Complete spatial randomness (CSR) is roughly defined as the lack 
of spatial interaction between the points in a given study area. Segregation is the pattern in which points of 
one class tend to cluster together, i.e., form one-class clumps. In association, the points of one class tend to 
occur more frequently around points from the other class. For convenience and generality, we call the different 
types of points "classes" , but the class can be replaced by any characteristic of an observation at a particular 
location. For exampl e, the pattern of spatial se gregation has been investigated for plant spec i es (Di ggle 1983), 
age classes of plants (|Hamill and Wrighd ()l986f )V and sexes of dioecious plants (jNanami et al 



In recent years, the use of mathematical graphs has also gained popularity in spatial analysis (jRoberts et al 
In spatial pattern analysis graph theoretic tools provide a way to move beyond Euclidean metrics for 
spatial analysis. For example, graph-based approaches have been proposed to determine paths among habitats 
at various s cales and dispersal movement distances, and balance data requirements with information content 
(|Fall et all |2007[ )). Although only recently introduced to landsca pe ecology, gra ph theory is well suited to 
ecological applications concerned with connectivity or movement ((Minor and Urban (2007)). However, conven- 
tional graphs do not explicitly maintain geographic reference, reducing utility of other geo-spatial information. 



Fall et al.l (|2007l ) introduce spatial graphs that integrate a geometric reference system that ties patches and paths 



to specific spatial locations and spatial dimensions thereby preserving the relevan t spatial inform ation. After a 



graph is constructed using spatial data, usually the scale is lost (see for instance, Su et al. ( 20071 )). Many con- 



cepts in spatial ecology depend on the idea of spatial adjacency which requires information on the close vicinity 
of an object. Graph theory conveniently can be used to express and communicate adjacency information allow- 
ing one to compute meaningful quantities related to spatia l point patte rn . Adding vertex and ed ge properties to 
graphs extends the problem domain to network modeling (|Keittl (120071 ) 1. IWu and Murravl (|2008[ ) propose a new 
measure based on graph theory and spatial interaction, which reflects intra-patch and inter-patch relationships 
by quantifying contiguity within patches and potential contiguity among patches. [Friedman and Rafskv (|l983l) 
also propose a graph-theoretic method to measure multivariate association, but their method is not designed 
to analyze spatial interaction between two or more classes; instead it is an extension of generalized correlation 
coefficient (such as Spearman's p or Kendall's t) to measure multivariate (possibly nonlinear) correlation. 

A new type of spatial clustering test using directed graphs (i.e., digraphs) which is based on the relative 
positions of the data points from various classes has also been developed recently. Data-random digraphs are 
directed graphs in which each vertex corresponds to a data point, and directed edges (i.e., arcs) are defined in 
terms of some bivariate function on the data. For exa mple, nearest n eighb or digraphs are defined by placing 



an arc between each vertex and its nearest neighbor. iPriebe et al.l (|200lh introduced the class cover catch 
digraphs (C CCDs) in R and gave t he exact and the asymptoti c distribution of the d omination number o f 
the CCCDs. IPeVinnev et al.l (l2002l ) . [Marchette and Priebd (|2003h . IPriebe et all (|2003al ) . [Priebe et al.1 (|2003bl ). 
and DeVinnev and Priebd (j2006[ ) applied the concept in higher dimensions and demonstrated relatively good 
performance of CCCDs in classification. Their methods involve data reduction (i.e., condensing) by using 
approximate minimum dominating sets as prototype sets (since finding the exact min imum dominating set 
i s an NP-hard problem in general — e.g., for CCCD in multiple dimensions — (see DeVinnev and Priebd 
(|2006D l. Furthermore the exact and the asymptotic distribution of the domination number of the CCCDs are not 
analytically tractable in m ultiple dimensions. For the dom ination number of CCCDs for o ne-dimensional data, 
aSLLN result is proved in DeVinnev and Wierman ( 20031 ) , an d this result is extend e d by Wierman and Xiand 
( 20081 ): fur thermore, a generaliz e d SLL N result is provided by ierman and Xiand (l2008|), and a CLT is also 
proved by IXiang and Wierman ( 2009f l. The asymptotic distribution of t he domination nu mber of CCCDs 
for non-uniform data in R is also calculated in a rather general setting (jCevhanl (|2008l )). ICevhanI (|2005[ ) 
generalized CCCDs to wh at is called proximity catch digraphs (PCDs). The fir st PCD family is intr oduced by 
Cevhan and Priebd (|2003l) : the parametrized version of this PCD is developed bv lCevhan et al.l (j2007t ) where the 



relative arc density of the PCD is calculated and used for spatial pattern analysis. Cevhan and Priebe ( 20051 ) 
introduced another digraph family called proportional edge PCDs and calculated the asymptotic distribution 
of its domination number and used it for the same purpose. The relat ive arc density of this PCD family is 
also computed and used in spatial pattern analysis ( Cevhan et al. ( 2006[ )). ICevhan and Priebe ( 20071 ) derived 
the asymptotic distribution of the domination number of proportional-edge PCDs for two-dimensional uniform 
data. 

The underlying graphs based on digraphs arc obtained by replacing arcs in the digraph by edges based on 
bivariate relations. If symmetric arcs are replaced by edges, then we obtain the AND-underlying graph; and if all 
arcs are replaced by edges without allowing multi-edges, then we obtain the OR- under lying graph. The statistical 
tool utilized in this article is the asymptotic theory of [/-statistics. Properly scaled, we demonstrate that the 
relative edge density of the underlying graphs of proportional-edge PCDs is a [/-statistic, which has asymptoti c 
normality by the general central limit theory of [/-statistics. For the digraphs introduced bv lPriebe et al. (2001), 
whose relative arc density is also of the [/-statistic form, the asymptotic mean and variance of the relative density 
is not analytically tractab l e, due to geometric diffic ulties encountered. H owever, for the PCDs introduced in 
Cevhan and Priebd (|2003l ). ICevhan et all (|2006h . and lCevhan et all (|2007[ l. the relative arc density has tractable 
asymptotic mean and variance. 

We define the underlying graphs of proportional-edge PCDs and their relative edge density in Section [21 
provide the asymptotic distribution of the relative edge density under the null hypothesis in Section 13. li and 
describe the alternatives of segregation and association in Section [221 We prove the consistency of the relative 
edge density in Section 14. 1[ and provide Pitman asymptotic efficiency in Section 14.21 We present the Monte 
Carlo simulation analysis for finite sample performance in Section [SI in particular, provide the Monte Carlo 
power analysis under segregation in Section 15.11 and under association in Section 15.21 We treat the multiple 
triangle case in Section [6l provide extension to higher dimensions in Section [6.41 We provide the discussion and 
conclusions in Section [71 and the tedious calculations and long proofs are deferred to the Appendix. 



2 Relative Edge Density of Underlying Graphs 
2.1 Preliminaries 

The main difference between a graph and a digraph is that edges are directed in digraphs, hence arc called arcs. 
So the arcs arc denoted as ordered pairs while edges arc denoted as unordered pairs. The underlying graph of a 
digraph is the graph obtained by replacing each arc uv Q A or each symmetric arc, {uv, vu} d A hy the edge 
(u, u). The former underlying graph will be referred as the OR-underlying graph, while the latter as the AND- 
underlying graph. That is, the AND-underlying graph for digraph D = (V,-4) is the graph Gand(i^) = (V,£and) 
where £and is the set of edges such that (u, v) S f and uv £ A and vu € A. The OR-underlying graph for 
D = (V,^) is the graph Go^{D) = (V,£or) where £or is the set of edges such that {u,v) E Ed- iS uv G A or 
vu G A. 

The relative edge density of a graph G = (V,£) of order |V| = 7i, denoted p{G), is defined as 

2\£\ 



p{G) 



n{n — 1) 



where | • | denotes the set cardinality function ( Janson et al.l l 2000l )V Thus p(G) represents the ratio of the 



number of edges in the graph G to the number of edges in the complete graph of order n, which is n(n — l)/2. 

Let {n,A4) be a measurable space and consider N : p{^), where p(-) represents the power set functional. 
Then given C fl, the proximity map Ny{-) associates with each point x € J7 a proximity region Ny{x) C fl. 
The Pi-region ri{-,N) : p{^) associates the region ri{x, Ny) := {z £ : x G Ny{z)} with each point 

X £ ft. li Xi, X2, ■ . ■ , Xn are fi-valued random variables, then the Ny{Xi) (and Ti{Xi, Ny)), i = 1,2, . . . ,n arc 
random sets. If the Xi arc independent and identically distributed, then so are the random sets Ny(Xi) (and 

r,iX,, Ny)). 

Consider the data-random PCD D with vertex set V = {Xi,X2, . ■ . , Xn} and arc set A defined by XjXj £ 
A <;=^ Xj £ Ny{Xi). The AND-underlying graph, Gand, of D with the vertex set V and the edge set £and is 
defined by {Xi,Xj) £ £and iS XiXj £ A and XjXi £ A. Likewise, the OR-underlying graph, Gor, of D with the 
vertex set V and the edge set £or is defined by {Xi,Xj) £ £or '^=> XiXj £ A ot XjXi £ A. Then {Xi,Xj) £ 
£and iff Xj £ Ny{X,) and X, £ Ny{Xj) iff Xj £ Ny{X,) and Xj £ Ti{X„ Ny) iff Xj £ Ny{X,) n ri(A:„ Ny). 



Similarly, {Xi,Xj) G £or iff Xj € Ny{Xi) U ri{Xi, Ny). Since the random digraph D depends on the (joint) 
distribution of the Xi and on the map Ny, so do the underlying graphs. The adjective proximity — for the 
catch digraph D and for the map T Vy — - comes from thinking of the regi on Ny(x) as representing those points 
in O "close" to x (jToussaintl (|l980l ) and I Jaromczvk and ToussaintI (|l992l )'l. 



2.2 Relative Edge Density of the AND-Underlying Graphs 



The relative edge density of Gand{D), the AND-underlying graph based on digraph D, is denoted as PundiD). 



For Xi ~ F, Pa,ndiD) is a [/-statistic, 



P.nd{D) = - -tEE 

n[n — 1) ^-^ ^-^ ■' 



where 



Land 
"'ij 



KUXr,x,-N) = i((x„ Xj) e fand) = € A) ■ i{x.,x, e A) 
i[x, e N{x,)) ■ i{Xj e iV(xO) = g n{x^) n Ti (x„ 7V)). 



is the number of symmetric arcs between Xi and Xj in D or number of edges between Xi and in G'and(^)- 
Note that /if"*^ is a symmetric kernel with finite variance since < hand{Xi, Xj \ N) < 1. Moreover, PandiD) is a 
random variable that depends on n, F, and N{-) (i.e., y). But E [pand{D)] only depends on F and N{-). Then 



(1) 



where E [hf^'^] = P{XiX2 G A , X2X1 <E A) ^ P{X2 G N{Xi) D ri{Xi,N)) = /^and(A^) is the symmetric arc 
probability. Note that Pa.nd{N) = P{Xj G 7V(X,) n Ti {X„ N)) for i ^ j. Furthermore, 



< Var [p,nd{D)] 



in — 1)2 



Var 



EE^r 



(2) 



Expanding this expression, we have 

Var [pand(i^)] 



2 -Var [hr^] + Gov [/^^f , 



71 (71 — 1) 



Let A,j be the event that {X^Xj <G A] = {Xj G N{Xi)}, then hf^'^ = I(Ay) • = I(Aij n Aj^). In 

particular, h^'^ = 1(^12) • 1(^21) = I(Ai2 n A21). Then 



Var [h 



and] 



E 



{h'^^) \ - (E [/.-^]) = E [(/.-'i) J - (Ma„d(A^))' . 
Furthermore, E [(/ifS'i)^] = E [(I(Ai2 n ^21))^] = E[(I(Ai2 n A21)] - A'and(iV). So 

Var [/i-'d] = AW(A^) - [A*and(iV)]' ■ 

Moreover, 

Gov [/7?f , /i-'^] = E [hlt.h^"] - E [/i-d] E [h'^^] 
where E [/if^^] ^ g [^a„d] ^ and, 

E ./i^f ] - E[i(Ai2 n A21) (i(Ai3 n A31)] - E[(i(Ai2 n A^i n A13 n A31)] 
= p(X2 G iV(Xi)nri(Xi,iV), X3 g a^(Xi) n ri(Xi,Ar)) 
= p({X2,X3}c7V(Xi)nri(Xi,7V)). 



Thus 



Gov [h\t,h\f] = P{{X2,X^} c Ar(Xi) nri(Xi,iV)) _ [^and(A^)]' 



2.2.1 The Joint Distribution of (/i?" /i?" 

By definition (/ii2"^, is a discrete random variable with four possible values: 

{hif,hr/)emo),{i,o),{QA),{i,i)}. 

Then finding the joint distribution of (/if 2'^, /iia*^) is equivalent to finding the joint probability mass function of 

f Land l,and^ 
(,"12 ; "IS )■ 

First, note that 

{hit, hlf) = (0, 0) iff /iff ^ hlf = iff 1(^12 n A21) = i(Ai3 n ^31) = iff 
i(X2 ^iV(Xi)nri(Xi,7V)) = i(X3 ^iV(Xi)nri(Xi,7V)) = i iff 
1(^2 e r(:y3) \ N{Xi) n ri(Xi, 7V)) = 1(^3 e ^(^3) \ iV(Xi) n ri(Xi, 7V)) = 1 iff 
i{{X2, X3} c T(3^3) \ [NiXi) n ri(Xi, TV)]) = 1. 

Hence Pi{hlf,hlf) = (0,0)) = PHX^.Xs} C r(J^3) \ [iV(Xi) n ri(Xi, TV)]). 

Next, note that {hfi'^, hff) = (1, 1) iff hf^'^ = hff = 1. So P ((/if^'', /if§'^) = (1, 1)) = E [hff.hff] . 

Furthermore, by symmetry P {{hff, hff) = (0, 1)) = P{{hfi'^, hf^'^) = (1, 0)). 

Hence 

pmr,htf)^ (0,1)) = p ((/iff, /iff) = (1,0)) 

= i [1- (P((/iff ,/ifr) = (0,0)) +P ((/iff, /iff) = (1,1)))] . 

2.3 Relative Edge Density of OR-Underlying Graphs 

The relative edge densi 
PoriD) is a t/-statistic, 



The relative edge density of Goi{D), the OR-underlying graph of digraph D, is denoted as Por (D). ForXi - P, 



where 



/i°; = /ior(^^, X,-N) = I((X„ Xj) e £or) = max(I(X,X, e ^), I(X,X, e ^)) 

= I(X, e Ar(x,)uri(x„Ar)). 

is the number of edges between Xi and Xj in Go^{D). Note that /i°j is a symmetric kernel with finite variance 
since < hoi{Xi, Xj-, N) < 1. Moreover, Pq^{D) is a random variable that depends on n, F, and TV(-) (i.e., 3^). 
But E[por(^)] does only depend on F and -/V(-). Then 

< npor{D)] = J_ 5] ^ E[/i-] = E [/ify = P(X2 e N{X,) U ri(Xi, TV)) 

where E [h!^^ = P{XiX2 G .4 V ^2X1 e A) which we denote as Poi{N) for brevity of notation. 
Similar to the AND-underlying case. 



, 7 Var /i?o H ^ ^ 

71 n - 1 ^ n n - 1) 



< Var[por(i^)] = -rVar [/i°^] + -) -( Cov[/i?^, /if^] 



Notice that /ij^ = max(I(Ai2), I(yl2i)) = I(Ai2 U A21). Then 

Var [/i?^] = E [(/ir2)'] - (E K^])' = Mor(iV) - [^^oAN)f 

since E [/i°2] = HofiN) and 

E [(/i?^)2] = E [(I(Ai2 U A2l))2] = E[I(Ai2 U A21)] = PoAN). 



Furthermore, 

where E [h°l] = E[/i°|^] = /ior(iV), and 

E [h°^2K3] = E[(I(^i2 U A21) (1(^13 U ^31)] = E[(I((^i2 U A21) n (Ai3 U ^31))] 

= p{X2 G 7V(Xi)uri(Xi,Ar), X3 e7V(Xi)uri(Xi,7V)) 
= p({X2,X3}ciV(Xi)uri(Xi,7V). 

So 

Cov[/i-, = P({X2,X3} C N{X,) U ri(Xi, iV)) - [Mor(iV)]2. 
Remark 2.1. Note that % = /if"^ + since ) + = I(Ay n A^i) + I(Aij U A^i). □ 



2.3.1 The Joint Distribution of (/i^'f , /i^"! ) 

Finding the joint distribution of {h°2, /iis) is equivalent to finding the joint probabihty mass function of {h°2, h™^), 
i.e., finding Piih°l,h1l) = for each e {(0, 0), (1, 0), (0, 1), (1, 1)}. 

First, note that 

ihT2,Kl) = (0,0) iff = h°l = iff 1(^12 U A21) = I(Ai3 U A31) = iff 
MX2 <^ N{Xi)uri{Xi,N)) = I{X3 <^ N{Xi)uri{Xi,N)) = 1 iff 

i{{X2,X3} c r(3^3) \ [NiXi) u ri(Xi, AT)]) = i. 

Hence P{{h'i'2, hfa) = (0,0)) = P{{X2,X3} C r(3;3) \ [N{Xi) UTiiXi, N)]). 

Next, note that ihr2,h'il) = (1, 1) iff h°,^2 = /i?3 = 1- P{{hT2,hT3) = (1, 1)) = E K^./i?^]. 

By symmetry P((/i?^, /i?^) = (0,1)) = P{{hT2,h°,l) = (1,0)). Hence 

P{{hT2,h1l) = {0,l)) = P{{hT2,h°l)^{l,0)) 

= ^ (1 - [PiihT2, hT^) = (0, 0)) + PiihT2,h1l) = (1, 1))]) . 



2.4 Proportional-Edge Proximity Maps and the Associated Regions 

Let 51 = and 3^3 = {yi,y2,y3} C be three non-coUinear points. Denote by r(3^3) the triangle (including 
the interior) formed by these three points. For r G [l,oo] define Np^{x) to be the proportional- edge proximity 
map with parameter r and Ti{x) := Fi (x, Np^) to be the corresponding Fi-rcgion as follows; see also Figurcs[T] 
andm Let "vertex regions" i?(yi), R{y2), Rivs) partition T(3^3) using segments from the center of mass of r(3^3) 
to the edge midpoints. For x £ ^(3^3) \ 3^3, let v{x) £ 3^3 be the vertex whose region contains x; x £ R{v{x)). If 
X falls on the boundary of two vertex regions, or at the center of mass, we assign v{x) arbitrarily. Let e{x) be the 
edge of r(3^3) opposite v{x). Let £{v{x),x) be the line parallel to e(x) through x. Let d{v{x),£{v{x),x)) be the 
Euclidean (perpendicular) distance from v{x) to £{v{x),x). For r £ [l,oo) let £r{v{x),x) be the line parallel to 
e(x) such that (i(f (x), -^^(^(0;), a;)) = rd{v{x),£{v{x),x)) and d{£{v{x),x),£r{v{x),x)) < d(v{x),£r(v{x),x)). Let 
Tj.{x) be the triangle similar to and with the same orientation as r(3'3) having v(x) as a vertex and £r{v{x),x) 
as the opposite edge. Then the proportional-edge proximity region Np^{x) is defined to be Tr{x) n r(3^3). 

Furthermore, let £,i{x) be the line such that ^!;(a;) nT(3^3) ^ and r d(yi, ^i{x)) = d{yi,£{yi, x)) for i = 1, 2, 3. 
ThenrUx)ni?(y,) = {z £ R{y,) : d(y„£(y„z)) > d(y„e,(x)}, for z = 1,2,3. Hence FU^) = ULi(rUa^)nP(y.)). 
Notice that r > 1 imphes x £ Np^{x) and x £ Tl{x). Furthermore, limr^cx^ Np^{x) = T(3'3) for all x £ 
Tiys) \ 3^3, and so we define N^{x) = TO^s) for all such x. For x £ 3^3, we define N^Eix) = {x} for all 
r £ [1,00]. Then, for x £ i?(y,;) \imr-too^i{x) = ^'(3^3) \ {yj,yfe} for distinct and k. 

Notice that Xi F, with the additional assumption that the non-degenerate two-dimensional probability 
density function / exists with support in T'(3'3), implies that the special cases in the construction of Np^ — X 
falls on the boundary of two vertex regions, or at the center of mass, or X g 3^3 — occur with probability zero. 
Note that for such an F, Ny{x) is a triangle a.s. and ri(x) is a convex or nonconvex polygon. 



Figure 1: Construction of proportional-cdgc proximity region, Ny ^(x) (shaded region) for an x e i?(y)l). 



2.5 Relative Edge Density of the Underlying Graphs of Proportional-Edge PCDs 

Consider the underlying graphs of the data-random PCD D with vertex set V = {^i, X2, ■ ■ ■ , Xn} and arc set 
A defined by (X^, Xj) e A Xj e Ny{Xi). Recall that (X„ Xj) e fand iff X^ e N'pj^iXi) n Ti (X^, A^^^) 

and {X,,X,) e for iff X, G A^^bI^,) U Ti (X,, TV^^). 

Let h,^4X,,Xf,N^pE) = I{X, € A^^s(X,) n (X,)) and /i°;(r) := hoAX^, Xf, N^pj,) = l{Xj G 

Np^{Xi) U (Xi)) for i 7^ The random variable /0^"'^(r) := Pand('^n; h,Np^) depends on n explicitly, and 
on F and Npp, implicitly. The expectation E , however, is independent of n and depends on only F 

and N-pE. Let /.and(r) := E [hlf{r)] and i.a„d(r) Gov [hlf{r), hlf{r)] . Then 




Figure 2: Construction of the Fi-region, ri~'^{x) (shaded region) for an a: G i?(y)l). 




(3) 




Gov [hlf{r),hlf{r)] <l/4. 



(4) 



A central limit theorem for J7-statistics (iLehmann (llQQgh l yields 



- /iand(r)) AA(0,4i.a„d(r)) 



(5) 



provided i'and('') > 0. The asymptotic variance of p!^'^{r), 4i'and(''), depends on only F and Np^. Thus we 
need determine only /^and('') and Vandi^) in order to obtain the normal approximation 



The above paragraph holds for = po^{Xn\ h, Np-^) also with is replaced by Pn{f); h^2^{r) and 

h^'^{r) are replaced by ^'^d /i°3, respectively. 

For r = 1, Afp^-'^(.T)nr5^=-'^(x) = i'(w(a;), x) which has zero R'^-Lcbesgue measure. Then we have E = 1)] 

E [hf:f{r = 1)] = Ma„d(r = 1) = P{X2 e N^p=e\Xi) n rr'(^i)) = O. Similarly, P{{X2,X^} C iV^lH^i) n 
r^=i(Xi)) = 0. Thus, zyand(r = 1) = 0. Furthermore, for r = oo, A^^,i°°(x) n ri=°°(x) = r(3^3) for aU 
X e T(J^3)\3^3- Then E [p^^'^ir = oo)] E [hlf{r = c^)] = /iand(r = c^) = ^(^2 G iV^^°°(Xi)nrr°°(^i) = 
P{X2 e TCJ^a)) = 1. Similarly, P{{X2,X3} C A^^i;°°(Xi) n r™(Xi)) = 1. Hence i.and(r = ^) = 0. 
Therefore, the CLT result in Equation ^ holds only for r e (l,oo). Furthermore, = 1) = a.s. and 

PT'^i'f = oo) = 1 a.s. 

For r = 1, Npl^^{x) U T\=^{x) has positive R^.Lebesgue measure. Then P{{X2,X3} C Npl^^{Xi) U 
r^=i(Xi)) > 0. Thus, i/or(r = 1) 7^ 0. On the other hand, for r = oo, N^=^°°{Xi) U r^=°°(Xi)) = T(3^3) 
for all Xi € r(3;3). Then E [p°'{r = oo)] = E [/i°^(r = oo)] = P(X2 G iV^^°°(Xi) ur^=°°(Xi)) = /Xor(r = oo) = 
P(X2 G T(3^3)) = 1. Similarly, P({X2,X3} C A^^^°°(Xi) U r^=°°(Xi)) = 1. Hence j/or(?- = oo) = 0. Therefore, 
the CLT result for the OR- underlying case holds only for r G [1, oo). Moreover p™{r ~ oo) ~ 1 a.s. 

Remark 2.2. Relative Arc Density of PCDs: The relative arc density of the digraph D is denoted as p{D). 

For Xi ~ F, p{D) is also shown to be a f/-statistic ( Cevhan et al.l ( 20061) ). 



where hjj = h(Xi, Xj; N) = I(XiXj G .4) = I(Xj G N(Xi))) is the number of arcs between Xi and Xj in D. 
Here 

< E [piD)] = ^ J2 E EN] = E [/ii2] /2. 
nln — 1 ^ — ' — ' 



Furthermore, Moreover, 



Gov [/ii2, /lis] - P({X2, X3} C iV(Xi)) - [E [/112]]' . 



Let hij{r) := h{Xi, Xj;Np^) = I{Xj G Np^{Xi)) for i ^ j and the random variable p„(r) := fe , iVpg ). 

Let ^(r) := E [pn(r)] and v{r) := Gov [/112(f), /ii3(7')]- A central limit theorem for C-statistics ( Lehmann ( 1999t) ) 
yields 

^^(p„(r)-/i(r))-^AA(0,KO) (7) 



provided i^(r) > 0. The explicit forms of asymptotic mean p{r) and variance v{r) are provided in Cevhan et al.l 
()2006l) . □ 



3 Relative Edge Density under Null and Alternative Patterns 
3.1 Null Distribution of Relative Edge Density 

The null hypothesis is generally some form of complete spatial randomness; thus we consider 

Ho-.x^'^^uiTiys)). 

If it is desired to have the sample size be a random variable, we may consider a spatial Poisson point process 
on r(3^3) as our null hypothesis. 

We first present a "geometry invariance" result which will simplify our subsequent analysis by allowing us 
to consider the special case of the equilateral triangle. Let p^^'^ir) := panA{n;U{T{y3)),Np^) and p^'{r) 

PoM-M{T{y^)),N-p„). 



Theorem 3.1. Geometry Invariance: Let = {yi,y2,y3} C be three non-collinear points. For i = 
1, 2, . . . , n let Xi^'^ F ~ U{T{y3)), the uniform distribution on the triangle T^y^). Then for any r G [1, oo] the 
distribution of p'^"'^{r) and p'^{r) is independent ofy^, and hence the geometry ofT{y^). 

Proof: A composition of translation, rotation, reflections, and scaling will take any given triangle To = 
r(yi,y2,y3) to the "basic" triangle Tf, = T((0,0), (1,0), (ci,C2)) with < ci < i, C2 > and (1 - cif + cl < 1, 

preserving uniformity. The transformation i/i : ^ given by (j){u, u) = (^u + ^^^^ v, takes Tb to the 

equilateral triangle Te = T ((0,0), (1,0), (1/2, \/3/2)). Investigation of the Jacobian shows that also preserves 
uniformity. Furthermore, the composition of (f> with the rigid motion transformations and scaling maps the 
boundary of the original triangle To to the boundary of the equilateral triangle Tg, the median lines of To to 
the median lines of Te, and lines parallel to the edges of To to lines parallel to the edges of T^. Since the 
joint distribution of any collection of the /if"'^(r) and h°j{r) involves only probability content of unions and 
intersections of regions bounded by precisely such lines, and the probability content of such regions is preserved 
since uniformity is preserved, the desired result follows. ■ 

Based on Theorem 13.11 for our proportional-edge proximity map and the uniform null hypothesis, we may 
assume that T(3^3) is a standard equilateral triangle with 3^3 = {(0,0), (1,0), (1/2, -\/3/2)} henceforth. 

In the case of this (proportional-edge proximity map, uniform null hypothesis) pair, the asymptotic null 
distribution of p'^^^ir) and /3™'(r) as a function of r can be derived. Recall that ^and('') = E [/ii2'^(?')] = P{X2 G 
N^pe{Xi) n r^Xi)) = /za„d(r) and /ior(r) = E [/i?^] = P(X2 S ^^^(Xi) U r^Xi)) = fi^rir) are the probability 
of an edge occurring between any two vertices in the AND- and OR-underlying graphs, respectively. 

Theorem 3.2. Asymptotic Normality: For r G (l,oo), 



and for r G [1, 00), 



where 



/7i (pr'(^) - Mand(O) / v/4z^a„d(r) AA(0, 1) 



Mand('') 



_J_ (-l+r)(5 r^-148 -232 r-|-128) 

54 r'-'{r+2)(r+l) 
]_ 101 r^-l-1302 1-^-732 r^-536 r+672 

216 r{r+2){r+l) 
1 r^-13 r^+30 r° + 148 r^-448 r'^+264 r^-|-288 r^-368 r+96 



(r-''+3r^-2+2 r)(-l+r)^ 
r*(r+l) 



r'^(r+2)(r+l) 



for rG [1,4/3), 

for rG [4/3,3/2), 

for rG [3/2,2), 

for r G [2, oo). 



(8) 



( 47 r"- 195 -1-860 r'*-846 7 



-108r^+720r-256 



108r2{V-|-2)(r-|-l) 
175 r^-579r'' + 1450r^-732 r^-536 r+672 
216r(r-|-2)(r+l) 
3 r^-7 r'^-30 r'^-|-84 r-^-264 r''+304 r^ + 144 r^-368 r-t-96 
Wr^r+T)(r+2) 

r--'+r'^-&r+2 
r*(r-|-l) 



for rG [1,4/3), 

for rG [4/3,3/2), 

for rG [3/2,2), 

for r G [2, oo). 



(9) 



11 



Va,nd 



i=i 
11 



(10) 



(11) 



where df^'^{r) and (r) are provided in Appendix Sections 1 and 2, and the derivations of jiandir) and i'and('') 
are provided in Appendix 3, while those of Poi-if) and Vorir) are provided in Appendix 4. 

Notice that //and('' = 1) = and linv-,oo /^and(?') — 1 (at rate 0(r~^)); and poi{r = 1) = 37/108 and 
liiHr^oo Porir) = 1 (at rate 0(r^^)). 

To illustrate the limiting distribution, for example, r = 2 yields 



V^(Prj£^_Mand(2)) 
\/4l/a„d(2) 



362880 n 
58901 



and 
rn 



(2)-^) ^AA(0,1) 



Figure 3: Result of Theorem 13.21 asymptotic null means /i(r), /iand(''), and ^ot{r) (left) and variances 
4i/and(?'), and Avoi{r) (right), from Equations ([8]), ([9|), and (fTO|) . pT|) . respectively. Some values of note: /^(l) = 
37/216, /^and(l) = 0, and ^or(l) = 37/108, lim^^cx) m(^) = linv^oo Aiand(r) = linv^oo Mor(r) = 1, i'and(r = 1) = 
and linv^oo zyand('') = 0, Vc„:{r 1) = 1/3240 and linir^oo '^or('') = 0, and argsup,,gj]^ z^(r) « 2.045 with 
snp,rg[i,oo] '^(^) ~ -1305, argsup^g[i 4i/and(?') ~ 2.69 with sup^gj^ ^oj i/and(?') ~ .0537, argsup^gf^ j^or(?') « 
1.765 with sup^gjj t'or(^) ~ .0318. 



and 



or cquivalently. 



V^(p-(2) - /io,.(2)) _ / 120960 n / „ 19 



P"(2)-^) ^AA(0,1) 



^4:.or(2) V 13189 V 24 

a„d/2\ ^P^™'' A/- (^ii ^^^^^ \ and «°'-('S>^ ''PS.'™ KT 
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By construction of the underlying graphs, there is a natural ordering of the means of relative arc and edge 
densities. 

Lemma 3.3. The means of the relative edge densities and arc density have the following ordering: fiand{f) < 
fi{r) < ^or{f) for all r e [l,cx)). Furthermore, for r ~ oo we have ^landir) ~ tJ-{r) = fJ-orir) = 1- 

Proof: Recall that /ia„d(r) = EipT'^ir)] = -P(^2 G A^^^ (Xi ) n r^Xi ) ) , fi{r) = E[p„(r)] = P{X2 € iV^^l^i))- 
and//or(r) = EiPnir)] = P{X2 G N-pE{Xi)UTl{Xi)). And iV^^(Xi)nrUXi) C iV^s(^i) ^ ^Pij(^i)urU^i) 
with probability 1 for all r > 1 with equality holding for r = oo only. Then the desired result follows. See also 
Figure [1 ■ 

Note that the above lemma holds for all Xi that has a continuous distribution on T(3^3). There is also a 
stochastic ordering for the relative edge and arc densities as follows. 

Theorem 3.4. For sufficiently small r, p'!^'^{r) <^'^ Pn{r) <^'^ PrT(^) as n ^ oo. 

Proof: Above we have proved that /Xand('') < p{r) < tJ-oi{r) for all r € [l,oo). For small r {r < r fa 1.8) 
the asymptotic variances have the same ordering, 4i^and(^) < '^{f) < Avorir)- Since pf^^{r), p„(r), p°/(r) are 
asymptotically normal, then the desired result follows. See also Figure [31 ■ 

Figures m and [5] indicate that, for r = 2, the normal approximation is accurate even for small n although 
kurtosis may be indicated for n = 10 in the AND-undcrlying case, and skewness may be indicated for n = 10 
in the OR-underlying case. Figures [6] and [7] demonstrate, however, that severe skewness obtains for some 
values of n, r. The finite sample variance and skewness may be derived analytically in much the same way 
as was 4fand(?') (and 4i'or('')]) for the asymptotic variance. In fact, the exact distribution of pf^'^ir) (and 
Pri {''')) i'^ principle, available by successively conditioning on the values of the Xi. Alas, while the joint 
distribution of h'l^ {^),hf^'^{r) (and h°l{r) , h™^{r)) is available, the joint distribution of {/if"''(r)}i<i<j<„ (and 
{/i°j'(r)}i<i<j<„), and hence the calculation for the exact distribution of Pn^'^{r) (and p™(r)), is extraordinarily 
tedious and lengthy for even small values of n. 

Let jnir) be the domination number of the proportional-edge PCD based on Xn which is a random sample 
fromZ/^(T(3^3)). Additionally, let 7i^"'^(r) and 7"(r) be the domination number of the AND- and OR-underlying 




Figure 4: Depicted are p^'"^{2) M (li, 3^^) for n ^ 10, 20, 100 (left to right). Histograms are based 

on 1000 Monte Carlo replicates. Solid lines arc the corresponding normal densities. Notice that the vertical 
axes are differently scaled. 




Figure 5: Depicted are p°'^(2) Af (if, j^^omon ) for n = 10, 20, 100 (left to right). Histograms are based on 
1000 Monte Carlo replicates. Solid lines are the corresponding normal densities. Notice that the vertical axes 
are differently scaled. 



graphs based on the proportional-edge PCD, respectively. Then we have the following stochastic ordering for 
the domination numbers. 

Theorem 3.5. For all r e [l,oo) and n > 1, 7°'Xr) <^'^ 7„(r) <-5^ 7n""^{r). 

Proof: For all x € T(3^3), we have N^^^ix) n F^ (x) C N],e{x) C N^e{x) U F^ (x). For X - U{T{y3)), 
we have N^si^) ^ F^ (X) C N],j^{X) C N^e{X) U Tl (X) a.s. Moreover, 7„(r) = 1 iff Af^ C iVp£;(X,) for 
some i] -/T'^ir) = 1 iff Xn C N^eI^^) ^1 i^^) for some i; and 7°'(r) = 1 iff X„ C N^si^^) U ^[ (X^) for 
some i. So it follows that ^'(7™'^(r) = 1) < P(7„(r) = 1) < Pi'Jn ir) = 1). In a similar fashion, we have 
PhT'^ir) < 2) < P(7„(r) < 2) < P(7°'(r) < 2). Since P(7„(r-) < 3) = 1 (ICevhan and Priebd (l2005l) l. it 



follows that P(7™(?') < 3) = 1 also holds as P(7„(r) < 3) < P(7™(?') < 3). Hence the desired stochastic 
ordering follows. H 

Note the stochastic ordering in the above theorem holds for any continuous distribution F with support 
being in T{y3). 



3.2 Alternatives: Segregation and Association 

The phenomenon known as segregation involves observations from different classes having a tendency to repel 
each other — in our case, this means the Xi tend to fall away from all elements of 3^3. Association involves 
observations from different class es having a te nden cy to attract one ano ther, so that the Xi tend to fall near an 
element of 3^3. See, for instance, Dixon ( 1994 ) and Coomes et al. ( 1999[ ). 



Figure 6: Depicted are the histograms for 10000 Monte Carlo replicates of pIq'^{1.05) (left) and Pio'^{5) (right) 
indicating severe small sample skewness for extreme values of r. Notice that the vertical axes are differently 
scaled. 



Figure 7: Depicted are the histograms for 10000 Monte Carlo replicates of Pio(l) (left) and /Oio(5) (right) 
indicating severe small sample skewness for extreme values of r. Notice that the vertical axes are differently 
scaled. 

We define two simple classes of alternatives, and Hf with e £ (O, V3/3), for segregation and association, 
respectively. For y S 3^3, let e(y) denote the edge of r(3^3) opposite vertex y, and for x £ r(3^3) let £y{x) denote 
the line parallel to e(y) through x. Then define T(y,e) = {x £ T(3^3) : d{)/,£y{x)) < e}. Let be the model 

under which X, - U{T{y3) \ Uyey,T(y, e)) and be the model under which X, - Z^(Uyey,T(y, ^/3/3 - e)). 
Thus the segregation model excludes the possibility of any Xi occurring near a yj, and the association model 
requires that all X; occur near a yj. The V3/3 — e in the definition of the association alternative is so that e = 
yields Ho under both classes of alternatives. 

Remark 3.6. These definitions of the alternatives are given for the standard equilateral triangle. The geometry 
invariance result of Theorem 13. II still holds under the alternatives and H^. In particular, the segregation 
alternative with e £ (O, \/3/4) in the standard equilateral triangle corresponds to the case that in an arbitrary 
triangle, S x 100% of the area is carved away as forbidden from the vertices using line segments parallel to the 
opposite edge where S = (which implies 5 £ (0, 3/4)). But the segregation alternative with e £ (a/3/4, V3/3) 
in the standard equilateral triangle corresponds to the case that in an arbitrary triangle, S x 100% of the 
area is carved away as forbidden around the vertices using line segments parallel to the opposite edge where 
(5 = 1 — 4 (1 — \/3e) (which implies S £ (3/4, 1)). This argument is for the segregation alternative; a similar 
construction is available for the association alternative. □ 

The asymptotic normality of the relative edge density under the alternatives follows as in the null case. 

Theorem 3.7. Asymptotic Normality under the Alternatives: Let iJ,and{r,s) be the mean and z^and('', e) 
be the variance of Pn^'^{r) under the alternatives for r £ [l,oo) and e £ (O, V3/3). Then under and , 



-y/n(p™''(r) — ^and{i', s)) — > A/'(0, 4 lyandir, e)) for the values of the pair (r, e) for which Vandif, e) > 0. A similar 
result holds for p°i(r). 

Proof: Under the alternatives, i.e., e > , Pn^'^{r) is a [/-statistic with the same symmetrie kernel hfj'^{r) as 
in the null case. Let Ee[-] be the expectation with respect to the uniform distribution under the alternatives 
withe G (0, V3/3). The mean Pundir, e) = Eg [/9™'^(r)] = [/ii2'^(?')] , now a function of both r and e, is again 
in [0,1]. The asymptotic variance, 4i^and('', s) = 4 Gov [hf2'^{r),hl^'^{r)'^, also a function of both r and e, is 
bounded above by 1/4, as before. Thus asymptotic normality obtains provided Vandir, s) > 0; otherwise p^^{r) is 
degenerate. Then under iff, i/a„d(r-, e) > for (r, e) in (1, VS/{2e)) x (0, \/3/4] or (1, Vs/e - 2) x (V3/4, 73/3), 
and under Hf, i/and('', e) > for (r, e) in (l,oo) x (O, ■\/3/3). Also under ijf, :^or(r, e) > for (r, e) in 
[l,y3/(2£))x(0, V3/4] or [1, V3/£-2) x (^3/4, ^3/3), and under ii^^ , i^orir,e) > for (r,£) in (1, oo) x (O, V3/3) 
or {1} x (0,^/3/12). ■ 

Notice that for the association class of alternatives any r G (l,oo) yields asymptotic normality for all 
e e (0, V3/3) in both AND- and OR-undeiiying cases, while for the segregation class of alternatives only r ~ I 
yields this universal asymptotic normality in the OR-underlying case, and such an e does not exist for the 
AND-underlying case. 

The relative edge density of the underlying graphs based on the PCD is a test statistic for the segrega- 
tion/association alternative; rejecting for extreme values of pf^'^{r) is appropriate since under segregation we 
expect p^"^(r) to be large, while under association we expect pf^'^{r) to be small. The same holds for p°^{r). 
Using the test statistics 

RT\r) = (p-'i(r) - Aiand(r)) /V4^^and(r), and i?-(r) = ^ (p°:(r) - Mor(r)) / ^/l^) (12) 

for AND- and OR-underlying cases, respectively, the asymptotic critical value for the one-sided level a test 
against segregation is given by 

z„ = $-i(l-a) (13) 

where $(•) is the standard normal distribution function. The test rejects for i?™'^(r) > Zq, against segregation. 
Against association, the test rejects for R^'^{r) < zi_q. The same holds for the test statistic R™{r). 



4 Asymptotic Performance of Relative Edge Density 
4.1 Consistency 

Theorem 4.1. The test against which rejects for R'^'^{r) > and the test against which rejects for 
Rf^'^ir) < o,re consistent for r G (1, oo) and e G (O, \/3/3) . The same holds for R^{r) with r G [1, oo). 

Proof: Since the variance of the asymptotically normal test statistic, under both the null and the alternatives, 
converges to as n ^ oo (or might be zero for n < oo), it remains to show that the mean under the null, 
A''and('^) = E[p™'^(r)], is less than (greater than) the mean under the alternative, /^and(^:e) = E [/^^^"^(r)] 
against segregation (association) for e > 0. Whence it will follow that power converges to 1 as n ^ oo. Let Pe(-) 
be the probability with respect to the uniform distribution under the alternatives with e G (0,\/3/3). Then 
against segregation, we have 

/iand(r) = Po(X2 G iVp£;(Xi) n r^(Xi)) 

= Po(X2 G iVpB(Xi) n rl-(Xi),Xi G VJyey^Tiy.e)) + Po{X2 G iVpB(Xi) n V{{X^),X^ G Tiy^) \ Uyey,T{y,e)) 
= Po{X2 G N'pEiXi)nrl{Xi)\X-, e Uyey,T(y,e))PoiXi G Uyey,T{y,e)) 

+Po{X2 G N^pe{Xi)\X^ G T{ys)\Uyey,T{y,s))Po{Xj G T{ys) \Uyey,T{y,e)) 
< Po{X2 G N^PE{Xi)nrl{Xi)\Xi G Uyey,T{y,e))pi + P,{X2 G iV^B(Xi) n rI(Xi)iXi G T{y3)\Uyey,T{y,e))p2 
= E„(/(X2 G 7V^£(Xi)nr5:(Xi))|Xi G Uyey,T{y,e))pi 

+Ee{I{X2 G N^pe{Xi) n rI(Xi))|Xi G T{ys) \ Uyey,T{y, e))p2 
where pi = Po(^i e ^yey,T{y, e)) and p2 = Po{Xi G ^(^3) \ ^yey,T{y, e)) = l-pi. Then 

Pand{r,e) > Pa.nd{r) — — = Pand{r)- 

P2 

Likewise, we have ^and(?', e) = E^ [/9f"'^(r)] < E[p„(r)] = Pandir), for association. 
The consistency follows for the OR-underlying case in a similar fashion. ■ 



4.2 Pitman Asymptotic Efficiency 



Pitman asymptotic efficiency (PAE) provides an investigation of "local asymptotic power" — local about iJo. 
T his involves the limit as e — > as well a s the limit as n ^ oo. A detailed discussion of PAE can be found 
in Kendall and Stuart! (I1979I) and lEedenl ^96^. For segregation or association alternatives with the AND- 

((Mand)('=H^£ = 0))' 



underlying graphs, the PAE is given by 



(r) 



where (/i. 



0) is the fc derivative with 



respect to e so that (Aiand)^'^H''j ^ = 0) 7^ but (^and)^'^ ^H'I'tS = 0) = for fc = 1, 2, . . .. Likewise the same 
holds for the OR- underlying case. Then under segregation alternative , the PAE is given by 



((A^and)^^(r,g^O))^ 
I^and(?') 



and PAEf,(r) = 



((Mo.r(r,£ = 0))^ 



since (/Xand)'(?', £ = 0) = and (/Xor)'(r, e = 0) = 0. Under association alternative is 



((Aw)"(r,£ = 0))^ 



and PAE;7^(r) 



((Mo.r(r,£^0)f 



since (/iand)'(T', e = 0) = (/ior)'(?', e = 0) = 0. Equations (fTO|l and (fTTjl provide the denominators; the numerators 
require a bit of additional work, but /Ltand(r, e) and fior{''',s) are available for small enough e, which is all we 
need here. See Appendix 5 for explicit forms of /^and ('',£) and /^or ('',£) for segregation and association, and the 
derivations of Hiind{r,e) and Horirje) are provided in Appendix 6. 

Let PAE'^(r) and PAE"^(r) denote the PAE score against the segreg ation and ass o ciatio n alternatives, 
respectively, for the relative arc density of the PCD based on Np^ (see Cevhan et al. ( 2006() more detail). 
Figure [5] presents the PAE as a function of r for both segregation and association in the digraph, AND, 
and OR-underlying graph cases. For large n and small e, PAE analysis suggests choosing r large for test- 
ing against segregation in all three cases and choosing r small for testing against association, arbitrarily close 
to 1 for the AND- and OR-underlying cases, but around 1.1 for the digraph case. Furthermore, in segrega- 
tion, PAEf,.(r) < PAE'^(r) < PAEfjj(j(r), suggesting the use of AND-underlying version. Under association, 
max (PAEf,j(j(r), PAE'^(r) < PAEfj.(r)) implying the use of OR-underlying version. 




' 1 1.5 2 2.5 3 3.5 4 




Figure 8: Pitman asymptotic efficiency against segregation (left) and association (right) as a function of r. Some 
values of note: PAE^(r = 1) = 160/7, PAE^^^{r = 1) = 4000/17, PAEf,(r = 1) = 160/9, lim^^oo PAE^(r) = 
lim,.^oo PAE^„d(r) = lim^^oo PAE£(r) = 00, and PAEf„d(?') has a local supremum at « 1.35. Also PAE^(r = 
1) = 0, PAE£d(r = 1) = PAE^,(r = 1) = 00, lim.^oo PAE^{r) = lim.^oo PA-^tdir) = lim.^oo PAE^^ir) = 0, 
argsup^g[]^ PAE'^(r) w 1.1, and PAE^,j^(r) has a local supremum at r = 1.5 and a local infimum at r w 1.2 

Remark 4.2. Hodges-Lehm ann Asymptotic Efficiency: Hodges-Lehmann asymptotic efficiency (HLAE) 



(jHodges and LehmannI ()1956l )) is given by 

HLAE(prW,£) 



and/ ^ .^ (Mand(?',e) -Ma„d('')) 



Unlike PAE, HLAE does not involve the limit as e ^ 0. Since this requires the mean and, especially, the 
asymptotic variance of p^'^{r) under the alternative, we avoid the explicit investigation of HLAE. HLAE for 
OR-underlying graphs can be defined similarly. The ordering of HLAE seems to be the same as that of PAE. □ 



Rema rk 4.3. Asymptotic Pov^rer Function Analysis: The asymptotic power function (jKendall and Stuart 
allows investigation of power as a function of r, n, and e using the asymptotic critical value and an appeal 
to normality. Under a specific segregation alternative , the asymptotic power function for AND- underlying 
graphs is given by 

nf„d(^,«,£) = 1-$ 

Under H^, we have 



v/twW" , (Mand (r) - Mand (''j e) ) \ 

\ VI^and(r,e) VI^and(r,£) / 

For OR-underlying graphs, the asymptotic power functions, Il^j.{r,n,e) and H^j.(r, n, e), are defined similarly. 
However it is not investigated in this article. □ 



5 Monte Carlo Simulation Analysis for Finite Sample Performance 

We implement the Monte Carlo simulations under the above described null and alternatives for r £ {1, 11/10, 6/5, 4/3, v^, 
3/2, 2, 3, 5}. 

5.1 Monte Carlo Power Analysis under Segregation 

In Figure [HI we present a Monte Carlo investigation against the segregation alternative for r = 1.1 and 

V 3/8 

n = 10 (left) and n = 100 (right). The empirical power estimates are calculated based on the Monte Carlo 
critical values. Let (3^^^ (pf"'^(?')) and /3,fj^ {Pn {''")) stand for the corresponding empirical power estimates for the 
AND- and OR-underlying cases. With n = 10, the null and alternative probability density functions for Piq'^{1.1) 
and p°q(1.1) are very similar, implying small power (10,000 Monte Carlo replicates yield empirical power values 
Pic {Pio'^) = 0.1318 and {pf„) = 0.0539). Among the 10000 Monte Carlo rephcates under Ho, we find the 
95*'' percentile value and use it as the Monte Carlo critical value at .05 level for the segregation alternative, and 
use 5*'' percentile value for the association alternative. With n = 100, there is more separation between null 
and alternative probability density functions in the underlying cases where separation is much less emphasized 
in the OR-underlying case; 1000 Monte Carlo rephcates yield P^^^ (pfoo) = 0-994 and (pll^^) = 0.298 where 
the empirical power estimates are based on Monte Carlo critical values. Notice also that the probability density 
functions arc skewed right for n ~ 10 in both underlying cases, while approximate normality holds for n = 100. 

For a given alternative and sample size we may consider optimizing the empirical power of the test as a 
function of the proximity factor r. Figure [TU] presents a Monte Carlo investigation of empirical power based on 
Monte Carlo critical values against H^^^ and H^^^ as a function of r for n = 10 with 1000 replicates. The 
corresponding empirical power estimates are given in Table [TJ Our Monte Carlo estimates of r* , the value of 
r which maximizes the power against , are ''^^g = 3 and G [4/3,3] in the AND-underlying case, and 

'^^3/8 ~ ^ ^"^^ ^\/3/4 ^ t'^/'^' ^] OR-underlying case. That is, more severe segregation (larger e) suggests 

a smaller choice of r in both cases. For both e values, smaller r values are suggested in the OR-underlying case 
compared to the AND-underlying case. 

For a given alternative and sample size we may consider analyzing the power of the test — using the 
asymptotic critical value — as a function of the proximity factor r. Let an(r) denote the empirical significance 
levels and I3n{r) empirical power estimates based on the asymptotic critical value. Figure [11] presents a Monte 
Carlo investigation of empirical power based on asymptotic critical value against H^^^ and H^^^ as a function 
of r for n = 10. The corresponding empirical power estimates are given in Table O In the AND-underlying 
case, the empirical significance level, a„=io(?'), is closest to .05 for r = 2 and 3 which have the empirical power 
/3io(2) = .3846 and /3io(3) = .5767 for e = ^/3/8, and /1io(2) = /3io(3) = 1 for e = ^/3/4. In the OR-underlying 
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Figure 9: Two Monte Carlo experiments against the segregation alternatives H^r- Depicted are kernel density 

V 3/8 

estimates of for n = 10 (top left) and n = 100 (top right) and p°''(l.l) for n = 10 (bottom left) and 

n = 100 (bottom right) under the null (solid) and alternative (dashed) cases. 

case, the empirical significance level, an=io{r), is closest to .05 for r = 2 — larger for all r values — which have 
the empirical power /?io(2) ~ .1594 for e ~ a/3/8, and /3io(2) = 1 for e = V3/4. So, for small sample sizes, 
moderate values of r is more appropriate for normal approximation, as they yield the desired significance level, 
and the more severe the segregation, higher the power estimate. Furthermore, the AND-underlying version 
seems to perform better than the OR-underlying version for segregation alternatives. 

5.2 Monte Carlo Power Analysis under Association 

In Figure [T^ we present a Monte Carlo investigation against the association alternative H^^^^ for r = 1.1 and 
n = 10 (left) and n = 100 (right). The empirical power estimates are calculated based on the Monte Carlo 
critical values Let (Pn"'^(^)) and P:^^ (Pn ('')) stand for the corresponding empirical power estimates for the 
AND- and OR-underlying cases. As above, with n ~ 10, the null and alternative probability density functions for 
p^o'^(l.l) and p°o(l-l) are very similar, implying small power — in fact, virtually no power — (10,000 Monte Carlo 
replicates yield the following empirical power estimates based on Monte Carlo critical values: (pfo'^) = 0.0 

and P:^^ {pfo) = 0.0). With n = 100, there is more separation between null and alternative probability density 
functions in the underlying cases where separation is much less emphasized in the AND-underlying case; for this 
case, 1000 Monte Carlo replicates ^ield the following empirical power estimates based on Monte Carlo critical 
values: P:^^^ (Pioo) = 0.009 and P:^^ (Pioo) 0.939. Notice also that the probability density functions are 
skewed right for n = 10 in both underlying cases, with more skewness in OR-underlying case, while approximate 



r r 

Figure 10: Empirical power estimates based on Monte Carlo critical values as a function of r against segregation 
alternatives with the AND-underlying case (top two) and OR- underlying case (bottom two); in both cases, we 
have -ff^/g (left) and H^^^ (right) for n = 10 and N,nc = 1000 Monte Carlo replicates. 

normality holds for n = 100 for both cases. 

In Figure [T31 we also present a Monte Carlo investigation of empirical power based on Monte Carlo critical 
values against H^^_^^ and ^^^^24 ^ function of r for n = 10 with 1000 replicates. The corresponding 
empirical power estimates are presented in Tabic [3] Our Monte Carlo estimates of r* are '^^^-^2 ~ ^ ^"^^ 
r*^^24 = 3 in both underlying cases. That is, more severe association (larger e) suggests a larger choice of r in 
both cases. 

In Figure [Ml we present a Monte Carlo investigation of power based on asymptotic critical values against 
^V3/i2 ^tvz/2i ^ function of r for n ~ 10. In the AND-underlying case, the empirical significance level, 

q!„-io(?'), is about .05 for r = 2 and 3 which have the empirical power /?io(2) « .2 with maximum power at 
r = 2 for e = -\/3/12, and /3io(3) = 1 for e = 5-\/3/24. In the OR-underlying case, the empirical significance 
level, an=w{r), is closest to .05 for r = 1.5 which have the empirical power /3io(1.5) w .45 for e = %/3/12, 
and /3io(1.5) = 1 for e = 5^3/24. So. for small sample sizes, moderate values of r is more appropriate for 
normal approximation, as they yield the desired significance level, and the more severe the association, higher 
the power estimate. Furthermore, the OR-underlying version seems to perform better than the AND-underlying 
version for association alternatives. The empirical significance levels, and empirical power /3,f (r, e) values based 
on asymptotic critical values under for e = •\/3/12, 5 •\/3/24 are given in TableHl 



71 = 10 and N„ic = 1000 AND-underlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 




0.02 


O.I 


.2 


0.28 


0.35 


0.42 


0.73 


0.97 




0.023 


0.048 


0.035 


0.044 


0.040 


0.036 


0.031 


0.039 


/?L(V3/8) 


0.043 


0.109 


0.096 


0.153 


0.128 


0.119 


0.211 


0.287 


/?^c(V3/4) 


0.000 


0.98 


1 


1 


1 


1 


1 
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71 = 10 and N„ic = 1000 OR-undcrlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 




0.48 


0.48 


0.53 


0.62 


0.68 


0.73 


0.95 


1.00 




0.030 


0.045 


0.049 


0.043 


0.037 


0.043 


0.034 


0.000 


/?L(V3/8) 


0.028 


0.045 


0.059 


0.107 


0.113 


0.109 


0.151 


0.000 


/?^,(^/3/4) 


0.145 


0.681 


0.958 


0.998 


0.999 


0.999 


1.000 


0.000 



Table 1: The Monte Carlo critical values, C^^, empirical significance levels, a^^^{n), and empirical power 
estimates, P^^, based on the Monte Carlo critical values under H^^^ and H^^^, Nmc = 1000, and n ~ 10 at 
a = .05. 



77. = 10 and N,nc = 10000 AND-underlying case 


r 


1 


11/10 


6/5 


4/3 


a/2 


3/2 


2 


3 


as{n) 


0.2272 


0.2081 


0.1777 


0.1467 


0.1042 


0.1228 


0.0761 


0.0784 




0.3014 


0.4273 


0.4518 


0.4259 


0.3600 


0.4187 


0.3846 


0.5767 




0.6519 


0.9985 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


77 = 10 and Nmc = 10000 OR-imderlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 


as{n) 


0.2901 


0.1939 


0.2033 


0.1146 


0.0947 


0.0831 


0.0380 


0.0000 


M{r, V3/8) 


0.3182 


0.2621 


0.3135 


0.2601 


0.2466 


0.2554 


0.1594 


0.0000 


/3^(r, V3/4) 


0.7069 


0.9310 


0.9958 


1.0000 


1.0000 


0.9999 


1.0000 


0.0000 



Table 2: The empirical significance levels, as{n)^ and empirical power values, f3^{r,e), based on asymptotic 
critical values under iJf for e = a/3/8, a/3/4, Nmc = 10000, and tt, = 10 at a = .05. 

6 Multiple Triangle Case 

Suppose ym is a finite collection of tti > 3 points in M?. Consider the Delaunay triangulation (assumed to 
exist) of - Let Ti denote the 7*'' Delaunay triangle, Jm denote the number of triangles, and Ci/(3^m) denote 

the convex hull of ym- We wish to investigate Ho : Xi '~ l^iCniym)) against segregation and association 
alternatives using the relative edge densities of the associated underlying graphs. The underlying graphs are 
constructed using the PCD D, which is constructed using Np^{-) as described in Section [2T4l where for Xi e Tj, 
the three points in ym defining the Delaunay triangle Tj are used as 3^[j] . We consider various versions of the 
relative edge density as a test statistic in the multiple triangle case. 

6.1 First Version of Relative Edge Density in the Multiple Triangle Case 

For Jm > 1, as in SectionHH let pf"^(r) = 2 |fa„d| /{n {n-1)) and p°/(r) = 2 \£or\ /{n (ti- 1)). Let £f'"^ be the 
number of edges and p^'^'^{r) be the relative edge density for triangle i in the AND-underlying case, and £™ and 
p°.^(r) be similarly defined for OR-underlying case. Let ti, be the number of X points in Ti for 7 = 1,2,..., Jm- 
Letting Wi = A{Ti)/A{CH{ym)) with A{-) being the area functional, we obtain the following as a corollary to 
Theorem [3^ 

Corollary 6.1. The asymptotic null distribution for P/'jf (?') conditional on ym for r e (1, oo) is given by 

A/^(p?"fW-Ma„d(r)) ^ U {OA^andir)) , (14) 

where Jiand{r) = Hand{r) (Z]/=i^l) andvandir) = 



^d{r) (E.f=i wf) + {lland{r)f 



Figure 11: The empirical size (circles joined with solid lines) and power estimates (triangles with dotted lines) 
based on the asymptotic critical value against segregation alternatives in the AND-undcrlying case (top two) 
and the OR-underlying case (bottom two); in both cases, H^^^ (left) and H^^^ (right) as a function of r, for 

71 = 10 and N,nc = 10000. 



with ^landiT) and Vandij') being as in Equations and I110\) . respectively. The asymptotic null distribution of 
PTni'"') ''^'^^^ ^ [IjOo) is similar. 

The Proof is provided in Appendix 7. By an appropriate application of the Jensen's Inequality, we see 
that J2i=i W — (SfcTi . So the covariance above is zero iff Va^ndir) ~ and J2i=i — {j2i=i ' 
asymptotic normality may hold even though i'and('') = 0. The same holds for the OR-undcrlying case. 

Under the segregation (association) alternatives with S x 100% where S = 4e^/3 around the vertices of 
each triangle is forbidden (allowed), we obtain the above asymptotic distribution of P/"^(7') with /iand('') being 
replaced by /^and ('',£) and :^and('') by !^and('", e)- The OR-underlying case is similar. 



6.2 Other Versions of Relative Edge Density in the Multiple Triangle Case 

Let ET'ir) ~ P^;" (r) . Then E^^r) = pj^ir), since E^^r) = ^ "^^^ J 1) and(^) ^ 

n(n — 1) n (n — 1) 

Furthermore, let S,™*^ := '"'f Pn^'^C^) where Wi is as above. So S^"'^ a mixture of pf'."'^(?')'s. Then 
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Figure 12: Two Monte Carlo experiments against the association alternative H'^^^^. Depicted are kernel 

density estimates of for n = 10 (top left) and n = 100 (top right) and p°''(l.l) for n = 10 (bottom left) 

and n = 100 (bottom right) under the null (solid) and alternative (dashed). 



since p™'^(r)'s are asymptotically independent, ^J^"'^(r), pf"^(r) are asymptotically normal; i.e., for large n their 
distribution is approximately M [Jiandif) i '^i^andi'i') I n) . A similar result holds for the OR-underlying case. 

In Sect ion [01 the denominator of pf'i^(?') has n(n— 1)/2 as the maximum number of edges possible. However, 
by definition, given the n^'s we can at most have a graph with J,„ complete components, each with order rii for 
i 1, 2, . . . , J,„. Then the maximum number of edges possible is nj :— J2i=i {ni — l)/2 which suggests another 

\f il y^'^™ icandl Jul T] 

version of relative edge density: pT/Jr) := Then pff„(r) = ^'^^ ' ' — ^ ^ pt'^'^ir). 

nt ' nt ^ 2rit w 



Since ^^4^ ^ foi' each i, and ^ = 1' Pn%ir) is a mixture of p'^firYs. 

Theorem 6.2. The asymptotic null distribution for P^^^j^{t) conditional on for r G (l,oo) is given by 



' {Pniir) - fiandir)) ^ (0, 4 i^„„,(r)) , (15) 



where fLand{r) = fiand{r) and Vandir) 



Vand{r) (Ejt'l wf) I (E/ri ' 



with \Xand(f) and Vand{r) being as 



in Equations (0) and hlO^) . respectively. The asymptotic null distribution of pjj ^^{r) with r G [l,oo) is similar. 
Proof is provided in Appendix 8. Notice that the covariancc i'and(?') is zero iff fundir) ~ 0, Under the 



Figure 13: Empirical power estimates based on Monte Carlo critical values against the association alternatives 
with the AND-underlying case (top two) and OR- underlying case (bottom two), in both cases, H^^^^ (left) 

and (''^S^^) ^ ^ function of r, for n = 10 and N^c = 1000. 

segregation (association) alternatives, we obtain the above asymptotic distribution of pY}'\^{r) with /Xand('') 
being replaced by /iand('", e) and J^and(?') by fc'and(''j£)- The OR-underlying case is similar. 

Remark 6.3. Comparison of Versions of Relative Edge Density in the Multiple Triangle Case: Among 
the versions of the relative edge density we considered, E'^^^r) = /3f"^(r) for all n > 1, and 5^'^'^ and ^^"^(r) 
are asymptotically equivalent (i.e., they have the same asymptotic distribution in the limit). However, pf"^(r) 
and pf}'\^{r) do not have the same distribution for finite or infinite n. But we have pf"^(r) = „('^„"l\-| Pj"'^„(?") 

and ptand(?') < Mand(»') = Mand(?'), sincc "'j < 1- Furthermore, since = nir^-'i) — ' J2i=i 

we have lim„,_^oo Var[Vnpf"^(r)] = {j2i=i'^i^ li^n.^oo Var[^pf"jf (r)] Hence i>and(?') > J'and(r-). Therefore, 
we choose pf^^i^) further analysis in the multiple triangle case. Moreover, asymptotic normality might hold 
for pf;;f (r) even if i^,nd{r) = 0. □ 

6.3 Power Analysis for the Multiple Triangle Case 

Let S'^"''(r) pf^nir) and 5"(r) := /o"„(r). Thus in the case of J,n > 1 (i-C, m > 3), we have a (conditional) 

test of Ho : Xi *~ Z^(C//(3^„i)) which once again rejects against segregation for large values of S^^'^{r) and rejects 
against association for small values of S*^"*^ (r) . The same holds for 5*°'' (r) . 

Depicted in Figures fT5l and [TBI are the realizations of 100 and 1000 observations, respectively, independent 



n = 10 and N,nc = 1000 AND-underlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 


5 


10 




0.0 


0.0 


0.02 


0.06 


0.08 


O.I 


0.24 


0.46 


0.68 


0.82 


"mc(") 


0.000 


0.000 


0.005 


0.030 


0.027 


0.037 


0.038 


0.043 


0.048 


0.041 




0.000 


0.000 


0.003 


0.045 


0.057 


0.077 


0.154 


0.136 


0.077 


0.055 




0.000 


0.000 


0.009 


0.051 


0.060 


0.081 


0.492 


0.964 


0.941 


0.396 


n = 10 and N„ic = 1000 OR- underlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 


5 


10 




0.26 


0.26 


0.28 


0.31 


0.3 


0.35 


0.6 


0.84 


0.95 


1.00 


"™c(") 


0.000 


0.000 


0.040 


0.045 


0.049 


0.042 


0.049 


0.044 


0.022 


0.019 




0.000 


0.000 


0.169 


0.227 


0.331 


0.328 


0.396 


0.163 


0.069 


0.032 




0.000 


0.000 


0.000 


0.352 


0.352 


0.612 


0.988 


1.000 


0.935 


0.344 



Table 3: Monte Carlo critical values, C^^^, empirical significance levels, a-^^{n), and empirical power estimates, 
B:^^. based on Monte Carlo critical values under H^^.^ and i?"* /;t,„,, Nmc = 1000, and n = 10 at a = .05. 

'"-^ Vo/12 ovo/ZA 



n = 10 and N„ic = 1000 AND-underlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 


5 


10 


aAin) 


0.7707 


0.3343 


0.1872 


0.0859 


0.0774 


0.0671 


0.0551 


0.0593 


0.0771 


0.1182 


'(3^{r, V3/12) 


0.7406 


0.2829 


0.1869 


0.1156 


0.1323 


0.1506 


0.2053 


0.1599 


0.1336 


0.1618 


/?;^(r, 5^3/24) 


0.7415 


0.2923 


0.1833 


0.1220 


0.1491 


0.1891 


0.5605 


0.9664 


0.9510 


0.6241 


n = 10 and N„ic = 1000 OR- underlying case 


r 


1 


11/10 


6/5 


4/3 


V2 


3/2 


2 


3 


5 


10 


aA{n) 


0.5194 


0.3935 


0.2302 


0.0920 


0.0834 


0.0665 


0.0759 


0.0980 


0.0708 


0.0193 


/3;^(r,V3/12) 


0.6293 


0.6258 


0.5661 


0.4318 


0.4247 


0.4346 


0.4343 


0.2624 


0.1421 


0.0336 


/?;^(r, 5^3/24) 


0.6315 


0.6340 


0.6259 


0.6265 


0.6279 


0.7480 


0.9900 


1.0000 


0.9649 


0.3505 



Table 4: The empirical significance level and empirical power estimates based on asymptotic critical values 
under for e = %/3/12, 5 V3/24, Nmc = 10000, and n = 10 at a .05. 



identically distributed according to the segregation with 5 = 1/16, null, and association with 6 = 1/4 (from left 
to right) for |3^,„| = 10 and Jio = 13. 

With n = 100 ; for the null realization, the p-value is greater than 0.1 for all r except r = 1, 4/3, ^/2 for 
both alternatives in the AND-underlying case, and for all r values and both alternatives in the OR-underlying 
case. For the segregation realization with 6 = 1/16, we obtain p < 0.018 for all r values except r = 1 in the 
AND-underlying case and p < 0.02 for all r values in the OR-underlying case. For the association realization 
with S = 1/4. we obtain p < 0.043 for r = 2, 3 in the AND-underlying case and p < 0.05 for r = 4/3, \/2, 1.5, 2 
in the OR-underlying case. 

With n = 1000, in the AND-underlying case under the null distribution, p > .05 for all r values relative 
to segregation and association. Under segregation with 6 = 1/16, p < .01 for all r values considered. Under 
association with 6 = 1/4, p < .01 for r <E {4/3, v^, 1.5, 2, 3, 5} and p > .05 for the other r values considered. 
In the OR-underlying case under the null distribution, p > .05 for all r values relative to segregation and 
association. Under segregation with S = 1/16, p < .01 for r <E {1.1, 1.2, 4/3, v^, 1.5, 2, 3, 5} and p > .05 for 
the other r values considered. Under association with S = 1/4, p < .01 for r G {1.1, 1.2, 4/3, a/2, 1.5, 2, 3} and 
p > .05 for the other r values considered. 

We repeat the null realization 1000 times for n = 100 and find the estimated significance level above 0.05 for 
the AND-underlying case relative to both alternatives with smallest being 0.12 at r = 2 relative to segregation 
and 0.099 at r = 2 relative to association. The associated empirical size and power estimates are presented in 
Figures [T71 and [T51 These results indicate that n = 100 (i.e., the average number of points per triangle being 
about 8) is not enough for the normal approximation in the AND-underlying case. For the OR-underlying case 
the estimated significance level relative to segregation is closest to 0.05 is 0.03 at r = 5 and all much different 
at other r values. The estimated significance level relative to association are larger than 0.25 for all r values. 
Again the number of points per triangle is not large enough for normal approximation. With n = 500 (i.e.. 



Figure 14: The empirical size (circles joined with solid lines) and power estimates (triangles with dotted lines) 
based on the asymptotic critical value against association alternatives in the AND-undcrlying case (top two) 
and the OR-underlying case (bottom two), in both cases, H^^^^ (left) and H^^^^^ (right) as a function of r, 
for n = 10 and N,nc = 10000. 

the average number of points per triangle being about 40), the estimated significance levels get closer to 0.05, 
however they still are all above 0.05, hence for moderate sample sizes, the tests using the relative edge density of 
the underlying graphs are liberal in rejecting Hq. The empirical power analysis suggests the choice of r = 2 — a 
moderate r value — for both alternatives in both underlying cases. Note also that AND-underlying case seems 
to perform better for segregation. 

The PAE is given for J,„ = 1 in Section [42l For J,„ > 1, the analysis will depend both the number of 
triangles as well as the sizes of the triangles. So the optimal r values suggested for the Jm = 1 case does not 
necessarily hold for Jm > 1, so it needs to be updated, given the J^m points. The conditional test presented here 
is appropriate when the are fixed. An unconditional version requires the joint distribution of the number 
and size of De launay triangles wh en y,n is, for instance, a Poisson point pattern. Alas, this joint distribution is 
not available I Okabe et al. ( 2000l )). 



6.4 Extension to Higher Dimensions 

The extension to R'' for d > 2 is straightforward. Let yd+i = {yi,y2, ■ • ■ jYrf+i} be d + 1 non-coplanar points. 
Denote the simplex formed by these d + 1 points as 6(3^d+i). A simplex is the simplest polytope in having 
d + 1 vertices, c? (d + l)/2 edges and d + 1 faces of dimension (d — 1). For r S [1, oo], define the proportional- 
edge proximity map as follows. Given a point x in &{yd+i), let y := argminj,gj;^^j^ volume(Qj,(x)) where 
Qy{x) is the polytope with vertices being the d{d + l)/2 midpoints of the edges, the vertex y and x. That 
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Figure 15: Realization of segregation (left), Ho (middle), and association (right) for |3^m| — 10 and n = 100. 
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Figure 16: Realization of segregation (left). Ho (middle), and association (right) for |3^,„| = 10 and n = 1000. 

is, the vertex region for vertex v is the polytope with vertices given by v and the midpoints of the edges. 
Let v{x) be the vertex in whose region x falls. (If x falls on the boundary of two vertex regions or at the 
center of mass, we assign v{x) arbitrarily.) Let ^p{x) be the face opposite to vertex v{x), and •q{v{x),x) be 
the hyperplane parallel to ^{x) which contains x. Let d{v{x)Tri(v{x),x)) be the (perpendicular) Euclidean 
distance from v{x) to ri(v{x), x). For r € [l,oo), let 'r]r{v{x), x) be the hyperplane parallel to ip{x) such that 
d{v(x),rir(v{x),x)) ~ r d{v(x),ri(v{x),x)) and d{r]{v{x),x),'r]r{v{x),x)) < d{v{x),rir{v{x),x)). Let 6r{x) be the 
polytope similar to and with the same orientation as 6 having v(x) as a vertex and •qr{v{x), x) as the opposite 
face. Then the proportional-edge proximity region Np^{x) := &r{x) n S(3^d+i)- Furthermore, let C,i{x) be 
the hyperplane such that C,i{x) n S(3^d+i) 7^ and r d{yi, Ci{^)) — d{yi, ij{yi,x)) for i ~ 1, 2, . . . , d + 1. Then 
Tlix) n i?(y,) = {ze i?(yO : d{y,,r,{y„ z)) > d(y„0(a:)}, for i ^ 1,2,3. Hence r[{x) = W^lliTUx) n i?(y,)). 
Notice that r > 1 implies x G Np^{x) and x S T\{x). 

Theorem 1 generalizes, so that any simplex S in can be transformed into a regular polytope (with edges 
being equal in length and faces being equal in volume) preserving uniformity. Delaunay triangulation becomes 
Delaunay tessellation in W'-, provided no more than 4 points being cospherical (lying on the boundary of the 
same sphere). In particular, with d = 3, the general simplex is a tetrahedron (4 vertices, 4 triangular faces 
and 6 edges), which can be mapped into a regular tetrahedron (4 faces are equilateral triangles) with vertices 
(0, 0, 0) (1, 0, 0) (1/2, V3/2, 0), (1/2, ^3/4, ^3/2). 

Asymptotic normality of the [/-statistic and consistency of the tests hold for d > 2 in both underlying cases. 



Figure 17: The empirical size (circles joined with solid lines) and power estimates (triangles with dotted lines) 
based on the asymptotic critical value for the AND-underlying case (top) and the OR-undcrlying case (bottom) 
in the multiple triangle case, in both cases, H^r- (left) and HA- (right) as a function of r, for n = 100. 



7 Discussion and Conclusions 



In this article, we consider the asymptotic distribution of the relative edge density of the underlying graphs 
based on (parametrized) proportional-cdgc proximity catch digraphs (PCDs), for testing bivariate spatial point 
patterns of segregation and association. To our knowle dge the PCD-based metho d s are the only graph theo- 
retic methods for te sting spatial patterns in literature (jCevhan and Prie bi (l2005l ). ICevhan et all (|2006l) . and 
Cevhan et all (120071 )1. The proportional-cdgc PCDs lend themselves for such a purpose, because of the geom- 
etry invariance property for uniform data on Delaunay triangles. Let the two samples of sizes n and m be 
from classes X and 3^, respectively, with X points being used as the vertices of the PCDs and y points being 
used in the construction of Delaunay triangulation. For the relative density approach to be appropriate, n 
should be much larger compared to m. This implies that n tends to infinity while m is assumed to be fixed. 
That is, the difference in the relative abundance of the two classes should be large for this method. Such an 
imbalance usually confounds the results of other spatial interaction tests. Furthermore, we can perform Monte 
Carlo randomization to remove the conditioning on . 

Previouslv. lCevhan et al. ( 20061 ) employed the relative (arc) density of the proportional-edge PCDs for testing 
bivariate spatial patterns. In this work, we consider the AND- and OR-underlying graphs based on this PCD; in 
particular, we demonstrate that relative edge density of these underlying PCDs is a J7-statistic, and employing 
asymptotic normality of [/-statistics, we derive the asymptotic distribution of the relative edge density. We then 
use relative edge density as a test statistic for testing segregation and association. 



Figure 18: The empirical size (circles joined with solid lines) and power estimates (triangles with dotted lines) 
based on the asymptotic critical value for the AND-underlying case (top) and the OR-undcrlying case (bottom) 
in the multiple triangle case, in both cases, H^r- (left) and HA- (right) as a function of r, for n = 500. 



The null hypothesis is assumed to be CSR of X points, i.e., the uniformncss of X points in the convex hull of 
y points. Although wc have two classes here, the null pattern is not the CSR independence, since for finite to, 
we condition on to and the areas of the Delaunay triangles based on y points as long as they are not co-circular. 

There are many types of parametrizations for the alternatives. The particular parametrization of the alterna- 
tives in this article is chosen so that the distribution of the relative edge density under the alternatives would be 
geometry invariant (i.e., independent of the geometry of the support triangles). The more natural alternatives 
(i.e., the alternatives that arc more likely to be found in practice) can be similar to or might be approximated by 
our parametrization. Because in any segregation alternative, the X points will tend to be further away from y 
points and in any association alternative X points will tend to cluster around the y points. And such patterns 
can be detected by the test statistics based on the relative edge density, since under segregation (whether it is 
parametrized as in Section 13.21 or not) we expect them to be larger, and under association (regardless of the 
parametrization) they tend to be smaller. 

Our Monte Carlo simulation analysis and asymptotic efficiency analysis based on Pitman asymptotic effi- 
ciency reveals that AND-underlying graph has better power performance against segregation compared to the 
digraph and OR-underlying version. On the other hand, OR-underlying graph has better power performance 
against association compared to the digraph and AND-underlying version. When the number of X points per 
triangle is less than 30, we recommend the use Monte Carlo randomization, otherwise we recommend the use of 
normal approximation as n — cxd. Furthermore, when testing against segregation we recommend the parameter 
r « 2, while for testing against association we recommend the parameters r S (2, 3) as they exhibit the better 
performance in terms of size and power. 
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Figure 19: Var [hf^'^ir)] (left) and Var [h"^^^)] (right) ftmction of r for r € [1,5]. 

APPENDIX 

Appendix 1: The Variance of Relative Edge Density for the AND-Underlying 
Graph Version: 

The variance term is 



ine variance term is 

Var [hlfir)] = ipZ^{r)Iir € [1,4/3)) + ipY:^{r)Iir G [4/3, 3/2)) + ^-d(r)I(r e [3/2, 2)) + ^^(^^j^^ ^ [2,oo)) 

U and/ \ (5r'^-153r^+393r-'-423r^-54r^+360r-128)(447r''-261r^+54r^+5r'^-153r^+360r-128) 

Where ^ [r) 29167^(7+2)^(7+1? ' 

and/ \ _ (101 r°-801 r-^ + 1302 r-^-732 -536 r+672)(1518 r^-84 r^-104 r+101 r'^+672) 
'''1,2 1';— 46656 r^(r+2)^(r+l)^ ' 

and/ \ _ (r'^-13r'^+30 r'' + 148 r-'^ -448 r^+264 r''+288 -368 r+96) (22 r'^ + 124 -464 r^+r**- 13 r'^+264 r^+288 -368 r+96) 
¥'l,3 ~ 64r«(r+2y(r+l)^ ' 

^lf{r) = ('■^+'■^-3 r^-3r^+6^.-2K3 .=^+3 .^-6 r+2) gee FigurcIH Note that Var and (r = 1) = and hm,^oo Vara„d( 
(at rate 0{r^^)), and argsup^gj^ ^o) Var and (j') ~ 2.1126 with sup Varand('') 



iieie (/^i 1 (/ ; 2916r-i(r- 

and/ \ _ (101r°-801r-* + 1302r-^-732r^-536r+672)(1518r^-8-t < - , -r-Lu± < -ou-l -rui^; 
'''1,2 1';— 46656 r^(r+2)^(r+l)^ ' 

and/ \ _ (r'^-13r'^+30 r'' + 148 r-'^ -448 r^+264 r''+288 -368 r+96) (22 r'^ + 124 -464 r^+r**- 13 r'^+264 r^+288 -368 r+96) 
¥'l,3 ~ 64r«(r+2)^(r+l)^ ' 

Var and (r = 

id(0 = -25. 

Moreover, 

11 

^.and(r) Gov [htf{r),h^f{r)] = ^C'd(r)I(XO 

i=l 

where 
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Appendix 2: The Variance of Relative Edge Density for the OR-Underlying Graph 
Version: 

The variance term is 
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Figure 20: i'and(r-) = Gov [hf^"^ (r) , hf^"^ (r)] (left) and i/or(r) = Gov [hf2{r),h1l{r)] (right) as a function of r 
for r e [1, 5]. 
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67536 r^** - 511355 r" + 1161830 r^® - 634128 r^^ - 3001568 r^* + 9512164 r^^ - 11014136 r^^ + 2344968 r" + 7126240 r^°- 
13850504 + 14466592 r** - 3823216 - 4018976 + 5155776 - 4633984 + 1959808 - 244480 - 3584 r - 1024) 

^10 (0 - 960(2r2-l)(r + 2)3(r-l)(r + l)2ri° " " + + " '"^^^^ + + 

166787 - 717366 r^^ + 1196092 r^'^ + 89468 r" - 5130844 + 12748688 r^^ - 11274744 r" - 12243496 7-^" + 
33980568 - 14886656 r** - 19910592 + 20667776 r'^ - 1262208 - 5402752 + 2217088 - 235776 - 2560 r - 1024) 

2 180r-*-48r^- 648r'' + 396r5 + 214r*-190r3 + 39r2-4r + l 
iiW - 15 (2r-2-l)(r + l)Vio 

and Ji = [1,2/V3), I2 = [2/^3,6/5), X3 = [6/5,^/5-1), X4 = [Vs - 1, (6 + 2 ^/2)/7), X5 = [(6 + 
272)/?, 4/3), Xe = [4/3, (6 + VT5)/7), X7 = [(6 + yi5)/7, 3/2), = [3/2, (1 + V5)/2), Xg = [(1 + 
y5)/2,l + 1/^2), Xio = [1 + l/%/2,2), Xii = [2,00). See Figure EQl Note that CoVor(r = 1) = 1/3240 
and linir^oo Vo^(r) ~ (at rate 0(r~^)), and argsup^gjj^ Vo^(r) k. 1.765 with sup:^or('') ~ .0318. 

Appendix 3: Derivation of ;Uand(^) and z/and(^) under the Null Case 

In the standard equilateral triangle, let yi = (0,0), y2 = (1,0), ys = (l/2, V3/2), Mc be the center of mass. Mi 
be the midpoints of the edges for i = 1,2,3. Then Mc = (1/2, a/3/6). Mi = (3/4,^3/4), M2 = (1/4,^3/4), 
A/3 = (1/2,0). Let be a random sample of size n from Z^(T(3^3)). For xi = (m, u), i'r(a;i) = r w+r a/3w — a/Sx. 
Next, let iVi lr{xi) n 63 and iVa 4(a;i) H 62. 

Derivation of yUand(^) in Theorem 13.21 

First we find /^and('') for r e (l,cx)). Observe that, by symmetry, 

/^and(r) = P(X2 G iv^B(Xi) H r^Xi)) = 6 ^(^2 s ivj(Xi) H ruxi),xi e t,) 

where Tg is the triangle with vertices yi, M^-, and A/c*. Let ^^(r, x) be the line such that r (i(yi, £s(r, a;)) = 
'^(yijei)i so is{r,x) = V3{l/r — x). Then if .ti e T, is above ls{r,x) then iVp£;(.Ti) = ^(3^3), otherwise, 

N^PEi^i) c T(3;3). 

To compute /iand(f), we need to consider various cases for Np^{Xi) and T\{Xi) given Xi = (a;,y) G Tg. 
See Figures [2T] and [22l For any x = (u,?;) G X(3^), r5^(a;) is a convex or nonconvex polygon. Let ^i{r,x) be the 
line between x and the vertex y^ parallel to the edge such that rd{yi,£,i{r,x)) = d{yi,£r{x)) for i = 1,2,3. 
Then ri{x) n -R(yi) is bounded by ^i{r,x) and the median lines. For x = {u,v), ^i{r,x) = — VSa; + [v + 

V3u)/r, £,2{r,x) = (w + \/3r (x - 1) + ^3(1 - u))/r and 6(^,^) = (\/3(r - 1) + 2 v)/(2 r). For r e 6/5, V5-1), 

there are six cases regarding r^(a:) and one case for Np^[x). See Figure [22] for the prototypes of these six 
cases of Fi {x,Ny). For the AND-underlying version, we determine the possible types of Np^{xi) fl r^(xi) 
for xi G Tg. Depending on the location of xi and the value of the parameter r, Np^{xi) n TKxi) regions are 
polygons with various vertices. See Figure [24] for the illustration of these vertices and below for their explicit 
forms. 

^ _ f V3y+3x n\ r< — f V3y-3r+3-3 x n\ n — ( ^/3y-(, r+3-3x \/3(-%/3y-3+3x) \ _ ( [V3r+V3-2 
- \ — 37 — '"J' '-'S - 1^ 3;^ ,0J, Cr3 - I ^ , g7 I, <-^4 - q;. 

y3(3r-^3+2\/3j^\ f~i _ ( {V3r~^+2y)V3 ^/%3r-3+2^^3y)\ ^ _ / ^3^+3^ \^{V3y+3x) \ ^ 
6r / V ' ^r ) ' — I 6r ' 6r j ' 

Pi = (l/2,V3/6 {2V3ry + 6rx~3)), and P2 = (-1/2 + (VSr y + 3 r a;)/2, - V3/6 (-3 + VSr y + 3rx)); 




Figure 21: The cases for relative position of (sif, x) with various r values. These are the prototypes for various 
types of Np^{xi). 



T fl/r, V3{2V3y+6x-3r)\ ( (-2 \/3y-6+6 a:+3 r)\/3\ / v^j^.g ^+3.3 ^{ir—Jiy-Z+Zx) 
^1 - I 1/^, yl-l- 11/^, 67 I , ^3 - I Yt ' 

^ _ ( 3r-3+2 ^f?,y \/3(3 r-3+2 V3iy) \ ^ ( r-3+2 \/3y \/3(3 r-3+2 \/3i;) \ ^ _ ( -r-+\/3j/+3a \/3(V3iy+3 a:-3 r) ^ ^ 

A^i = (V3ry/3 + ra;,0), N2 ^ (V3ry/6 + rx/2, V3 (V3y/6 + 3x) r), and 

= (V3r y/4 + 3 r x/4, ^/3 (^/3y/12 + 3 x) r) ; and = (^ x/3'--y+3^-x-V3,+3.-3+3. ^ (v/3.^,+3.^.+^^^,-3.+3-3.)^/3 



and Q2 



2 ^/3r^y+6r^x-3r+3-2 ^/3y \/3(3r-3+2 \/3y) 
6r ' 6r 



Let ,f3^(ai, a2, . . . , a„) denote the polygon with vertices ai, a2, . . . , a„. For r e [l, 4/3), there are 14 cases to 
consider for calculation of /iand(?') in the AND-underlying version. Each of these cases correspond to the regions 
in Figure [26l where Case 1 corresponds to Ri for i = 1,2,3,4, and Case j for j > 1 corresponds to Rj+3 for 
j = 1, 2, . . . , 14. These regions are bounded by various combinations of the lines defined below. 

Let £am{x) be the line joining yi to Mc, then £a7n{'x) = V3x/3. Let also ri{x) ~ \/3 (2 r + 3 a; — 3) /3, r2{x) — 

^/3/2 - V3r/3, r^ix) ^{2x-2 + r) ^/3/2, nix) = V3/2 - x/3r/4, r,{x) = r^ix) = - ^^-^3+^"^ , 

, , (l+r'^x-r-x)V3 , , (r'^x-l+x)V3 , . ^r'^x-\)^^i , , (-2 r+2+r^x) \/3 , , 

riKx) = ^2+1 — , ''sW = ;23T^ — , rg[x) = -^^3+^ — , rio[x) = ,4+^2 — , ''iila;) = 

/'_2 r+2 — 2 :r+r^a:)\/3 / 

— -^^ :p2^ — , ri2(a::) = — {2x — r) v3/2, and ri^{x) = — (—1 + x) v3/3. Furthermore, to determine the 

integration limits, we specify the x-coordinate of the boundaries of these regions using Sk for k = 0, 1, . . . , 14. 
See also Figure [26] for an illustration of these points whose explicit forms are provided below. 

so = 1 - 2r/3, si = 3/2 - r, = 3/(8 r), S3 = ''^'/'^' ^ S4 = 1 - r/2, = ^^i^, - 1/(2 r), 57 - 
2T2F^' ^8 = ^-^^t^+Vr^^ -^9 = 1/ (r + 1), sio = ^^Hl^, -sii = 3r/8, s,^ - ^13 = 3/2-5r/6, 

and si4 = r - 1/2 - r^/8. 

Below, we compute P{X2 £ Np^{Xi) n ri(Xi),Xi e T^) for each of the 14 cases: Case 1: 



Jo J B2 Jo 



A{^{ Gi,Ni,N2 ,G6)) 

(r - 1) (r + 1) (r-^ + l) 
64 r6 



where A(^(Gi, M'l, A^2, Ge)) = ^/3/36 (\/3j/ + 3 2:) 



^(^i;+33:)^ 



36 r2 



y, = (l/2,\/3/2) 



V, = (1/2, \/3/2) 




M3 



y?^, 0) yi = (0,0) 



Figure 22: The prototypes of the six cases of (x) for a; G for r e [1, 4/3). 



Case 2: 



P{X2 G iv;B(Xi) nn(Xi),Xi G T,) = 



yjss Jrr,{x) Jsfi Jo J 



A{g^{Gi,Ni,P2,M3,Ge)) 



dydx - 



(9 + 23 r* + 24 + 24 + 13 r + 3) (r - 1)* 



96 r6 (r + 1)^ 

where yl(^(Gi, A^^i, P2, M3, Gej) = ^ITi ■ 

Case 3: 



P{X2 G a^;b(Xi) n r5:(Xi),Xi g t,) 



/ / +/ / +/ / 

\Jss J r-r{x) Jag JO J ai2 JO J 



A{^{Gi,G2,Qi,P2,M3,Ge)) ^, , 



324 r" - 1620 - 618 + 4626 r*^ + 990 r'^ - 2454 + 2703 - 5571 r"* - 3827 + 1455 + 3072 r + 1024 

7776 (r + 1)^ r6 

where yl(^(Gi, G2, Qi, P2, Afs, Ge)) = - [V3(-4 VSr y - 12 a; + 4 + 4 - 12r + 9r^ + 12ra: + 4r^y^ - 12xV^ ■ 
24 r''x + 12 r^x^ + 8 r"*^?/ a; + 12 a;^ + 12 r^x + 6 - 8 r^^i/ + 4 + 4 ^r^y) j / |^24 



Case 4: 



I rs(x) J 31^ 1-3(2:) 



[512 + 138240 + 3654 r^^ - 255 + 43008 - 12369 - 86387 r'' - 193581 + 148224 - 100608 + 94802 r"^ - 

35328 r^^]/ [7776 (r" + 1)%^] 

where ^(^(Gi, Mi, La, Qi, P2, M3, Ge)) = -[^(6x + 3r2 - 2 VSy + 2V3r^y + 2r'^ y'^ - 4r-^V3y + 4V3y x + 2r^y'^ + 
4 r*^/3y a; - 6 - 12 r^i + 6 r*x'^ + 6 r^x - 3)1 / [l2 1 



Case 5: 

Pix,eNMX^)nr,ix,),x,en)-^J^^ + J^^ y^^^^J ^^^^^^^^^^ dydx- 

(177 - 648r^ + 570r^ - 360 r''^ + 28 r'' - 24r^ + 174r2 + 72r + 27) (-12r + 77-^ + 3)^ 

7776 (r2 + ifr^ 

where yi(^(Gi, A/i, L2, Qi, P2, M3, Gej) = -■ 



Case 6 



r-^ /"^-("'^ ^«o™(-)^ r=B rrsWX A(^(Gi, Ml, Pi, P2, Ma, Ge)) ^^.^^ _ 

; A{T{y3)y 

[137472 r^*^ - 952704 r'^'^ + 2792712 - 5116608 r^^ + 7057828 - 7725792 r^^ + 7022682 r^^ - 5484816 r" 



+ 



3631995 - 2213712 + 1213271 - 578976 + 292518 r'' - 101952 + 36612 r"* - 11664 + 3051 - 1296 r + 243 

[(l5552rVl)'(2r2 + l)'/] 

where yl(3^(Gi, Mi , Pi, P2, M3, Ge)) = '-. 

Case 7: 

P^Y P /V Inr-rY 1 Y f=Tl- T'^^V T" /""'^'^ A(^(Gi, Mi, L2, Qi, P2, Kh, Ge)) _ 

P(X2GiVp,(Xi)nri(Xi),XiGr.)-^^ |^^^J ^^.^^^^^^ dydx- 

4 (100 r" - 408 r^" +454 - 564 + 283 - 108 r'* - 34r-^ + 204 r--* -r'' + 132r^ + 26r + 24) (2r - 1)^ (r - 1)^ 

243 (r2 + l)=*r3 (2r2 + 1)^ 

where ^(^(Gi, Mi, L2, Qi, P2, M3, Ge)) = x + ?, r'^ - 2 VSy + 2 VSr^y + 2r'^ y'^ - 4r^V3y + 4V3y x + 2r'^y'^ + 

4 r*^/3y x-6 - 12 r^s + 6 r^s^ + 6 r^x - 3) j / [l2 r^j . 
Case 8: 

[(-2 + r-) (2369 - 11342 + 29934 - 50340 r** + 54056 - 51824 r'' + 48320 r'' - 20864/ - 640 

- 1280 + 512 r + 1024) j / [l5552 r**] 

where A(^(Gi, G2, Qi, iVs, Mc, M3, Ge)) = -^V3{4V3r''y-12x-12r + 5r^ + 12rx + 4y^-12x^r^ + 4r''y'' + r^ + 
2r''x/3yx-4r='V3y + 6- 12 r^x + 3 r*a;^ + 12 + 12r^x - 4 ^/3r y + 4 \/3y)] / [24r2] . 



Case 9: 



S12 rr2(x) rsi3 rr2(x) 



'.10 ireW -^ = 12 Jr3{^) / ^(^(3^3))^ 



A(.9>{Gi,M^,L2,Qi,N3, Mc, M3,G^)) 

A{T{y,)y 

(49 - 168 + 354 - 528 r"^ + 236 r'' - 96r^ - 224 + 384 r + 64) (-12 r + 7 + 4)'' 

15552 r6 

where A(^(Gi, Mi, La, Qi, A'^a, Mc, M3, Ge)) = - [^/3(8^/3ya; + 4^/3r2y + 12 a: + 2r2 - 122:^^ - 4r3^/3^/ - 12r^x + 
3 r^x^ + r'^y^ + 2 r'^V^y a; + 12 r^a; - 6 - 4 VSy + 4 y^) j / ^24 r^j . 
Case 10: 

P{X2 e iV^.B(Xi)nn(Xi),Xi e T,) = 

/"-loC.)^ rsrs rri2(^) ^ ri/2 |■r^2(.)\ Mi, La, Qi, iVa, L4, L5, M3, Gs)) ^^^^ _ 

6144 + 195456 r-*^ + 324 r" - 76720 - 801792 + 217856 r + 946432 - 239904 - 275328 + 39408 r** - 11849 

31104 r3 

where A(£P(Gi,Afi,La,Qi,Af3,L4,L5,M3,G6)) = - [^/3(4^/3r2y + 8 v^ya; + 4r2 y^^mV3ry-4:r''VSy-2Ay^ + 12x + 
24r- er^ - 12 xV^ - 12 r^x + 3 r'^x^ + 12 r^x + 20 + 2 r^'^i/ x + r"* - 24)j / [24r-^j . 
Case 11: 

P{X2 G iV;.B(Xi)nn(Xi),Xi G Ts) = 

) /-.i* |^r^2(x)\ A{^^{Gi,Mi,L2,Qi,Q2,L5,M3,Ge)) , , 
I 'II A(T(V ))2 ^ " 

(r - 1) (1080 r^'^ + 1080 r^^ - 17820 - 540 r^^ + 65394 r^^ - 46926 r" + 105435 - 261765 r'' + 229286 - 180586 r^+ 



101638 r-*^ + 40774 r'^ - 46112/ + 24448 - 20224 + 10496 r - 6144)] / |^10368r'^ (2r^ + l)''* 

where^(,^(Gi,Mi,L2,Qi,Qa,L5,M3,G6)) = - [\/3(6x + 3r2-4r-2x^?/-4y2-6xV^ + 2r''^yx + 4^?/x-2r2 
4 r'^VSy + r^y'^ - 12 r^x + 3 r'-'x^ + 12 r^x - 6 + 4 V^r'^y + 2 ^/3y)] / \l2 r^] . 
Case 12: 

P(y aN^ (Y\nr^(X\ Y aT\ f'^' T'^^"' yl(^(Gi, Ga, Qi, jV3, L4, L5, M3, G5)) 

P(Aa G Npe{Xi) n Ti Xi),Xi G T^O = / / s^a = 



(49 r** - 204 + 476 r--* - 768 - 8 + 768 r - 288) (-6 + 5rf 

777672 

where A(^(Gi,Ga,Qi,Ar3,L4,L5,M3,G6)) = - [\/3(-12x + 12r■-3r■2 + 12rx-20^/3ry-12xV^+4^/3r2y-12r3x + 
3r''x^ + 28 V3y + 12x2 _^ Ur'^x - 12 - 20 + 4 - 4 r^\/3y + r* + 2 r*\/3?/x)j / [247-^ 
Case 13: 

ni/2 rrw(x) A{^{Li,L2,Qi,N3,Li,Ls,L6)) 



rL/2 r 

P{X2 G n^pe{Xi) n r^(Xi),Xi g t,) = / / 

J s^A J r 



^ri,{.) A{T{y,)Y '^^'^^ 

(4r-^ + 8r-^ - 37 r'^ - 58r* - 84r^ + 168r2 +336r - 352) (-2 + r) (r^ + 2r - 4)^ 

384 (r + 2)^2 

where A(^(Li, La, Qi,iV3,L4,L5,L6)) = - [\/3(-4 r^x/3y-8 \/3r y + 12 x + 24 r-8 \/3y x- 12 + 24 r x- 24- 12 xV^ + 
4 VSr^y - 32y2 - 12r^x + 3r'*x2 + 20 ^y- 24x2 _^ 12 r^x + 2 r^x/Sy x + r'' y^ + 4r2 y^) ^|]24r2j. 
Case 14: 

(135 r" +675r^° - 1350 - 9450 + 702 39150 + 24272 - 47432 r* - 135040 + 57088 + 204800 r- 

134144) (r - 1)] / 1^10368 (r + 2f r^] 



where L2, Qi, Q2, 1-s, ^e)) = - ^V3{-4:r^^y + 4:^/3r y + r* +6 x - A^Sy x + 2r'^^/3y x + 12r x - Ar^x - 

6 + 4 ^/3r^^/ - 12 r^x + 3 r^s^ + 2 v^Sj/ - 12 a;^ + 12 r^s - 6 - 8 - 2 2/^)] / ["^^ 



Adding up the P{X2 G Np^{Xi) n r']^(Xi), Xi e T^) values in the 14 possible cases above, and multiplying 
by 6 we get for r G [1, 4/3), 

(r - 1) (5 - 148 + 245 - 178 - 232 r + 128) 
''""^^'^^ " 54 r2 (r + 2) (r+1) " 

The /Xand('") values for the other intervals can be calculated similarly. For r = 00, /iand('') = 1 follows trivially. 
Derivation of Uand{r) in Theorem 13.21 

By symmetry, P{{X2,X3} C iV^B(Xi) n r^Xi)) = 6P({X2,X3} C A^^i;(Xi) n ^(^1), Xi G T,). 

For r G [6/5, V5— l), there are 14 cases to consider for calculation of ;^and(?') in the AND-underlying version: 
Case 1: 

P({X2,X3} c iVp,(X0 n r,(x,), e t.) = ^ | j ^^^^^^^^^^ dydx = 

(r^ + l)'(r + l)'(r-l)' 
384 rio 

(v^H + 3a-)^y3 



where ^(^(Gi, iVi, iVa, Ge)) = ^/3 {y/Zy + 3 x) rV36 3^;^- 

Case 2: 



P{{X2,Xs}cN^PEiX^)n^l{X^),X^eTs)=i / +/ / 

\J J r^{x) J ss Jo 



/ + 



(5 + 38 r + 137 + 320 + 552 r"* + 736 + 792 + 640 + 407 + 178 + 35 r^") (-1 + r)^ 

960 (r + 1)^ 

where yl(3^(Gi, A^i, P2, M3, Gs)) = 247^ -. 

Case 3: 

P({X2,Xs} c iV^E(Xi)nr5:(Xi), Xi e r,) = 

- /^^^ r\ r r'] ^(^(^-■^//^^-^^-^-^^.Ga))- 

- [l7496 - 122472 r + 139968 r" + 524880 r ^® - 553095 r^^ - 595971 r " + 368826 r^^ - 724758 r^^ - 543876 r " - 
1416996 r"' + 1646470 + 92870 r** + 523048 r'^ - 768368 - 1729902 - 1434990 r"* + 122185 + 941941 + 

573440 r + 114688] / [2099520 (r + 1)^ r"^] 

where A(^(Gi, G2, Qi, P2, M3, Ge)) = -[V3(4^/3r2y - Sr^VSy + 4r'y'^ +4r*y^ + 4y^ + Sr^V3yx + 6 - 12 s:^^ - 
12 a;- 12r- 24r^a; + 12r'*a:^ + 9r^ + 12 r x - 4 V3r y + 12 + 4 + 12 r^a;) j / ^24 r^] . 
Case 4: 

P({X2,X3} c Ni,E{Xi)nrl{Xi), Xi e r,) = 

=5 |-.io rsi2 rre(^)\ A{3^{Gi, Mi, L2, Qi, P2, M3, Gs))' ^^.^^ _ 

[32768 - 409264128 + 1455989508 r^^ + 680709729 r** - 4423680 + 155509 + 22889801 + 202936917 
6011901 + 1060982949 r^'^ - 614739456 + 240330993 r^^ - 56097792 - 77783040 - 999857664 + 
1299257316 - 1461851136 - 1407624192 r" + 1414729905 r" - 1352392704 r^^j / [2099520 (r^ + l) ^ r"^] 



where A{^{Gi, Mi, L2,Qi, P2, M3, Gg)) = ~\V3{-G x^r'^ - 3 + 6 x - 12 r^x + 6 r*x'^ -4r^ VSy + 4 V^y x + 4 r-^V^y x + 
2 r'^y^ + 3 + 2 %/3r^y - 2 VSy + 2 r^y^ + 6 r^i) j / |^12 r^j . 
Case 5: 

[(35361 r'^^- 229392 r^^ +602820 r" - 858384 r^^ + 778848 r^^ - 460368 r" + 277740 r " - 258768 + 160594 - 62256 r^- 
5892 - 17712 + 19224 r"* + 11664 + 5076 + 1296 r + 405) (-12 r + 7r-^ + 3)^j / ^699840 (r^ + l)^] 



where yl(^(Gi, A/i , Pi, P2, M3, Gsj) = jaT^ - 

Case 6: 

p({X2,X3} c iv;B(Xi)nr5:(Xi), g t,) = 

ea„^(x) rsr rta,„(x} 



'1), G T,) = 

.3 l-ia„Ax) ^ rsr rl^,^(x) ^ rss rrs(x)\ Aj^jGl, Mi, Pr, P2, Ms, Ge))\ ^^^^ _ 
^ . „, , J S3 Jv2{x) Js-r Jr2(x) ) ^(^(J^s))^ 

- [3645 - 17496 r + 5003898912 r^® + 31646646384 r^*^ + 110098944 r^" - 1090803456 r^^ - 14630751360 r^^ + 66339 r^- 

99072645696 r^^ + 79269457632 r^'' + 66073158 - 4870743552 - 168073488 + 535086 - 262440 - 1737936 r^- 

18592416 - 107383563504 r^^ - 41219053272 r^^ + 58981892347 r^** - 78265758888 r^^ + 95887286866 r^°+ 
109053166552 r^^ + 5500548 + 466565130 - 1070573040 r" + 2380992104 r^^ + 9191633420 r^"* - 16312513248 r-^^+ 

26801184917 - 54759787776 r^^j / [l399680 (r^ + 1)'^ (2r^ + l)%"'j 

where yl(^(Gi, Mi , Pi, P2, M3, Gej) = ^ 1275 -■ 

Case 7: 

^({^2,^3} c N'pE{Xi)r\Tl{Xi), Xi G r,) = 

/ rs^ nMx) ^ ..10 /rgWX A{3^{Gl,Mi,L2,Ql,P2,M3,G&)f ^ 

\J s-j J rs{x) Jss Jr2(x) ; A{T{y3)r 
[4 (162576 r^^ - 1083456 r^^ + 3368016 - 6969888 + 11578088 r^** - 15664080 r^^ + 18796852 r^*^ - 19984824 r^'^+ 
19534445 r-"-18170472 r^^ + 15507752 r^^-13150464 r"+9987958 r^°-7448736 r^+5016464 r**-2991768 r^+1857485 r®- 
749160 r-^ + 481804 r'* - 96720 + 76160 - 4032 r + 4320) (2r - 1)^ [r - 1)^]/ [32805 (r^ + l)%® (2r^ + l)^] 

where A(^(Gi, Mi, L2, Qi, P2, M3, Ge)) = -^V3{-Gx''r'^ - 3 + 6x - Ur'^x + 6r'^x'' - Ar^VSy + 4V3y x + Ar^VSy x + 
2 r'^y^ + 3 + 2 VSr^y - 2 V3y + 2 r^y^ + 6 r^x) j / [l2 r^j . 



Case 8 

P({X2,X3} c iV^B(Xi)nr^(Xi), Xi e r,) = 

^ A(J^'(Gi,G2,Qi,A^3,Mc,M3,G6))' ^ ^ 

A(T(3^3))=' '^^"'^ " 



S13 fraix) rl/2 rr2(x)'^ 

/ +/ / 

S12 Jreix) J si3 J re(x) 



- [-458752 + 811008 + 329205504 - 582626304 r^^ - 489563136 - 65536 - 168708096 r'' - 57883680 r"+ 
18009258 r-^** - 3623400 r^^ + 352563 + 41502720 r'^ + 659111904 - 761846400 r" + 725173376 r^^ + 409477188 r"- 

254829600 r^^ + 135968852 r^'^] / [8398080 r^°] 

where ^(^(Gi, G2, Qi, A^s, Mc, M3, Ge)) = -[^(-12x^2 - 12s:- 12 r- 12 r^a; + 3 r^'x^ + 4 \/3r22/ + 5 + 12 r x + 
12 + 2 r'^VSy x + 4 r'^y'^ - 4 r^VSy + 6 + 4 + r'^y^ + 4 VSy + 12 r^x - 4 VSr y) j / [24 r^] . 



Case 9: 

P{{X2,X3} c iv;E(Xi)nr^(Xi), Xi e r,) = 

170226 r" - 392112 r^'^ + 680784 r^^ - 1040256 r" + 1385628 r^" - 1337760 + 816224 - 253824 
469088 - 1029888 + 820992 - 488448 + 190976 + 49152 r + 8192) (- 12 r + 7 + 4) ^ j / ^8398080 r^°] 

where A(^i2i(G'i, Mi, La, Qi, iVs, Mc, Mg, Ge)) = - [\/3(-12 xV^ _ 6 + 12a; - 12 r^'a; + 3 r*x'^ + 2r^ +2r*V3yx + r^y'^ + 
SVSyx + Ar^y^ - AVSy + ^VSr'^y + 12 r^a; - 4 VSy) j / ^24 r^j . 
Case 10: 

p{{X2,Xs} c iv;B(Xi)nr5:(Xi), Xj e r,) = 

/ ..14 .rio(.) ^ ..13 .ri,(.) ^ .1/2 rrM-)\ {G I , Ah , L2 , Q 1 , Ns , L4 , L5 , M; , G6))\ y_^^ ^ 

\J sig Jr2{x) J Bi4 J r2{x) J S13 J r^ix) ; A(r(^3))« 

1^4423680 - 4627454976 + 511684992 r" + 2163142656 r'^ - 660127744 - 31555584 r + 3534520320 + 7647989760 
7785504 r^^ - 1313880 + 19683 r^* - 7240624128 - 1511047552 + 1204122240 - 796453824 - 282583320 r 

107804736 - 30362052 r^^j / [l6796160 r®] 

where A{.0»{Gi, Ah, L2, Qi, N3, L4, L5, Ms, Ge)) = - [\/3(-16 \/3r ?/ + 20 \/3y - 24^^ - 12 3:^^ + 12a: + 24r - 12r='a; + 
3 r^x'^ -Gr'^ -24 + 4 VSr^y + 8 VSy x-4 r^y/iy + 4 r^y^ + r^y^ + 2 r■'*^/3y a; + 12 r^a;) j |^24 r^] . 
Case 11 



P({X2,X3} C AfpE(Xi)nn(Xi), Xi G Ts) = 

/ rsil rla^x) ^ rs,„ rr^2(x) ^ rs^, .ri2(r.)\ yl( ^ , Ml , La , Ql, Qa, L5, M3, Gfi))' ^^^^ 

- [(r - 1) (-1474560 + 8847360 r + 111456 r^^ + 111456 r^'' - 27738112 + 23311152 r^^ - 167184 r^* - 808889416 r*- 
2228253688 r^'V 366739256 r^-207619072 ^+98557952 r^+397199360 r''^+802401664 r^-34733448 r^^ -624736557 r"+ 

400615470 r^** - 134938386 r^^ + 39014136 - 18026064 r^^ - 640058432 + 407655352 - 1227078728 r" + 
1996721576 r^^ + 2033409092 - 1681870468 r^'^ + 1064030499 r^'' - 2842128 r'^^)]/ [l866240 (2 + l) ^ r'^j 

where^(^(Gi,Mi,La,Qi,Qa,L5,M3,G6)) = - [^(4 ^r^y + 4 a; — 2 r^y^ — 4r^ VSy - 4 - 4 VSr'^y x-6 x'^r^ + 
6 s - 12 r^a; + 3 r*a;^ + 3 + 2 r'^VSy x + r'^y^ + 2 VSy + 12 r'^x [l2 r^] . 

Case 12: 

P({Xa, X3} c N^MX.) n r^(Xi), x. e r.) ^ £^ AonG^^G^^ ^ 

[(2322432 - 7554816 r + 9510912 + 1046068 r** - 558720 + 2444224/ - 5799360 - 2134656 - 1608672 r'^+ 

2169696 + 216300 r^" - 55440 + 7095 r^^) (-6 + 5 r)^] / [4199040 r"] 

where A(,^(Gi, Ga, Qi, iV3, L4, L5, M3, Ge)) = - [\/3(-12a:V^ - 12a; + 12 r - 12 r^a; + 3 r'^a;^ - 3?-^ + 12 r a- + 28 + 
12 a;^ - 20 + 12 r^a; + r*y^ + 4 r^y^ - 4 V3y + 2 r''^/3y a; + 4 VSr^y - 20 \/3r y - 12)] / [24 r^j . 
Case 13: 

P({Y y\rN^ (y^nV^(Y\ y aT\ f'^' A(^(Li,La,Qi,iV3,L4,L5,L6))% 
P(|Aa,A3| C iVpB(Ai) nri(Ai), Ai e Ts) = / / Afrrr^j \\3 dyda; = 

- [(9 r^" + 36 - 132 r^^ - 576 + 164 + 2512 + 4976 r** - 1536 - 13888 r'' - 17536 r'^ - 3072 + 79360 + 

9216 - 120832 r + 61440) (-2 + r) (r^ + 2 r - 4)^] / [768O (r + 2)'S-''j 



5r^ + 



J^and(r) 

2832544 



where La, Qi,A''3, 1.4,1/5,^6)) = - \y^{r^y'^ -9>V^r y-^y/^y x + 4, r^y^ ~ A r^y/?,y -2,2 y'^ + 2 r^^/^y x~l2 

12x + 24r -12r^x + 3 r^x^ - 12 + 4 ^r^y + 24 r a: - 24 a;^ - 24 + 20 ypiy + 12 r^x) j / ^24 r^j . 
Case 14: 

P({X2,X3} c iV^B(Xi)nr5:(Xi), e r,) = 

[(r - 1) (3483 r + 24381 r^^ - 34830 r^^ - 529416 r - 265680 r^^ + 4274208 r + 4999320 r - 15227352 r " - 
25751336 r^" + 19466488 + 62834064 r" + 17452256 r'' - 53339200 r'' - 117114624 - 51206656 + 270430208 r' 

58073088 - 296222720 r + 122159104) j j [l866240 (r + 2)^ r"*! 

where 4(^(Li, L2, Qi, Q2, is, Le)) = - [^(-4 a; - 2 r^y^ + 4 V3r y - 4 r^^y - 8y^ -4^r^yx - 6a;V^ +6a; - 

12 r^x + 3 r--*x^ + 4 ^r^y + 12 r- X- - 12 x^ + 2 r^'^y x + r^'y^ + 2 + 12 r^x - 6)]/[l2r^]. 

Adding up the P{\Xi^X-^ C N^p^(X\) n r5^(Xi), X\ e Ts) values in the 14 possible cases above, and 
multiplying by 6 we get for r € [6/5, ~ l): 

= - [219936 r-3041936 r2-30889822 r^+18084672 ri3+27137438 r^+2364868 r^+2305864 r3-4168820 r^- 
'^+486 r2i-118850 r"-45155 ri^-269 ri^+3402 r^o+lllOlieO r^+24604048 r^°-43009544 r"+8770788 r^^ 
13736295 ri^ + 2751855 r^^ + 443518 ri^ + 49152 j 116640 (r + 2)^ (2 + l) (r + 1)^ . 

The fand('') values for the other intervals can be calculated similarly. 

Appendix 4: Derivation of /ior(^) and i^or(^) under the Null Case 
Derivation of fj.or{r) in Theorem 13.21 

First we find fioiif) for r £ l,oo). Observe that, by symmetry, 

/.orW = P{X2 e n^pe{Xi) u riiXi)) = 6P{X2 e iVJ(Xi) u T[{Xi),x, e t,). 

For r € [1,4/3), there are 17 cases to consider for calculation of lyoiir) in the OR-underlying version. Each 
Case j correspond to Ri for i = l,2,...,17in Figure Case 1: 

P{x, e iVp,(xo u r,(xo,x. e r.) = ^ +J^^ j ^^.^^3^^^^^ dydx = 

^r'-4r/9 + l/3 

where A{&>{A, Mi, Mc , Ivh)) = ^/3/12. 
Case 2: 

P(X2 e iV^.B(Xi)un(Xi),Xi G r«) = 

/ .SI /TiC.)^ .S3 rr^{..) ^ rs, rr,(.) ^ fs, rr,[.)\ A{3^ {A, Ml , L2 , Pj , Mq , M3)) ^ ^ 
VJsq •'0 J si J -/ S3 -/ 

(r - 1) (1817 r'' - 7807 + 14157 - 14067 r"^ + 7893 - 2475 + 405 r - 27) 

86476 

where A{3^{A, Mi, L2, L3, Mc , Ms)) = — ^ '-. 



Case 3: 



[Js, Jo A. Jo J A{T{y,)Y 

(iSr'' -4r-^ +4r- - 1 - 2r^) (r- 1)" 



where G2, G3, M2, Mc, M3)) = ^ '-. 

Case 4: 

P(X ^N^ (X\,^r^(X \ X ^T\-( r f'-^''\ r r° A ^ ( ^ Ml , La , L3 , L4 , L5 , M3 ) ) 

P(x.eiVp,(X0ur,(X0,XieT.)-^^ y^^^^^ + j^^ j^^^, dyd.- 

(9 - 72 r + 192 - 192 r» + 76 r") (4 r - 3 + \/3)^ (4 r - 3 - 

10368 r6 

where ^(^(^, Mi, L2, L3, ^4, is, M3)) = — ^ j^T^ ^. 

Case 5: 



/ rr,{x) rs, rr,(x)\ A{^?» {A, G2 , G3 , M2 , Mg , P2, N2)) , , 

[Js, +4 io J AiTiy,)r ''^^"^ 



^1 + 2r + 6r-^ - er'^ + 22r^ + 17r^) (r - 1) 
96 {r + lf 



where G2, G3, Ma, ALc, ^^2, A'^a)) = [\/3(-2y2 ~4V3y + 4V3yx - 6 + 12 a- - 6x^ + 7r^ - 4r^^/3^/ - 12r^2; + 



Case 6: 

P(X2 G iVpB(Xi)ur^(Xi),Xi G T,) = 

rs9 /-raW /■=i2 r-3W /■1/2 r^e(^)\ A{3^(A,Ni,Qi,Ga,M2,Mc,P2,N2)) 

'si2 -^0 



dydx 



81 - 189 + 561 - 45 r'' - 1894 - 18 r"* + 1912 + 224 - 384 r - 128 



1296 (r + l)''r4 



where yl(£52i(yl,iVi,Qi,G3,M2,Mc,P2,iV2)) = [V3(4r-y2-4 V3y + 12x + 13r-12 + 18r=*a;^ + 12rx--12rx^-8^r2y + 
4 ^/3r y - 24 r'^x + 12 ^r^y a; + 6 y^) j / ^24 rj . 



Case 7: 

P{X2 G iVpB(Xi) U r5:(Xi),Xi G Ts) = 

p /T^C:.)^ rs^o /TsWX yl(^(yl,Afi,Qi,L3,Mc,P2,jV2)) ^^^^,_ 

; Ainy^w 

128 - 1536 r - 302592 r'^ + 11753 r^^ + 346171 - 28416 + 8384 + 69760 / + 220201 r*^ - 135936 - 305664 r V 

186683 r^" - 69120 r^^]/[l944 (r^ + l)'V^] 

where A(^(A,Ari,Qi,L3,Mc,P2,iV2)) = [^(-4 VSr ?/ + 2 VSr^y- 12 a;- 12 r + 8 + 12 r x-6 a;^r^ + 2 -4 a; + 
3 j/2 _ 4 r'^^y _ 12 r^'a + 9 r^'a^ + 4 V^y + 6 r*^y x + 6 + 6 r^a; + 6 + 2 / [l2 r^j . 
Case 8: 

p(x P (X )nr^(x ) X ^T) - ( r r^^\ r r^^^\ a{,^{a,n^,p^,L2,L3,Mc,P2,N2)) 

P(X2GiVpMXi)uri(Xi),XiGr.)-^^ y^^^^J ^^^^^^^^^^ d.rfa- 

(895r** - 2472 + 3363 r-** - 2880 r'^ + 2220 r'' - 1296r=* + 675r2 - 216r- + 27) (-12r- + 7r'' + 3)^ 



7776 (r2 + 1) 



where A{.^{A, Ni, Pi, L2, L3, Mc , P2, N2)) = [y3(4r*j/2 + 8r^^yx + Vlr^x'^ - 4r-^^j/ - Vlr'-^x -Ay^ry - 12r + 
12 r s + 8 + 3 1/^ + 6 ^/3y - 6 VSy a; + 9 - 18 x + 9 a:^) j / [l2 7 



Case 9: 



P{X2 G iv^.£;(^i)uri(Xi),Xi e r,) = 



S3 ria^(^) ^ rsr ^ p ^rs(.)\ ^(^(^, jy^ , ^ ^2 , L3, L4, L5, P2, iV2)) ^^.^^ _ 

- [355328 r^** - 2204160 + 6591792 r^*^ - 13254912 r^'^ + 20639832 r^* - 26417664 r" + 28578916 r^^ - 26760576 r" + 
21960774 - 15877152 + 10180620 - 5753232 r'^ + 2856483 - 1222128 + 438777 - 128304 + 28107 r^- 

3888 r + 243] / [7776 (r^ + l) ^ (2 r V l) r'* j 

where ^(^(^, A^i , Pi , La, L3, 1-4, L5, ^^2, ATz)) = [^(l8 + 4^ry - 18a- - 24r + 12 + 12ra: - 6\/3y + 8r*V3ya: - 
12 r^x + 12 r*x'^ + 9 + 15 + 4 r*?/^ - 4 r^^y - 6 V3y x)] / [l2 r^j . 
Case 10: 

PiX2eNMX^)uriiXi),XieT.)-^J^^ J^^^^^ + J^^ J^^^^J ^^^^^ dydx - 

[8 (288 r^^ - 864 r " + 1486 r - 1896 + 2056 - 1608 r'' + 1189 - 654 + 317 r'' - 132 + 44 - 12 r + 2) (2 r - 1)^ 

(r- if]/ [243 (rVl)'(2r^ + l) 



3 4 

r 



where iVi , Qi, L3, ^4, L5, P2, iV2)) = |^V3(4 \/3r y + 2 ^/3r2y - 8 ^/3y - 12 a: - 24 r + 12 + 12 r a: - 6 a: W 15 - 

12 r'-^x + 9 r-'^a;^ + 6 a;^ + 6 r^a; + 6 r'^VSy x + 2r^y^ - A VSy x + ir^y'^ -A r'^VSy + 14 y^) j / [l2 j . 
Case 11: 



P{X2 € N'pe{Xi) U r!(Xi), Xi G T,) = / / + / / 



A{3^{A,Ni,Qi,G3,M2, N3,N2)) ^^^^_^ 



1536 - 6528 + 133834 r** - 48240 + 95616 - 20736 - 158976 - 200064 + 196680 r** + 7107 r ^° 



15552 r-i 



where ^( iVi, Qi, Gs, M2, iVa, iVa)) = [V3(4ry^ + 12a: + 9r - 12 + 9r-^a:^ + 12 rx- 12 ra;^ - A^Sr^y + A^Sry + 
6 VSr^y x + 3r'^y'^ - 12 r^a: - 4 \/3y)] / ^24 rj . 



Case 12: 



nsi3 /-^W /"'IS /■'■2(=^)\ A(^(A,iVi,Qi,L3,7V3,iV2)) 



P{X2 G N^.e{Xi) U rI(Xi),Xi G Ts) = / / + / / ) -V- "-'7'7;'";'--'---^^ dyrfa: = 

(147 - 504 r'^ + 530 - 336 + 876 r* - 1056 + 896 - 384 r + 64) (-12 r + 7 + 4) ^ 

15552 r6 

where A{g^{A, Ni,Qi, L3, N3, N2)) = [^(4y^ - 8 VSy a; - 24 a; - 24 r + 8 ^y + 12 + 4 ^r^y + 6 r^'^Sj/ a; + 24 r a; - 
4r^V3y + 3r'*y2 - 8y3ry- 12 - 12 r^a: + 9 r*a:^ + 12a;^ + 12r-^a; + 4r^ y^ + 12) j / ^24 r^] . 
Case 13: 

P(X2 G iv;B(Xi)ur5:(Xi),Xi G Ts) = 

.rioW ^ ..13 .ri,(.) ^ .1/2 .ri3W\ ^(.^(4, jV, , , L3 , iV3 , iV2 ) ) ^ 
fill J r2{x) J sn J r2{x) J B13 J r^ix) 

1024 - 12288 r- + 295680 + 1053 - 197140 r** + 626688 - 100864 - 1294848 - 686528 r** + 1282560 + 



114336 - 



30930 r-"'j/ [31104 r''] 



where A{.^{A,Ni,Qi,L3,N3,N2)) = ^V3{'iy'^ - SV^y x - 24x - 24r + sVSy + 12r'^ + 4V3r'^y + er-^VSy x + 24r x - 
4r^V3y + 3r^y^ - 8^/3ry- Ux'^r^ - 12 r^x + 9 r^'x^ + 12a:^ + 12r^a; + 4r^y^ + 12)j/|^24r-^ 



Case 14: 



P{X2 G iv;B(^i)uri(Xi),Xi G T,) = 



- [(r - 1) (1512 r" + 1512 r^'^ - 16740 r^^ + 540 r^* + 84078 - 83538 r^^ - 164835 r^^ + 401085 r^" ~ 487872 r''+ 
535728 - 463124 + 335596 r'^ - 197440 + 64640 - 7936 ~ 1792 + 5632 r - 512) j / ^5184 (2 + l) r"*! 

where A{^{A, Ni,Qi, L3, L4,Q2, N2)) = [^(-6x - 12r + 6r'^ + 6rx + 2 ^r^y - _ 2^ya: + r"* + 5^2 _ 
2 r^x %/3y + 2 r''\/3y x + 2 ^/3r y - 2 r^^/3y - 3 x^r^ - 6 r^a: + 3 r^x^ - 2 ^/3y + 3 a:^ + 6 r^s + 6) j / ^6 r^j . 
Case 15: 

isi3 "''■2(3:) "^1--' IJ'ajJ 

(147 - 612 r'* +980 - 768 + 744r - 288) (-6 + 5r-)^ 
7776 r- 

where A(.0^{A, Ni,Qi, L3, L4,Q2, N2)) = [^/3(4ry2 + 12a: + 9r - 12 + 9r^x'^ + 12rx - 12 r x^ - 4 ^r^y + 4 V3r y + 
6 V3r^y a; + 3 y^ - 12 r^s - 4 \fZy) \ I [24 r] . 



dydx = 



Case 16: 

/■'■lo(^) A{g^{A,Ni,Qi,L3,N3,N2)) 



P(X2 G iV^.£;(Xi)ur5:(Xi),Xi G T,) = / / 



(I3r-** + 52r'^ + lOr*^ - 184 + 60 + 624 - 48r^ - 832r + 448) (-2 + r) (r^ + 2r - 4)^ 

384 (r + 2)^ r2 

where ^(^(^, A/'i, Qi, L3, A^'a, A*'2)) = [\/3(4y^ - 8 \/3y a: - 24 a - 24 r + 8 ^y + 12 + 4 \/3r^y + 6 r''\/3y a: + 24 r a; - 
4r3^/3y + 3r''y2 _ 8^/3ry- 12 - 12 r^a; + 9 r^x^ + 12^^ + 12 r^a; + 4 y^ + 12)]/[24r^] . 
Case 17: 

[(189 1323 r" + 1026 r - 10692 - 14364 + 51732 64664 - 183952 - 153504 + 398080 + 124928 - 

470528 r + 197632) (r - l)j / [5184 (r + 2)-'] 

where A{3^{A, Ni,Qi, L3, N3, N2)) = [\/3(-6x- 12r + 6r^ +6 r a + 2 ^/3r^y-r^ y^ -2 \/3y a; + r'' y^ + 5y^ -2r-^x ^/3y + 

2 r''^/3y a; + 2 ^/3r y - 2 r■^^/3y - 3 x'^r^ - 6 r^a; + 3 r*a;^ - 2 \/3y + 3 a^ + 6 r^a + 6)] / [6 r^j . 

Adding up the P{X2 G Np^{Xi) U r5^(Xi), Xi G T,) values in the 17 possible cases above, and multiplying 
by 6 we get for r G [1, 4/3), 

_ SeOr^- 195r5- 256 + 720r- 846r3- 108r'^ + 47r^ 
~ 108 r2 (r + 2) (r + 1) ' 

The for(^) values for the other intervals can be calculated similarly. 
Derivation of Uorir) in Theorem 13.21 

By symmetry, P{{X2,X3} C N^pe{Xi) U T{iXi)) = 6P{{X2,X3} C N^p^iX^) U T{{Xi), X^ G T,). For 
r G [6/5, a/S — l), there arc 17 cases to consider for calculation of iyor{r) in the OR- underlying version (see also 



Figure 126)) : Case 1: 

4 o 4 

r - — r + 1/9 



81 ' 27 



where A{d^{A, Mi, Mc , M3)) = 1/12 \/3. 
Case 2: 

p({X2,X3} c iv;B(Xi)ur5:(Xi), Xi g r,) = 



- [(r - 1) (119155 r" - 845345 r^" + 2724777 r'' - 5206743 + 6475257 - 5454855 r*^ + 3155193 r'^ - 1249479 r-'*+ 

332181 r-' - 56619 + 5589 r - 243)] / [25920 r"^] 

where ^(^(yl, Ml, L2,L3,Mc,M3)) = — ^ '-. 

Case 3: 

' -5 /-SB pW\ A(^(A,G2,G3,M2,Mc,M3))2 



P({X2,Xs} C iV;B(^i) Un(Xi), Xi G T,) = / / + / / 



215 - 136 r'^ - 56 + 172 - 55 r* - 60 + 66 - 24 r + 3) (r - 1)" 

2880 rio 



where Ga, G3, M2, Mc, M3)) = ^- — — '-. 

Case 4: 

P({X2,X3} C iV;E(Xi)ur^(Xi), Xi G Ts) = 

r^a,„W /-sa r-5(-)\ A(^(A,Mi,L2,L3,L4,L5,M3))' 



dydx ■ 



X) J A{T{ys)r 

37072 - 195072 + 453120 r'^ - 589248 + 460728 / - 217728 + 60480 - 9072 r + 567) 



4r-3 + Vs)^ (^4r-3-\/3)^]/[l866240r^°] 



u r r r r VsU V3r y + 9 -24 v + 12 r x + 15 -6 Vay-6 VSy x + 18-18 x + 9 x^) 

where 4(^(4, Mi, L2, L3, L4, Ls, M3)) = — ^ ^. 

Case 5: 

P({X2,X3} C N;,E{Xi)url{Xi), Xi G Ts) = 

^6 /■'■^W /-.g /ttWX A(^(A,G2,G3,M2,Mc,P2,iV2))' J 

+ / / I MTiMf '"''"^ 



■5(3;) J SQ Jo 



(3 - 12 r - 15 + 84 + 18 r'* - 232 + 130 + 504 r'' - 108 r** - 288 + 623 + 920 r" + 373 r^^) (r - 1)^ 

2880 r-io (r- + 1)^ 

where G2, G3, M2, Mc, P2, Af2)) = [\/3(-2y^ - 4 x/3y + 4 x/3y x - 6 + 12 a: - 62;^ + 7 r'^ - Ar^y/^y - 12r^x + 



! r^V^y x + 12 r'^x'^ + 4 r^y^) j / ^24 ? 



Case 6: 

p({X2,X3} c iv^B(Xi)ur5:(Xi), Xi g t^) = 



rso rr^i^) ^ rsri rr^i^) ^ rre(x)\ A{-S^{A, Ni, Qi, G3, M2, Mg , P2, N2)f ^ ^ 

Js, Jrrix) Jo Js,2 Jo J A(T{ys))^ ^ ^ 

19683 r^^ - 59049 + 83106 + 167670 r^^ - 211626 r" + 344466 - 142614 - 2573586 - 128853 r''+ 
3465675 r'^ + 1103824 - 1473304 - 730880 + 107776 + 158720 r + 31744j / [l049760 (r + l)"^ r'*] 



where A{^{A,Ni,Qi,G:i,M2,Mc,P2,N2)) = ^V3{4ry^ + 12 x + 13r + 12 rx-4V3y -12 + 4 V^ry-SV^r^y + 18 x'^r^ 
12 r + 6 r^y'^ - 24 r'^x + 12 VSr^y x)] / ^24 rj . 
Case 7: 

p({X2,X3} c N;,E{Xi)url{Xi), Xj e t,) = 

I l^s, ^r,(.) ^ 1^,0 ^r,i.) ^ ^.1. iVi, Qi, L3, Mc, Pa, iVa))' ^^^^ _ 



- 1^6144 - 1 10592 r - 310846464 + 2127553557 r + 570050560 - 5031936 r'^ +936960 + 19526656 + 147203072 r 
7627473 + 1419072042 - 762467328 + 288811029 r^^ - 68327424 - 59166720 - 923627520 r^ + 
1340817105 r"^ - 1765251072 r^^ - 2350015488 r^^ + 2339575338 r^* - 2016377856 r^^j / ^262440 (r^ + l)%"^j 



r X 



where A{.0»{A, Ni,Qi, L3, Mc , P2, N2)) = [^(-4 V3ry + 2 V3r^y-6a;V^- 12x- 12r- 12r^; 
6 + 6 r*^/3^/ x + 2 r^y'^ - 4 \/3?/ 1 + 3 r^'y^ - 4 r^\/3y + 4 \/3y + 2 + 6 r^s: + 6)] / [l2 r^j . 
Case 8: 

p{{X2,X3} c iv;B(Xi)ur5:(Xi), Xi g t,) = 

[(426497 r^^ - 2443992 r^^ + 6726107 - 11753232 r^^ + 15220771 r^^ - 16367448 + 15754449 7-^" - 13773024 r^- 
10839672 - 7552440 + 4592889 - 2374272 + 1018899 - 344088 + 81891 - 11664 r + 729) 

(_12r + 7r2 + 3)']/ [699840 (r^ + l)' r"^] 

where A( J«(y4, A^i, Pi, L2, L3, Mc, P2, A'^a)) = [\/3(-4r^^y - 12r^2; + Sr-'^ya; + 12r'*a;^ + 4r'^y'^ - 4^ry - 12r 
12 r 2: + 3 + 6 \/3j/ - 6 \/3j/ a: + 8 + 9 - 18 X + 9 a:^)] / [12 r^] . 
Case 9: 

p({X2,X3} c iv;B(Xi)ur^(Xi), Xj e r,) = 

/ rss ria^(x) ^ rsr rl^^(x) ^ p rrs(x)\ A{^{A,Nr,Pj,L2,Ls,L4,U,P2,N2)f 

; A{T{y3)y 

- [l5309 - 367416 r + 60475010560 r^^ + 437704472832 r^*^ + 1444872192 r^" - 13250101248 r^^ - 185909870592 

4148739 - 2027754648576 r^^ + 1397612375040 r^* + 20429177589 - 677278256112 - 49656902904 
159963012 - 30005640 - 681714144 ~ 7515142416 r - 3097406755584 r^^ - 2609245249920 r^^H 
3051035360256 r^* - 3315184235136 r^^ + 3337272236928 + 2631941507968 r^^ + 2435971806 r' + 
109069315047 r^" - 218273842152 r" + 400534503738 r^^ + 1059615993384 r^* - 1538314485120 r^^ + 

2076627064432 r^^ - 845838600192 r^^] j [l399680 (r^ + l) * (2 r- + l) r 

where iVi, Pi, L2, 1-3, ^4, L5, P2, iV2)) = [\/3(l8 - 18a; - 24r - 12r^x + 12r'*a;^ + 12r^ + 12r-a; + 4 ^3; 

4 r^^y + 4 r'*y^ - 6 x + 8 r^^VSy a; + 9 + 15 j/^ - 6 ^y)] / [l2 r^] . 
Case 10: 

P({X2,X3} civ;B(Xi)ur^(Xi), Xi G r,) = 

/ ..e A(^(A,iVi,Qi L3 L. L.,P2,iV2))- 

[64 (12 -U4r + 924 - 683328 r^^ + 112976 r^" + 757211 - 10554918 r^^ - 1513230 + 16242 / - 4320 - 51372 
344988 r'^ - 4867848 r^^ - 18583080 r" + 16493828 r^^ - 12883116 r ^ ^8668124 r^" + 2177536 r^^ + 141366 r** + 2774371 
4692510 r"+7331714r^^ + 14002613r"-16948218r-^^ + 18708475r^'^) (r - 1)^ (2r - l)^jy/[32805 (r^ + l)^ (2r^ + l)% 



27 , 

r + 



r' + 



>ry 



where A{^{A, Ni, Qi, L3, L4, L5, P2, = [^(2 VSr'y + 15-6 x^r'^ - 12 2: - 24 r - 12 r^x + 9 r^x^ + 12 + 12 r x - 
8 VSy + 6x^ + 6 r^VSy x + Uy'^ -4 VSy x + 2 r'^y^ - 4 r■^^/3y + 3 r^y^ + 6 r^s + 4 ^/3r y) j ^ |^12 j . 
Case 11: 

p({X2,X3} c iv^E(Xi)ur5:(Xi), e r,) = 

- [-253952 + 1529856 r' + 601574256 - 385780320 - 776518272 + 7803648 r"* - 70917120 - 396524160 
209710080 r*^ +869661288 r^" -845940960 r"+668092108 r"+147067614 r^^-32610600 r^''^+3173067 r^*^] / [8398080 r*^] 

where iVi, Qi, G3, Ma, iVa, iV2)) = [^(4 r + 12 x + 4 1/ + 9 r - 4 + 12 r x - 12 + 9 + 6 ^r^y a; - 

12 r - 4 ^r'^y - 12 r'^x + 3 r^y^)] j [24 rj . 
Case 12: 

P({X2,X3} C iV;B(Xi)un(Xi), Xi G T,) = 

[(64827 r^^ - 444528 r^'^ + 1223334 r^^ - 1793232 r^^ + 1839416 r^^ - 2003712 r" + 2286224 - 2421504 + 3095088 r*- 
4428288 r + 5889152 - 6093312 + 4557056 - 2340864 + 774144 - 147456 r + 12288) 

(- 12 r + 7 r V 4) ^ j / [8398080 r 

where A{&>{A, Ni, Qi, L3, N3, N2)) ^ [V3(-12 xV^ - 24 a; - 24 r - 12 r^x + 9 r^x'^ + 4 y^ - 8 ^/3r y + 6 r'^VSy x + 8 ^/3y + 
12 + 24 r X + 12 x^ - 8 VSy x + 4 r'^y'^ - 4 V3y + 3 r'*y^ + 4 V3r^y + 12 r'^x + 12)] / [24 r^j . 
Case 13: 

P({X2,X3} CiVj.B(Xi)un(Xi), Xi G T,) = 



Xi G Ts) = 

/ /-^- ^ rsi, p2W^ /-i/a p2(:.)\ ^(^(4, iVi, Qi, L3, iVa, iV2))' ^^.^^ _ 

[l96608 - 3538944 r + 8927944704 r'^ - 1883996112 r^^ - 9492593152 r** - 146866176 + 29196288 + 220250112 7-"- 
4486594560 r*^ + 213597 - 259250904 r^'^ + 69124752 - 10683306 r^^ + 864387072 + 5220357120 r^- 

1081136256 + 602097408 r" + 2223664128 r^^ - 1509638512 r^"* + 716568768 r^^j / [l6796160 r**] 

where y4( ^^'(.4, A^i, Qi, L3, A^s, A''2)) = [^/3(-12 xV^ - 24x - 24r - 12 r=^x + 9r'*x2 +4y2 - 8 ^/3ry + 6 r^^Syx + 8 V3y + 
12 + 24 r- X + 12 x^ - 8 VSy x + 4 r'^y'^ - 4 r^VSy + 3 r^y'^ + 4 V3r^y + 12 r^x + 12)] / [24 r^j . 
Case 14: 

p({X2,X3} c iv;E(Xi)ur^(Xi), Xi G r,) = 

/ ..11 rta^{x) ^ ..10 .ri2(a.)^ ..14 /-l^WX A( ^ ( A, iVl , Ql , L3 , L4 , Q2 , iV2 ) ) ' ^^^^ _ 

- [(r - 1) (-16384 + 278528 r + 215136 r^** + 40176 r^'' + 215136 r^^ - 3381264 r^^ - 2301952 - 99212040 r^^- 

25050384 r^^ - 312101312 - 7215869272 r^^ - 147586784 - 42770432 r"* + 12591104 + 114049024 + 345810944 r^+ 

55914462 r^^ - 2082969096 r^'' + 43443459 r^® + 826941555 r^^ - 641846754 + 209930616 r^^ - 232963072 
1311322268 - 3191747236 r" + 5434516904 r^^ + 7756861008 - 6865898928 r^^ + 4727296416 r^* + 26115696 r^") j y/ 

[466560 (2r^ + l)%-^j 



where A{.^{A,Ni,Qi, Ls, L4,Q2, Ni)) = ^V3{-Sx^r^-6x-12r~6r-^x + 3r'^x'^+2V3ry + 6r'^ + 6rx + 3x'^-2V3y- 
2 VSr'^y x + 2 r^V^y x + 2 VSr^y - r'^y'^ + 5 - 2 r^VSy + r'^y'^ - 2 ^/3^/ a; + 6 + 6 r^x) j ^ r^j . 
Case 15: 

[(63855 - 498960 + 1650060 - 3036960 r'^ + 3703292 r*^ - 3657696 + 3268368 - 2419200 + 1550448 r^- 

725760 r + 155520) (-6 + 5 r)^j / [4199040 r^] 

where iVi, Qi, G3, M2, iVs, iV2)) = [\/3(4ry^ + 12 x + 4 VSr 1/ + 9 r - 4 ^/3^/ + 12 r 2; 12 + 9 x\-^ + 6 ^r^y a; - 

12 r x^ - 4 y/2>r'^y - 12 r'^x + 3 r^y^)] / [24 rj . 
Case 16: 

- [ (293 r ^'^ + 2344 r + 4662 r " - 9088 r - 32320 r + 42976 r " + 175408 r - 1 19680 - 544144 r** + 372352 
1216512 r'' - 882688 - 1564672 r'^ + 1373184 + 924672 - 1314816 r + 380928) 

(-2 + r) (r-V 2r- - 4)^j / [23040 (r + 2)''^ r"] 

where A{^{A, Ni,Qi, L-j„ N-a, N2)) = [V3(-12 xV^ - 24x - 24r - 12r^x + 9r''x^ +4j/^ - 8 \/3r y + 6r''^y x + 8 + 
12 + 24 r X + 12 x^ - 8 VSy x + 4 r^y^ - 4 V3y + 3 r'*y^ + 4 V3r^y + 12 r^x + 12)] / [24 r^j . 
Case 17: 

p({X2,X3} c iv;B(Xi)ur5:(Xi), G r,) = 

A(^(A,jVi,Ql,Z,3,L4,Q2,iV2))' 



S14 ptam{x) fl/2 flam{x)^ 

+ 



[(6723 r^" + 73953 r^^ + 213678 r - 433512 - 2873232 r^® + 627264 r^^ + 20218896 + 5675184 r^^ - 97577924 r^^- 



39916108 r" +343932568 r ^° + 108508576 - 906967296 r** - 96480192 1702951296 - 293251072 r"^ - 1994987520 

981590016 + 1118830592 - 1135919104 r + 287604736) (r - 1)] / [466560 r"* (r + 2)^] 

where Ni,Qi, L-i, Ni, N2)) = [V3(-3 x^r^ - 6 x - 12 r - 6 r^x + 3 r^'x^ + 2 \/3r y + 6r-^ + 6rx + 3x^-2 v^y - 

2 V3r^y x + 2 r'*\/3y x + 2 \/3r^y - r^y^ + 5 y^ - 2 r^Viy + r^y"^ - 2 VSy x + 6 + 6 r^x)] / [g r^j . 

Adding up the P{{X2, X3} C iVp£;(Xi) U r5^(Xi), Xi e T^) values in the 17 possible cases above, and 
multiplying by 6 we get, for r e [6/5, \/5 — l), 

,y„,{r) = 413208 r+3070468 r2-74801558 /+75243552 ri3-4883958 r^+14541630 r^+28880-11254002 

3667716 r^ + 64360782r^+ 13122 -3300900 r^^+lSeOUr^^- 175011 ri^ + 62825r2" + 1458 r22- 19812000?'*^+ 
99831906 r^" - 120628524 r" + 33155180 r^^ - 67685050 r^"^ + 5055135 r^^ + 11053023 r^^ ^ 116640 r*^ (r^ + l) 



The i'or(^) values for the other intervals can be calculated similarly. 



l)(r + 2)' (r + lf 



Appendix 5: The Asymptotic Means of Relative Edge Density Under 
Segregation and Association Alternatives 

Let Mfnd('"'^) ^^^^ M£md(^'^) ^^^^ mcans of relative edge density for the AND-underlying graph under the 
segregation and association alternatives. Define ^) ^^'^ A'or(''' ^) similarly. Derivation of /ifnd('^' ^) involves 



detailed geometric calculations and partitioning of the space of (r, e, x) for r G [1, oo), e G [O, V3/3) , and x G Tg. 
See Appendix 6 for the derivation of n(r,e) at a demonstrative interval. 



/if^jj(r, e) Under Segregation Alternatives 



Under segregation, we compute Mfnd(^' ^) ^^'^ Mot(''i ^) explicitly. For e G [O, \/3/8) , /ifnd(^' ^) = Si^i ^^f"'^('', G 
li) where 

, s (r - 1) (5 + 288 r^e* + 1 152 r^e* - 148 + 1440 r^e'* + 245 - 178 + 576 r^e^ - 232 r + 128) 

^37, (r, e) = = 5 ^ 

54r2(2e-l)2(2e+l)2(r + 2)(r + l) 



^T"^ (r, e) = - [ll52 r^e"* + 101 r'' +3456 r^e* - 801 + 1302 r'^ + 1152 r^e"* - 732 - 3456 r^e* - 536 r - 2304 re* + 672j / 

[2I6 (r + 2) r (16£* - + l) (r + l)j 



zus'"^ (r, e) = - 1^-3 r V 128 r*e* + 384 r^e" + 39 r V 128 r'^e'^ - 90 r*^ - 444 - 384 r^e* + 1344 - 256 r^e" - 792 

864 + 1104 r - 288j / [24 r* (16 e'* - 8 + l) (r + 1) (r + 2)] 



with the corresponding intervals Xi 



16 r^e'' + 16 r®e* - 3 - 16 r^e" - 3 - 16 r*e'^ + 9 + 9 - 18 r + 6 
3 (r + l)r-4 (4e2 - 1)^ 

1, 4/3) , I2 = [4/3, 3/2) , I3 = [3/2, 2) , and X4 



2,00 



For e G 



0,V3/8), ^liir,e) = Eti ^) e X,) where 



zji (r, e) = 1 47 r - 195 r + 576 r- e"* - 288 re + 860 r - 846 r- + 1728 re* - 864 r - 108 - 576 r + 1152 re + 

720 r - 256] / [1O8 (16 re'' - 8 re^ + r - 16 + 2 + 32 e"*) (r + 1)] 



175 r^ - 579 r" + 1450 r^ + 1152 r^e" - 576 r^e^ + 3456 r^e* - 1728 r^e^ - 732 r^ + 2304 re* - 536 r - 1 152 re^ + 



672 



:] / [2I6 (r + 2) r (2 e - 1)' (2 e + 1)^ (r + 1)] 



zu3'"^{r, e) = - [27 r* - 63 r'^ - 270 r** + 1728 r^e^ - 384 r^'e* + 1024 e^ V3r^ - 1152 r^e* + 576 r^e^ + 756 r^ + 1536 r*e^\/3- 
2376 r* - 6912 r*e^ - 2560 VSe^r'^ + 2304 r^e* + 2736 r^ + 1 152 r^e^ + 1296 r^ - 3072 r^e^V3 + 1536 r^e* + 6912 



3r^e* + 6912r^e^-3312 
864] / [72 r" (r + 1) (16 re* - 8 re^ + r - 16 e^ + 2 + 32 e' 



r^ 

,4\" 



t:<7r'^(r, e) = - [-18 - 48 r^e* - 48 r*e* + 72 r*e^ - 144 r^e^ - 9 r* - 32 r^'e* - 144r='e^ + 72 r^e^ - 9 r^ - 32 r^e* + 54 rH 

64 r^e^ V3 + 64 ^eV^ j / [9 r* (4 e^ - l) ^ (r + 1)] 

with the corresponding intervals Xi are same as before. 



/i^jjj(r, e) Under Association Alternatives 

Under association, wc compute fi^ndir,e) and fii{r,e) explicitly. For e G [O, (7 - 3 \/T5) /12 w .042), 
^^tAr, e) = Eti e) G where 

<;r'^(r,e) = - [-128+768 r-''\/3£^+360 r+8640 £"+5760 e^+393 r*-54 r2+6912 r^e^+S r'^-lSS r^-423 r^-4608 r*^e^ + 
6912 ^r^£^+1728 eV-3072 ^£^-7776 r-^£''-864 r''£''-2592 r'^£''-18144 £''r^+12960 £''r--576 r®£^-3456 r-^£^ + 1728 r'^e^- 



7776 r'^e^ - 12096 r^£^ 



]/[6 (ys + 6ey (^-6 £ + {r + 2) {r + l)j 



'i^'^ir, £) = [-672 r+20736 £*+13824£^-1302 r-''+536 r^-lOl rV801 r^+732 r^-3072 r'* ^£^ + 18432 r*^£-'-9216 VSre^ 
19968 \/3r^£^+4608 ^/3r^£^+31104 ^£''+4608 r^£^-17280 ^£^58752 £*r^~6912 £^+3456 r-'*£'*+10368 r^£''+72576 £V^ + 



31104 £'*r + 2304 r'^e^ + 17280 r^£^ - 6912 r^£^j / [2I6 {r + 2) (r + 1) (-1 + 12 £^) 



?r'*(r,£) = [9(r-*-13r-''+30r'^-192r''£^ + 1152r''£'' + 148r'^+3456r^£*-576r'^£^-448/+2688r*£*-128r''£^ + 1152£*r^ + 



264 + 768 r^£^ + 512 r^e^ + 768 £''r^ + 288 - 368 r + 96)] ^ ^8r^ (^-6e + Vtj (^Vs + 6e^ [r + 1) {r + 2)] 

and/ X 9(r^ +6r- + r'' - 3r^ - 3r^ - 2 + 144r^£'' + 144 r*£'' + 48 £'^r^ + 48 £''r2 - 24r^£^ - 24r''£^ + 32 r-^£^ + 32r^e^) 

?4 {r,£) = — 2 2 

r4(r + l) (-y3 + 6£) (y3 + 6£) 

with the corresponding intervals are same as before. 

For £ e [0, (7 x/3 - 3 VTE) /12 « .042), /x^,.(r, e) = ^^i ^r(^> ^) e X,) where 

?r(r,£) = [-256+720 r-13824£*-9216 £^+860 /- 108 r^+47r'^ -195 r''-846r^ + 12096r"*£'*+6912r^£^ + 1152/£^ + 
31 104 e*r'^ - 6144 VSe^ + 3072 ^/3£^ - 6144 r^VSs^ + 13824 VSr'^e^ + 4608 \/3r^£''^ + 13824 e'^r- 10368 r-^£* + 57024 e*r^ - 
20736£''r- 2304r'^£^ - 17280 r^£^ - 3456 r®£^] / [l2 (r + 2) (^-6£ + \/3)^ (Vs + Gey r'^ (^ + 1)] 



^"(r, £) = - [-672+579 r''-1450 r^+536 r+20736 r''£*+32832 r^£^-114048 £*r^-7488 r^£^+8064 £^-175 r^+6912 r^e^+ 
4608 r^£^ - 24192 e^r'^ - 76032 e^r + 12288 VSr^e^ - 9216 r-*\/3£^ + 4608 VSr'^e'^ + 732 - 6144 VSr^e^ - 9216 VSe^- 

19968 \/3r£^ - 27648 £^j / [2I6 r (r + 2) (r + 1) (-1 + 12 £^)^] 

?3''('-,£) = - [9(96+384 r*£^ + 192r-®£^-2304/£''-30r''-1152r-^£*+84r''+576r''£^+3r**-7r'^-368r+304r^ + 144r^- 

3456 /^e* - 264 r"*)] / [8 r* (r + 2) (r + 1) (Vs + 6 £) ^ (^-6 £ + \/3) ^] 

or, , 9(-6 r + r* + 7-^ + 2+ 144 r5£'* + 144 r''£''- 24 r5£2- 24 r''£2) 

?4 ('", £) = 2 2 

r4(r + l) (-6£ + ^ (v^ + 6£) 

with the corresponding intervals X^ are same as before. 

Appendix 6: Derivation of /ifnd(^'^) A^frl^'^) 

We demonstrate the derivation of i^sif, e) for segregation with e e [O, v^/S) and among the intervals of r that 
do not vanish as £ 0. So the resultant expressions can be used in PAE analysis. 



Derivation of /Ltfj^j(r, e) 

By symmetry, 

^li,^{r,e)^P{X2 e N^pE{X,,e)nriiX,,e)) =6P{X2 e N^{X,,e)nTl{X,,e),X, e T,\T(y,e)). 

Let q{yi, x) be the line parallel to and crossing r(3^3) such that d{yi, q{yi,x)) = £ for i = 1, 2, 3. Furthermore, 
let T, := T{y3) \ U^j^,T{y„e). Then q{y,x) = 2e - V^x, q{y2,x) = Vix - V3 + 2e, and g(y3,x) = \/3/2 - e. 
Now, let 



^1 -'Z(y,2:)nyy^= (2£/V3,0j , V2 = q{y2,x) = (l-2e/^/3,0j , 

^3 = '7(y2,a:)nyiyi= (l-e/V3,e), 1^ = (/(Vs, a:) n yiW = (l/2 + e/V3, V3/2 - , 

^5 = (Z(y3, a:) n yyi = (l/2 - e/VS, V3/2 - e) , l/g = <z(y, a;) n yy? = (e/VS, e) . 
See Figure m 

The points Gi, for i = 1, 2, . . . , 6, Pi, for i — 1,2, Li, for i = 1, 2, . . . , 6, Ni, for i = 1,2, 3, Qi, for i = 1, 2 
and the lines ri{x), for i = 1, 2, . . . , 11 are as in Appendix 3. 

so = -2r/3 + 1, .si = -r + 3/2, S2 = 3/(8r), .S3 = 1 - r/2, .S4 = 2(2?^+!) ' «5 = = l/(2r), 

= 1/(2 r), ss = --2,£_W+2.^ ^ --^-ff;t3zi, ^lo = l/(r + l), = .12 = 

si3 = '"'6(r'2+;f'^ si4 = 3r/8, si5 = r - rV8 - 1/2 

liix) = 1/3 V3 (-3 a; + 2 £ V3) , ^2 (a;) = -1/3 " ' , ^ " ^("^'-1) , ^ (x) = 1/3 ^/3 (-3 x + 2 e V3r) 
gi ^ 1/2 eVS, 92 - 2/3 eV3, qs = -1/4 ^^±f^, 94 = 3/4 r-i, = l/2e V3r, and qs = 2/3eV3r 

Then r(y,£) = r(y,0i,Q6), T(y2,e) = T(g2,y2,Q3), and r(y3,e) - r(Q4,Q5,y3), and for e e [0,V3/4), 
is the hexagon with vertices, Qi, i = 1, . . . ,6. So we have A(Te) = — £^\/3 + \/3/4. 

For r £ [1, 4/3), since e small enough that ^2(2;) n = 0, then iV(a::, e) C for all x G Tg \ T(y, e). There 
are 14 cases to consider for the AND-underlying version: Case 1: 

P{X2 G iV^.B(Xi,e)nn(Xi,e),Xi G \r(y,e)) = 



4e* (-3r2 + 2 + r-«) 



where ^(.^^^(yi, iVi, iVz, Ve)) 



9 (4e2 _ ly 



9(4e2_i)^ 

Case 2: 

P(X2 GiV];B(Xi,e)nrI(Xi,e),Xi GT, \r(y,e)) = 

256 e"*r^^ - 256 e*r^ - 9 + 9 



/ /-^a /-fa^W^ /■ = 2 /-^a^W^ /-SB ^TsCa^jX ^ ( ^ (Gl , A^l , , ) ) ^^^^ 

\J97 J ti(x) J qs Jo J 32 Jo 



where 4(^(Gi, iVi, iV„ Ge)) = ^ (-^^^+1/4 v^H.+TS.) ^(.^-i) 



576 rS (4e2 - ly 



Case 3: 



P(x.Giv;;.(x.,.)nri(x.,.),x.Gr.\T(y,.))^ f f p+ P f'^') ^(^(^.,^^2.M3,g.)) ^^^^ ^ 

yjsil Jrs(x) Jaa Jo J ^Uej 

9r^ - 13r-* - Mr'' + 30r'^ - 22r^ + 22r* - 6 - 10 + r + 3 
96 (4£2_ 1)2^6 (r + 1)^ 

1 Af,T^r^ Ar r, »r 2 (- ^^+1/4 y^) = (_ 12 r3y + 2 r* ^3 + 12 r^y a:- 12 v^a: + 6 r^* x/Sx^ +3 v/Sr^ +2 v^y Vl2 y a;+6 VSa^^ ) 

where A(.;3^(Gi,iVi,P2,Af3,G6)) = ^ 9^2(4^2_i)^ ■ 



Case 4: 



P(X2 eiV;B(Xi,e)nrl(Xi,e),Xi GT,\r(y,e)) 



[243+7022682 r^^-1296 r+36612 r*-952704 r-^'^+137472 r^®-578976 r-''+7057828 r"-5116608 r-^"+2792712 r^^-7725792 r"- 
5484816 + 3631995 r^" - 2213712 + 1213271 + 3051 - 11664 - 101952 + 292518 r^^ j [l5552 (r^ + l) ^ 

(2r^ + l)%'' (4e^-l)'] 

9 ( 4 e-^ — 1 J 

Case 5: 

313 /■'"9(a;) /■S12 rro{x) 



P{X, e iV];,(Xi,e)nrl(X„e),X, g TAT(y,^)) = I / / + / / I ^(■nGi,M,,L, Q,,P,,M,,Ge)) ^y^^ ^ 

\J B4, Jrs(x) Jai3 J r2(x) J ^Uej 

- [4(400 r^^-2832 r-"+8012 r-"-13608 r^^+16350 -14292 r^°+8677 r^-2442 r*-1963 r^+3288 r-'^-2751 r'^+1710 r''-743 

288 - 118 r + 24) j / [243 (2 r'^ + lf (r^ + l)^ (16 - 8 + l)] 

where yl(^(Gi, Mi, La, Qi, P2, M3, Ge)) = (-e^\/3 + 1/4 (--9 + 42^y3; - 45x-2 + 36x - 15y^ + 21 r'^ y'^ + 

2r*2/*-12r'*x^i/^ + 12r*y^a: + 18a;^-12y3y + 42y^a;-24r^ y^-&r'^yy/?,x + A^y^x + 12yx^y/?, + bAr'^x^ + 4,r*y/lf + 
12 r'^y y/?,x - 12 r^VSi^y + 18 r^x'^ + 6 r* - 36 r*s:=* + 18 r^x'^ - 18 r'^^/?,x^y + 12 r^VSaj^y + 12 r'^x^^/2,y - 4 r^^y^^i + 
12r^y^- 45r2s^ + 9r2 - 12 r=^y ^/3 - 4 V^y^ - l^r'^y'^x + Ur^y^x - A2y x^VS + Gr^VSy^ + 2r'^y^ ~ 2Ay^x^ - 
18r2a^* -36 r^a;=^ -36 r^'xH- 72 r-^a:^ -2 [3 r-2 (-^y- 3 + 3x) {-y ~ + V^x) (4e^ - l)^j. 

Case 6: 

P{X2 G N^PE{Xi,e)nrl{Xi,e),Xi € T, \T(y,e)) = 

/ ..10 rr,(x) ^ rs, rr.ix) ^ .1/2 rr,(.)\ yl( ^ (g^ , ^2 , Ql , ^2 , M3 , Ge) ) ^^^^ ^ 

324 r" - 1620 - 618 + 4626 r** + 990 - 2454 + 2703 - 5571 - 3827 + 1455 + 3072 r + 1024 

7776 r-6 (r + 1)^ (16 _ 8e2 + 1) 

where ^(^(Gi, G2, Qi, P2, M3, Ge)) = - [2 {-e^Vs + 1/4 Vs)'^ (-9 V3r'^ - 24 ^/3^ s - 21 r^y - 8 r'^VSy'^ + 24 r'^VSx^ - 

3 VSa; + 24 r ?/ + 24 y 2; - 24 VSx^ - 8 VSy^ - G VS + IS VSx - 4 y^ + 12 VSr + 12 VSr x'^ + 4 VSy^r ~ 18y + 12 r'^x^y - 
24 v^r^a:^ + 8 V^r^ i/^ + 12 r'^x^y/S - 24 y r a: - 4 j/^ + 24 r^y - 4 y^^Sa; - 12 a;^y + 12 r^a:^^ - 12 r^a;^\/| + 4 y^ VSx - 

Ar^VSy'^ - 24r''ya; + 24r^^/3a; - 12r*y3x^ - 4 r--* y^ + 4r^ y^VSa; + 12 a;^^)j / [3 (-V3y- 3 + 3a;) (4e^ - l)^j. 
Case 7: 

P(X2 GiV;B(Xi,e)nn(Xi,e),Xi GT, \r(y,e)) = 

/ ..14 ..12 .ri2(.)^ ..1. .ri2(.)\ 4(^(Gi,Mi,L2,Ql,Q2,L5,Af3,G6)) ^^^^_ 

[l080r"-18900r^^ + 17280r"+65934r-^^-112320r^^ + 152361r"-367200r^''+491051r^-409872r**+282224r'^-60864 
86886 + 70560 - 44672 + 30720 - 16640 r + 6144j / [l0368 (2 r V l) ^ (16 e"* - 8 + l)] 

where A(^(Gi, Mi, L2, Qi, O2, L5, M3, Ge)) = [4 (-e^^ + 1/4 ^/3)'' (-18 + 24 x/3y x - 54a:2 + 54a: - ey^ + 21 y^ + 

r'*y* - 6r^a:^y^ + 6r''y^x - 4y^ + 18a:^ - 6^/3y + 42y^x - 24 r^y^ - 18 r^y^Ss + 12^y^a: + 12 y a:^^/3 + 72 r^a:^ + 
2 r*y3 y^ + 6 r*y ^x - 6 r* VSx^y + 9 r^x^ + 3 r"* y^ - 18 r*a;^ + 9 r^x* + 12 v/3a;^y + 12 r^x^ y^ + 18 r^y V3 + 18 r^x - 
81 r^x'^ + 9 - 12 r^y ^ - 4 r^VSy^ - 24 y^a; + 12 y^a; - 30 y a;^^ - 2 y* - 36 y^x^ - 18 r'^x^ - 36 r^a;'' - 36 r^a; + 
72r^a;2 - 6 ^/3y^)] / [3 (VSy + 3- 3a;) (-y - \/3 + VSa;) (4 5^ - l)^]- 



Case 8: 



lag J re(x) 

81 r^^ + 2048 + 384 r* ~ 810 r"^ + 1296 - 3072 + 96 r*^ 



15552 r6 (16 £4 -8 £2 + 1) 



where A(^(Gi, 6*2, Qi,iV3,Mc, Ma, Gs)) = - [2 (-£2^/3 + 1/4^3) (-5^r2-24^ra:-17r2y-8r2^2/2^24r-2V3a:2 
7 ^/3a; + 24 r y + 24 1/ x - 24 ^a;^ _g^j^2_g^_|_-]^g _ 4 _|_ -j^g + 12 x'^ + ^VSy^r - 18y + 3 r^x^y - 
12 v^r^a:^ + 4 + 3 r''a;^x/3 - 24 y r x- - 4 + 12 r^y - r* y^VSx - 12 x^y + 12 r^x'^y - 12 r^x^ + 4 V^x - 

r^Vs - 6 r-'y X + 12 r'-^VSx - 3 r^'^Sx-^ - r^y^ + 4 y^^Sa; + 12 x^^S)] / [s (-^Sy - 3 + 3 x) (4 £2 - l) 



Case 9: 



"13 fr2{x) rail rrs{x) 



P{X2 G N-pE(Xi,e)nrUXi,s),X^ G TAT(y,£)) = / / + / / ' V,^ ^' "^^ dyrfa: 



r5(a;) ^313 Jr5(a;) 



^(i3^(Gi,Mi,Pi,P2,M3,G6)) 



A{T,) 



- [243+8673 r-^^-1296r+23571r-''-119712/-61488r-^^+169716r^''-246672r^+216121r** + 1404r^-3888r^-35424r^ + 



48816 r' 



y[7776r« (4£2-l)'(r2 + l)'] 



where A(^(Gi,Afi, Pi, P2,M3,G6)) = ^ ^, ^ 

9 r-^ ( 4 e-^ — 1 J 

Case 10: 

P(X2 G iv^.B(Xi,£) nr'i(A'i,£),Xi g t, \r(y,£)) = 

" ^'■^"'"^ , yl(i3^(Gi,Mi,L2,Qi,A^3,L4,L5,M3,G6)) , , 

324 r" - 6949 + 7248 + 26896 - 24960 + 2160 - 259200 + 645760 - 552960 r2 + 155648 r + 6144 

31104 r3 (16£'i - 8£2 + 1) 

where A(^(Gi, Mi, 1,2, Qi, Ar3, -L4, L5, Mg, Ge)) = [2 (-£2^+ 1/4^)^ (-72 - 24^yx - 144x2 _ ^44 3. ^ _,. ^gg a; + 

24 y2 + 72 r + 30 y'^ + r*y*-& r^x'^y'^ + 6 y'^x - 24 j/* + 36 x"^ + 12 ^/3y + 84 y'^x - 24 y2 12 r'^y ^/3x + 56 \/3 y'^x + 
24t/x^^/3 + 108r2x^ + 2r''^2/^+6r^y ^/3x-6r^^/3x2y + 9r''x2+3r■* j/2 _ ig ^4^3 ^4^4 _ gg ^2^^2^^ -^2 r^^x2y + 
24 r^x'^y^y - 8 r^VSy^x -72ry'^ + 96r y'^x + 72 r2x - 126 r2x2 _ 18 r-2 + 72 r x2 - 12 r^y ^3 - 4 r^^j/^ + 48 r y VSx - 
48 r y x2 - 36 r2 y2x + 12 y2x - 12 y x2^ + 12 r2^y3 + 4 r2 y"* - 120 y^x2 - 36 r2x* - 36 r^x^ - 36 r^x + 72 r^x^ - 
28V3y^ - 16r^y^)j/[3r2 (V3y + 3-3x) {-y - Vs + VSx) {4e'^ - l)^j. 

Case 11: 

P(X2 G iV^,(Xi,£) nn(Xi,£),Xi G T. \r(y,£)) ^ r r'-' ^(^(^.^2,Q.iV3,L.,L.,La)) ^^^^ ^ 

4r^2 _^ ig^ii _ gg^io _ 260 +372 + 1248 r'^ + 112 r*^ - 2624 r'^ - 8256 + 12288 + 13568 r2 - 27648r + 11264 
384 (16 r2£4 - 8 r2£2 + r2 + 64 r £* - 32 r £2 + 4 r + 64 £-* - 32 £2 + 4) r2 

where A{.^{Li, L2, Qi, N3, L4, L5, Le)) = [2 (-£2\/3 + 1/4 ^/3)^ (-72 + 24 \/3y r - 72 ^/3y x - 216x2 - 72xr+ 180x + 

72 r + 24 r2 y2 + y* - 6 r''x2 y2 + 6 r"* y2x - 32 y^ - 72 x^ + 180 x^ + 12 VSy + 36 y^x - 24 y2 + 24 r2y ^x + 56 VSy^x + 
24 y x^^ + 108 r2x^ + 72 r x^ + 2 r^'^Sy^ + 6 r*y VSx - 6 r^y^x'^y + 9 r-^x2 + 3 r"* y2 - 18 r^x^ + 9 r^x^ - 36 r2^x2y + 
12 VSx^y + 24 r^x^^^y -8r^V3y^x-24:ry^ + 72r y^x - 12 r'^y ^/3 + 108 r2x - lUr'^x^ - 36r^ - 72r x'^ - 12 r^y^/i- 
4r^V^y^ +48 ry^x- 72 ryx2y3-36 y22._|_ ^^2 y2x + 36 y x2v^+ 12 r-2^y=' + 4 r-2 y* - 72 y2x2 -36 r-2x'* -36r^x^ - 
36r^x + 72r-^x^ -44V3y^ - 8 r %/3 y^)] / [3 r2 (\/3y + 3-3x) (-y - \/3 + \/3x) (4 £2 - l)^]- 

Case 12: 

P(X2G7V];,(Xi,£)nn(Xi,£),XiGTAr(y,£))=( / / +/ / ''A A{nL.,L, Q, Q2,L,,U)) ^^^^__ 



Ar(y,£))= / / +/ / 

r®-71704r^-87608 

134144j j 1^10368 (r^ + 4 r + 4) (I6 £'' - 8 £ V l)] 



[135 r^^+540 r^^-2025 r-^°-8100 rVl0152 r*+38448 r^-14878 r®-71704 r^-87608 r''+192128 r^ + 147712 r^-338944 r+ 



where ^( J^(Li, La, Qi, Q2, ^5, ^e)) = " (-e^^+ 1/4^)^ (-18 + 12 r - 90 + 36 + 54 2; - 18 j/^ + 18 + 

r* y*-&r*x'^ y"^ + &r* y'^x-d.y* -ZQ a;* + 90 a;^ - 6 + 18 y^x- 24r^ - 12 r^y 73^+ 12 ^3 y^x + 12 y a;^^+72 r^x^ + 
36rs^ + 2r'*V3y^ + 6r'*y ^x-6r''^x^y+9r-*a;^ + 3r-* y^-18r'*a;^ + 9r'*a;'' + 12r^\/3a;^y+12r^x-^ y^ + 24 r y^- 12 r y^x + 
12r^y V3+36r^x-90r^x^-72rx^-12r^y\/3-4r^^y3-12r-yx^^/3-24r2 y^x-+12r^ y^x-6 y x2y3-2 y*-12y^x^- 
18 r^x'' - 36 r^x^ - 36 r^x + 72 r^x^ - 14 y^ + 4 r VSy^*) j / [3 (-^y - 3 + 3x)(-y-V3 + VSx) (4 - l) 

Case 13: 

P(X2 eiV^.B(Xi,e)nn(Xi,e),Xi eT, \r(y,e)) = 



'fl3 J rs{x) Jsii J rsix) J ai2 J r^ix) 



3654 r^^ - 35328 + 94802 - 100608 - 255 + 138240 - 193581 + 148224 - 86387 r'' + 43008 - 12369 + 512 

7776 r6 (r-2 + 1)^ (16 £-1 - 8 £2 + 1) 

where ^(^(Gi, Mi, La, Qi, P2, M3, Ge)) = " [4 (-£^\/3 + 1/4 \/3)^ (-9 + 42^/3yx - 45x2 + 36x - 15y2 + 21 r^y^ + 

2 r^'y'' - 12 r-'x^ y2 + 12 r*y2x + 18 x^- 12 ^y + 42 y^x- 24 y^ -6 r^y ^x + 4 y3y^x + 12 y x^y3 + 54r2x^ + 4r'*^y^ + 
12 r*y VSx - 12 r"* VSx^y + 18 r^x'^ + 6 r'' y^ - 36 r*x^ + 18 r^x"^ - 18 V^x^y + 12 r^^x^y + 12 r^x^VSy - 4 y^x + 
12 r^y \/3 - 45 r^x^ + 9r^ - 12 r^y ^/3 - 4 r^^/3y^ - 18 r^y^x + 12r-^y2x - 42y x'^VS + Gr^VSy^ + 2r^y* ~ 2Ay^x'^ - 
18r2x* -36r^x^ -36r=*x + 72r^x2 -2x/3y^)]/[3r2 (-v^y- 3 + 3x) {-y - + V^x) {4e'^ - l)^]. 

Case 14: 

P(Xa . N^MXus)nrUX.,s),X. , TAT(y,£)) ^ ( H f'^^ + /^^^ r'^') MnC. M., L.Q.N. Ma, M ^ 

49 r-^2 + 124288 + 50688 + 384 r" - 3562 r^" + 13440 - 36948 t/** + 27648 - 86016 - 1024 - 89088 + 160 r"^ 

15552 r6 (16 £-1 - 8£2 + 1) 

where A{^{Gi, Ml, L2,Qi,N3,Mc,M3, Ge)) = - [2 (-£^^3 + 1/4 ^/3)^ (-18+84 ^/3yx-90x2+72x-30y2+38r2y2 + 

r^y* - 6 r'*x2 y^ + 6 r^y^x + 36 x^ - 24 ^y + 84 y^x - 24r^ y^ - 4 r^y ^/3x + 8 \/3 y^x + 24y x^^/3 + 108 r^x^ + 2 r*V3 y^ + 
Gr'^yVSx - Gr'^VSx^y + 9r''x2 + 3r''y2 - 18r*x^ + 9r''x^ - SGr^VSx'^y + Ur'-^VSx^y + 24r^x^V3y - Sr'^VSy^x + 
16 r^y \/3 + 24 r-^x- 102 x'^ + 6 r'^ - 12 r^y ^ - 4 r■^^/3y^ - 36 y2x + 12 r^y^x - 84y x^x/S + 12 r2^/3y^ + 4 y* - 
48y2x2 -36r2x'' -36r='x^ - 36 r^x + 72 r-'^x^ - 4 \/3 y^) j / [3 (-^/3y-3 + 3x) (-y - \/3 + \/3x) (4 £2 ^ l)^]' 

Adding up the P{X2 € Np^{Xi,e) n ri(Xi, e), Xi e \ r(y, e)) values in the 14 possible cases above, and 
multiplying by 6 we get for r € [1,4/3), 

g (r-l) (5 + 288 r^e^ + 1152 r^e^ - 148 + 1440 r^e^ + 245 + 576 r^e^ - 178 - 232 r + 128) 

/^and('^,e) - 54 r2 (2 + r) (2£ - 1)' (2e + 1)' (r + 1) ' 

The Mfn(j('"'^) values for the other intervals can be calculated similarly. 
Derivation of yufj.(r, e) 

For r £ [1,4/3), there are 16 cases to consider for the OR-underlying version: Case 1: 



P{X2eN^,E{Xi,e)url{Xi,e),XieTs\T(y,e)) = 



; ri^am{x) rSQ r^am{x) PSi P^amK 

/ +/ / +/ / 

Jll(x) Jq2 JO Jso Jri(x) 



A{^{Vi,Mi, Mc, M3, Ve)) ^ _ 6e^ - Ar^ + 12r - 9 
A{T,)2 27(4£2-l) 



where 4(^(Vi, Mi, Mc, M3, H)) = ~ ' ^"''ffX'if 
Case 2: 

P{X2 GiVpB(Xi,£)un(Xi,£),Xi GTs\r(y,£)) = 

rsi rri(x) /-ss rr2(x) rs^ rr^(x) rsn /■r5{a;)\ 

Jsg Jo J SI Jo Js5 Jo J S3 Jr3(x) J 



A{T,Y 

-2304 r^£ ^432 r--21960 /-27+9624 r ^5952 r®£ ^288 r*£Vl824 r-**£2-1817 r*-2880 r ^10368 r^+28224 r^-5760 r V- 



21964 r*^ 



]/[864r« (l6£4-8£' + l)] 



where Afi, La, L3, Mc, M3, Ve)) = - [-27 + 12 e'^r'^x^ + 36 V3y r + 108 a; - 162 s:^ _ 108 a: r - 8 e^V^r'^yx + 

108 a; - 54 + 36 r - 5 - 3 y"* - 27 a;* + 108 - 36 ^y + 108 y^x + 10 r^y + 12 ^y^a; + 36 y a;^^ - 36 r a;^ + 
36 ry2 - 36ry2a: - 10 r^y ^ + 30 r^a; - ISr^a;^ - 15r^ + 108 rx^ - 72 ry^x + 36 r y + 12 r^e^ - 108 yx^^- 

54y^x^ - 12^y=' + 4r^y^ +4eV^y2 - 24 eV^x + 8 e^^r^yj / [4r2 (-^y - 3 + 3x) (-y - ^/3 + ^/3x)j. 

Case 3: 



- [-3456 r^e^ + 1296 r - 65772 + 26880 + 9216 r^s^ + 432 r^'e^ + 3072 r^e^ - 4864 - 8640 + 31104 + 83808 r^- 

9216 r'e^ - 63744 r** - 8l] / [2592 r*^ (l6 e"* - 8 + l)] 

where A(^(Vi, Mi, La, L3, 1/4, Lg, Mg, H)) = - [-54 + 12 eV^x^ + 36 \/3y r + 54 x/3y x - 189 x^ - 180 x r - 8 e^v^r^y x + 

162 X - 27 y^ + 72 r - 9 y^ - 15 y" - 27 x'^ + 108 x^ - 18 V3y + 108 y^x + 18 r^y VSx + 36 ^y^x + 36 y x^ ^3 - 36 r x^ + 
12ry^x - 18r-^y\/3 + 54r^x - 27r^x^ - 27 + 144rx^ - 48ry^x + 12ryx^^+ 12r^£^ - 72yx^^- 90y^x^ - 
24 ^y^' - 4r 73 y^ + 4e^r'^ y^ - 24£V^x + 8 VSr^yj / [4 (- V3y - 3 + 3 x) (-y - ^ + ^x) j . 

Case 4: 

PiX. . N^MX.,s)urUX.,s),X. . TAr(y,.)) ^ £"^£ ^(.^(1-.C.,(.3.M.,M.,M3,1-.)) ^^^^ ^ 

-8 r + 24 + 56 r'^ - 32 r-^ - 13 + 32 r^e^ + 32 r^e^ - 128 r^e^ + 192 r^e^ - 128 r^e^ + 64 - 92 r*^ + 1 

96 r6 (4e2 -1)^ 

where ^(^(Vl, G2, G3, Ma, Mc, M3, Ve)) = - ^y''+<^y-<^y^+3 ^g:.+3 v^r^e^ ^ 

Case 5: 



P{X2€N'pE{Xi,e)url{Xi,£),Xi€Ts\T{y,£))=[ / +/ / 

Jr^ix) Jae Jo J 



A{^{Vl,G2,G3,M2,Mc,P2,N2,V6)) , , 

-1 + 32r='e^ + 5r + 34r-* + 15r' +64r''e^ - 32r*e^ - 32r"e^ - 17r'' + 29r" - 3r-^ - 17r^ - 2r^ - 64r-''e2 - 43r'' + 32r^e^ 

96 (r + l)^(4e2-l)V6 

where ^(^(Fi,G2, 6*3, M2,Mc,P2,iV2,V6)) = - [2 (-£^^/3 + 1/4 ^3)^ (2 V3y^ + 12y-12yx+6^/3-12^x+6^x2- 
7^/3r2 + 12r-^y + 12r^^/3x -4r^\/3y^ - 24 r^'yx - 12 r''\/3x2 + 8 \/3r^e^) j / [9 (4e^ - l) 
Case 6: 



P{X2 G N'pe{Xi,£)VJTI{Xi,£),X^ G T, \r(y,e)) = 

t.r.(x) ^ rs, rt^^ix) ^ ..13 rrs{^)\ {Vl , N, , , L2 , L3 , L4 , L5, P2 , iV2 , Vg)) ^^^^^ . 



243-28578916r' -1147392 r'^e -1344384 r-£ -1734912 r"e -304128 e-" +989424 r-"e -10368 re^+3888r-438777 r + 



[- 

2204160 - 355328 r^* + 5753232 + 39312 r^e^ + 1296 /e^ + 296208 r^e^ - 20639832 r^"* + 13254912 r^^ - 6591792 r-" + 
1693728 r^^e^ + 1507392 r-^e^+637056 r^'^e^+26417664 r^^+26760576 r"-21960774 r^°+15877152 r^-10180620 r**-28107 r^ + 
128304 r^ + 1222128r^-120960r-''e^-2856483r®+92160r^®e^-563328r''e^jy^ [7776 (r^ + l)^ (l6e'' - 8e^ + l) (2r^ + l) 

where A(^(Vi, iVi, Pi, L2, L3, ^4, is, P2, iV2, 14)) = " [4 (-£^^3 + 1/4 \/3)^ (-54 + 12eV2x2 + 36\/3yr- + 54^/3yx 
189 x^ - 180 X r - 8 £^^/3r^y x + 162 x - 27 y^ + 72 r - 12 y^ - 4 r** y** + 24 r'^x^ y^ - 24 y^x - 15 y^ - 27 x'' + 108 x 



3 

2„ 



18V3y + 108y^x + 24r^ y^ + 24r2y^/3x + 36^y^x + 36yx^\/3 - 36rx^ - 8r■^^/3y^ - 24r''y^x + 24r''^/3x'y 
36 r*x^ - 12 y^ + 72 r^x^ - 36 r^x"* - 12 VSx^y + 12 r y^x - 24 r^y ^3 + 72 r^x - 36 r^x^ - 36 + 144 r x^ + 12 r^y v% + 
4 r^A/3 y^ - 48 r y VSx + 12 r y x^Vs - 12 y^x + 12 r^e^ - 72 y x^ ^3 - 90 y^x^ + 36 r^x^ + 36 r^x - 72 r^x^ - 24 ^3 y^ - 
4 r V3 y^ + 4 eV^ y^ - 24 eV^x + 8 VSr^y)] / [s (-\/3y - 3 + 3 x) (-y - v/3 + V3x) (4 - l) 



Case 7: 



\ny,e))=[ / +/ / 



P{X,eN},,{X,,e)urliX,,e),X,eTATiy,e))=\ I / +/ / | MnV., N,,Q,,L,M, L,, P,, N,, V,)) ^^^^ ^ 



-[8(-2-55766r^^-864r-^^e^-4104r"e^-3024r"e^+3690r^°e^-108r''^e^+24r-1833 r''+21576 r-''+342 r®e^ + 18r*e^ + 
1710 r**e^-20056 r^'*+6912 r-^^-1152 r^'^+SSie r-^^e^ + 1800 r"e^+288 ^^^^£^+38376 r"+65532 -63642 r-^''+52020 r-'-'-36277 r*- 
142 + 576 + 4848 - 864 r^e^ - 10994 r'^ - 2700 r^e^)j / ^243 (2 r'^ + if (r^ + l) ^ (16 e"* - 8 £^ + l)] 

where yl(^(Vi, iVi , Qi, L3, ^4, -Ls, P2, iV2, V'e)) = - [4 (-e^^+ 1/4^)^ (-36 VSr^ - 180 V3r x + 12 ^r^e^ _ gg^a^ _ 

30 ?/2 _^ 18 _^ 54 VSa; + 108 r y + 54 y a; - 135 y3a;2 - 9 73 - 45 ^/3 + 72 r^y a; + 126 - 60 + 72 + 

144 v^r x^ - 18 y + 18 - 36 x^ VSr + 36 xV^Vs - 3 r^Vs + 54 r^x^y - 72 ^r^x^ + 24 ^/3r^ + 54 r'*a:^^/3 - 

UAyrx- 36r^x^y - 12y^r + 96y^x - 18r^y^ - 18a;'* ^3 + 12r^y^ + 36r^y - ISr^y'^VSx - lOSx^y + uVSy^rx- 
12 r^VSx + 54 r^x^y - 54 r'^x^y/3 + 72 y^VSx - 9 r^'V^y^ - 54 r^y x + 36 r^V^x - 27 r^^x"^ - 2 y'^r'^^fl - 18 r^y^ + 
18 y^\/3s+72 a:^^/3-72 y^^/3a;2+24 eV^y-27r*y3x'*+12 y^a;-36 r^a:^y-14 ^"73+72 a;^y+36 x^r y+18 y2a;^ + 
4^2 y^ + 12e2^r2s2 _ 24e2^r2a; - 24e2r2y a;) j / [3 (-^y - 3 + 3 a;)^ (4 - l)^]. 

Case 8: 

P(X2 G iV^.E(Xi,e) ur5:(Xi,e),Xi G T, \r(y,e)) = 

[-81 r V 189 - 561 r'^ + 1008 r^e^ + 45 - 432 r^e^ + 1894 - 3120 r^e^ + 18 - 144 r'*e^ - 1912 + 2304 r^e^ - 224 + 

768 r^£^ + 384 r + 128] / [l296 (16 - 8 + l) (r + I)''*] 

where yl(^(Vi, iVi, Qi, G3, M2, Mc, P2, Af2, Ve)) = [2 (-£"^+1/4^" (-^r a; - 24 r^y - 8 r^^y^ + 24 r'VSa;^ - 

24 v'Sa; + 36 r^y x + 8 e^VSr a: - 8 e^VSr + 25 r y + 24 y a: - 12 ^/3a:^ - 4 \/3y^ - 8 y - 12 ^/3 + 24 A/3a; + 13 VSr - 
24 A/3r a;^ + 8^3 y^r- - 24y + 12 a:^\/3r - 18 xV^a/3 + 18 ^r^x-^ + 6 \/3r^ y^ - 18 r^x^y + 4 y^r + 6 y^ -4V3y^rx + 
6r='^/3xy2 - 12xVy)j/[3r (^/3y + 3-3x) (4£2-l)^]. 

Case 9: 

P(X2 G iV^.B (Xi , e) U (Xi , e) , Xi G \ r(y , e) ) = 

S4 ^rg(a;) ^si4 Jrg{a;) ^si2 J riQ(x) ; ^(r,)2 

[-512+81297 r*^ + 55296r"e^-51264r"'e^ + 6912r^£^ + 6144r + 72576 1512 r**-798720 r^-35424 r''£^-9216 r'^e^- 
45792 r**e^ - 83538 - 17280 r + 18252 r " - 51840 r^^e ^6912 r"e^ + 167616 r " - 565920 +888957 r^" - 1023600 r^ + 
998852 r** - 7424 r^- 6144 r-^- 262080 r^ + 41472 r^e^ + 533036 rV82944 r^e^j / [5184 r'* (2r^ + l)^ (l6e'' - 8e^ + l)] 

where A{.^{Vi,Ni,Qi,L3,L4,Q2,N2,V6)) = - [8 {-e^VS + 1/4^3^ {-18V3r'^-90V3rx+6V3r^e'^-54r'^y-15r^V3y'^+ 

27r2v/3x^ + 18r^^x + 54ry + 54yx-63^x^-9^y^-18y3 + 54r2yx + 54y3x-24y-'+36\/3r + 72 V3rx^-18y + 
9 x^r^ \/3 - 18 x^ V3r + 18 xV^ ^3 - r*^ y* + 18 r*x^y - 36 VSr^x^ + 12 VSr^ y^ + 18 r*x^x/3 - 72 y r x - 18 r^x^y - 6 yV + 
36 y^x - 9x'*V3 + 6 y^ + 18 r=*y - 6 r*y^^x - 72 x^y - 6 \/3r'^ y^x^ + 6 \/3 y V x - 6 r^^/^x y^ - 36r^x^%/3 + 36 y^V^x - 
3r* 73 y^- 18 r*yx + 18 r=*\/3x-9r''^x^+y''r2^-6r* y^ + 12r^ y^^/3x + 36 x^y3-30y^^/3x^ + 12eV^y-9^''^/3x''- 
5y''\/3+36x=*y+18xVy+6r''\/3y^x^+2e^\/3r^ y2+6 e^\/3r2x^-12 £2\/3r^x-12 eV^y x) j [3 {V3y + 3-3xf (4e^ - l)^]- 

Case 10: 

PiX2 G ;v;.(X„.) un(X„.),X. G T. \r(y,.)) ^ /^^^ /^^'^' ^(^(K.A^.Q. M2,iV3.iV2.V^a)) ,^,^ ^ 

2496 r"* + 1728 r^e^ - 4608 r'^e^ + 512 - 81 r^" + 270 - 2176 + 3072 r^e^ - 1080 

5184 (16e4 _ Se^ + 1) 

where A(i32i(Vi, iVi, Qi, G3, Af2, A^a, Af2, Ve)) = -[2 {-e^V3 + 1/4^^3^ {3 V3r x - I2r^y - 4r^V3y' + Ur^VSx^ - 
12r^V3x + 18 r^yx + 8e^V3rx - Se^VSr + 21 ry + 24 yx - 12 v^Sx^ - 4 VSy^ - Se'^ry- 12 ^3 + 24 VSx + 9 ^/3r - 



Sr^VSxy^ - 12 J/)] / [3 r {-VSy -3 + 3x) (4£^ " l)^]- 
Case 11: 

P{X2 e N'pE{X-L,e)url{Xi,e),Xi eTs\T{y,e)) = 

/ rsr, ^ rs,, rrs(.)\ ( , iVi , A , La, L3 , Mc , P2 , jV2 , Ve)) ^^^^ _ 

-43855 r^^ + 14112 r^^e^ +271488 r" -48384 r"e^ - 746553 +81792 r^°e^ - 117504 r^'e^ + 1230336 - 1404177 + 
123840 r®e^ + 1236528 r'' - 89856 r^e^ - 901350 + 58752 r'^e^ - 20736 r^e^ + 550800 - 276453 + 2592 r^'e^ + 104976 - 

26649 + 3888 r - 243] / [7776 r'' (16 e'' - 8 + l) (r^ + l)'^j 

where A{^{Vu Ni, Pi, L2, L3, Mc, P2, N2,Ve)) = " [4 (-£^^3 + 1/4 (-27 + Ue^r'^x'^ + 36V^yr + 108 - 

162 _ 108 X r - 8 VSr^y a; + 108 x - 54 j/^ + 36 r - 8 - 4 r"* y"* + 24 r^x^ - 24 r"* y^a; - 3 - 27 x* + 108 x^ - 
36^3^ + lOSy^a; + 24 y'^ + Wr^yVSx + uVSy^x + 36yx^V3 - 36rx^ - Sr^VSy^ - 24r^y^f3x + 24r-*^a-^y - 
36 r*a;2 - 12 t/^ + 72 r'^x^ - 36 r^x^ - 12 + 36 r - 36 r y^x - 16 r^y ^3 + 48 r^a; - 24 r^x^ - 24 + 108 r + 

12 r^y \/3 + 4 v^Sj/^ - 72 r y ^/3a: + 36 r y x^^- 12 y^x + 12 r^e^ - 108 y x^v^- 54 y^x^ + 36 r^x^ + 36 r^x - 72 r^x^ - 
12 ^ y^* + 4 r V3y^ + 4 eV^ y^ - 24 £ V^x + 8 e^^r^y)] j [3 (-^y - 3 + 3x)(-y-^/3 + ^x) (4 - l) ^ j . 

Case 12: 

- [5184 r-**e^-71424 r*'e^ + 138240 r^e^-73728 r'*e^-1053 r^^ + 16230 r"^-17856 r-^-68908 r* + 104448 r-''+276688 r''-916608 r^ + 

1032192/ - 516096 + 80128 + 12288 r - 1024j / ^31104 (16 - 8£^ + l)] 

where ^(i3^(Vi,A/'i,Qi,L3,Af3,iV2,V6)) = - [2 (-e^VS + 1/4 \/3)^ (-36 \/3r2-216 \/3r x+24 \/3r2e2-108 r2y-48r2^/3y^ + 
72r2^x2 + 36r-2^x + 216ry + 432yx-216^x2-72 V3y^-36x/3 + 72r2yx + 144y3x-48y-' + 72^r + 216^rx2 + 
72^/3yV - 144y + 36x^r^^ - 72 x^VSr + 36 xV^^S - 3 r^VSy"* + 54 r*x^y - 72 \/3r^x^ + 24 ^/3r^ y^ + 54 r'^x-'^x/S - 
432yrx - 36 r^x^y + 24 yV + 48y^x - 36 r^y^* - 36x*v^+ 12 r^y^ + 36 r^y - ISr^y^VSx - 432 x^y - ^2^/3y^rx - 
12 r^VSx y2 + 108 r^x^y - 108 r^x^* ^3 + 144 y^^x - 9 r^^fly^ - 54 r^y x + 36 r^VSx - 27r* VSx^ - 4 y^'r^ ^3 - 18 y^ + 
36 r^y^ySx + 144 x^^ - 72y^y3x^ + 48eV2y - 27r''v^x* + 24 r^y^x - 72 r^x^y - 4y'^^/3 + 144 x^y + 216xVy + 
ISr^ySy^x^ + 8 e^ySr^y^ + 24 e^ysr^x^ -48 e^ySr^x- 48 eV2yx)j/ [3 (-^Sy - 3 + 3x)^ (42^ _ ij^j. 

Case 13: 



P(A2 G Afpisl-^i.e) uri(Ai,e),Ai e Ts \r(y, e)) = / / dydx = 

[-13r^^-78r^V42r"+892r^''+64r-^eV220 r^-4952 r-^+384 r**e^-768 r^-3072 r^eV 18048 r**- 3136 r^- 2048 r^eV 
8192 r'^e^ - 39296 + 20992 + 4096 r^e^ + 41984 - 8192 r^e^ - 48128 r + 14336] / [384 (I6 r^e"* - 8 r^e^ + + 96 r^e"*- 

48 r^e^ + 6 + 192 r e"* - 96 r + 12 r + 128 e"* - 64 + 8)r^] 

where^(^(Vi,Ari,Qi,L3,Ar3, A^2, 14)) = -[2 {-e'^VS + l/^Vsf {~36V3r^-216V3r x+24V3r''e''-108r^y-48r'^V3y^+ 

72 r^ySx^ + 36 r^y3x + 216 r y + 432 yx- 216 73x^-72 73 y^- 36 73 + 72 r^yx + 144 A/3x-48y^ + 72 + 216 ^rx^ + 
72^/3yV - 144y + 36x'*r2 73 - 72 x^VSr + 36xV^\/3- 3 r^^Sy* + 54 rVy - 72 VSr^x^ + 24 ^/3r^ y^ + 54r*x^^- 
432 y r X - 36 r^x^y + 24 y^r + 48 y^x - 36 y^ - 36 x''^/3 + 12 y^ + 36 r^y - 18 r*y^^/3x - 432 x^y - 72 ^Sy V x - 
12 V3x y^ + 108 r^x^y - 108 r^x^%/3 + 144 y^^x - 9 r^VSy'^ - 54 r^y x + 36 r^VSx - 27 r''^/3x^ - 4 y''r^^/3 - 18 r'' y^ + 
36 r^y^^x + 144x^^3 - 72y^^x^ + 48eV2y - 27r''^x* + 24 r^y^x - 72 r^x'^y - 4y'^V3 + 144 x^y + 216xVy + 
18r■^^/3y^x2 + 8 e^VSr^y^ + 24 e^VSr^x^ -48 e^VSr^x- 48 eV^yx)]/ [3 (-^/3y - 3 + 3 x)^ (4£^ - ^Y]- 



Case 14: 



P{X2 e N'pE{Xi,£)url{Xi,e),Xi G%\T{y,e)) = 

S14 •''•12(2;) Jsis Jrio(a;) 

-[-189r^^-1134r^^+297r" + 11718r^''+864r^e^+3672r^-66096r*+5184 r**£^+2592 r^-32832 r^e^+248616 r*^ 

30448 r^-33408 r^e^+76032 r'*e2-551584 r'*+273152 r^+55296 r^e^+595456 r^-73728 r^e^-668160 r+197632] / ^5184 



3 



+ 6r^ + 12r- + 8) (16 e"* - 8 + l) j 



where A{3^{Vi,Ni,Qi,L3,L4,Q2,N2,V6)) = - [8 {-e'^VS + I/aVs) {~18V3r^-90VSrx+6V3r^e^-54:r^y-15r'^V3y^+ 

27 r^VSx^ + 18 r^ySa; + 54 r y + 54 y a: - 63 VSa;^ - 9 \/3 - 18 73 + 54 r^y a; + 54 ^/3x - 24 y3 + 36 %/3r + 72 ^/3r - 18 y + 
9 a;*r^y3 - 18 a;^ v^r + 18 x^r^VS- r'^V^y'^ + 18 r^x^y - 36 VSr^x'^ + 12 V3r^ + 18 r'^x^x/S- 72 y r x - 18 r^x'^y - 6 y V + 
36 y^x - 9 + 6 y'^ + 18 r^y- 6 r^y'^VSx - 72 x^y - 6 v^r^ y^x^ + 6\/3yVx-6r^ a/5x y^ - 36 r^x^ ^/3 + 36 y^VSx - 
3r* V3 y^- 18 r*yx + 18 r^V3x-9r*y3x2+yVy3-6r* y^ + 12r2 y2y3x + 36 x^ V^- 30 y^VSx^ + 12 e^r^y- 9 r^^V^x^ - 
5 y*^+36 x3y+18 xV y+6 r* VS y2x2+2 £273^2 y2+6 e'^y/3r^x'^-12 e^y/Sr'^x-U eV^y x) j / [3 (-^y - 3 + 3 x)^ (4 _ 1) 

Case 15: 



P{X2eN^PE{Xi,e)url{Xi,e),XieTs\T{y,e)) = 

= 11 /■r2(:.)^ /■.12 /-2(.)^ rso rre(.)\ A{.^{Vl,Nl,Ql,L3,Mc,P2,N2,V6)) _ 



[4536r^^e^-11753r^^-13824r-"e^+69120r"+23976r"^e^-186683r"-34560r^e^+305664 r^+35496 r^e^-346171 
27648 r^e^+302592 r^+17208 r'^e^-220201 r®-6912 r^e^+135936 r-^+1152 r''£^-69760 r'*+28416 r^-8384 r^+1536 r- 128] 

[l944r« (r2 + l)'(l6e^-8£' + l)] 

where A{^{Vi, Ni,Qi, L3, Mc , P2, N2,V6)) = - [4 (-£^\/3 + 1/4 ^/3)^ [-24, y/?,r'^ ~ 108\/3rx + 12y/i,r'^e^ - 66r2y - 

26 y2 + 30 r^ySx^ + 30 \/3x + 108 r y + 216 y X - 108 ^x^ - 36 V3 y^ - 18 V3 + 48 r^y X + 72 ^/3x - 24 y^ + 36 %/3r + 
108y3rx^ + 36y5yV - 72y + 18 x^'r^v^ - 36 x^^/3r + 36 xV^Vs - 3 r''^/3 y"* + 54 r-'x^y - 12y/^r^x^ + 24 ^Sr^ y^ + 
54r'*x^y3 - 216yrx - 36 r^x^y + 12 yV + 24y'''x - 18 r^y^ - 18x^v^ + 12 r^y^ + 36 r^y - 18r*y^^/3x - 216 x^y - 
36 VSy V X - 12 r^VSx y^ + 54 r^x^y - 54 r^x^VS + 72 y^^/3x - 9 VSy^ - 54 r^y x + 36 r^^Sx - 27r*y3x^ - 2 y^r^v^ - 
18 r'' y3 + 18 y^ v^x + 72 x^^/3 - 36 y^VSx^ + 24 eV^y - 27 r*^x* + 12 y^x - 36 r^x^y - 2 y''y3 + 72 x^y + 108 xV y + 
18 r-'^y^x^ +4 e^^r^ y^ + 12e2^r^x2 - 24e^^r^x - 24eV2y x)] y' [3 (^y + 3 - 3x)^ {Ae'^ - l)^] • 

Case 16: 

PiX2.NWiX.,e).VliX.,e),X.eT.\ny,e))4 f f'^ /^^^ f'l -^(^(^^■^^■Q^^3^^3,iV.,l^.)) ^^^^ 



S12 J re{x) Jag Jr3(x) I 

147 r^^+55296r^e^-12288r-+351872r-*-142080r-'^+1024- 73728 r-®e^-9216r'*e^+576 r**e^-1152 r-"+7018 r-^''-20352 r'' + 
51188 + 64000 - 190464 - 414720 + 27648 e' + 305920 r*^! / [l5552 r" (I6 e"* - 8 + l)] 



where A(^(yi,Ari,Qi,L3,iV3,A*'2,V6)) = - [2 {^-e' ^^l + Ij A^PS) [^m^r^ -21% ^/Ir x+2A^/It^ -1^8 r'y- 18 ^y^ + 
72 r^^x^ + 36 V3x + 216r y + 432 y X- 216 VSx^ - 72 ^3 y^- 36^+72 r^yx + 144 V3x-48y3 + 72 + 216 ^rx2 + 
72 VSy V - 144 y + 36 x^'r^^ - 72 x^ ^Sr + 36 xV^^ - 3 r'^^fZy^ + 54 r*x^y - 72 v^r^x^ + 24 VSr^ y^ + 54 r^x^^fZ - 
432 y r X - 36 r^x^y + 24 y^r + 48 y^x - 36 y^ - 36 x"* ^3 + 12 35 ^3^ _ ^4^2^^ _ ^2^ _ 72 ^Sy V x - 

12 r^^x y^ + 108 r^x^y - 108 r^x3^/3 + 144 y^VSx - 9 r^'^y^ - 54 r^y x + 36 r^^x - 27 r''^/3x^ - 4 y^r'^^- 18 y^ + 
36 r^y^^x + 144 x^v^^ - 72y^\/3x^ + 48eV2y - 27r''^/3x** + 24 y^x - 72 r^x^y - 4y^^/3 + 144 x^y + 216xVy + 
18 r-*^y2x2 + 8 e2^/3r2y2 + 24 e^^/Zr'^x^ ~ 48 e^^r^x - 48 eV^y a:)] / [3 (^/3y + 3 - 3 x) ^ (4 - l) . 

Adding up the P{Xi S A''p^(Xi, e) U T\{X\^ e), -'^i G 7s \ r(y, e)) values in the 16 possible cases above, and 
multiplying by 6 we get for r € [1,4/3), 



47 - 195 + 576 r^e^ + 860 - 288 r^e"^ - 864 r^e^ - 846 + 1728 r^e'' - 108 + 1152 r^e^- 



576r2e2 + 720r-256 



108 (2 + r) (l6 e"* - 8 + l) (r + 1) 



The /xfj.(r, e) values for the other intervals can be calculated similarly. 

For r = oo, it is trivial to see that ^(r) = 1. In fact, for fixed e > 0, /i(r) = 1 for r > \/3/{2e). 
Remark 7.1. Derivation of y^and(^'^) ^^'^ t^ori''^j^) is similar to the segregation case. 

Appendix 7: Proof of Corollary I6.lt 

Recall that S'™'^(r) = pY\f{r) is the relative edge density of the AND-underlying graph for the multiple triangle 
case. Then the expectation of 5,|^"'^(r) is 

But, by definition of iV];^(-) and r^(-), if and X2 are in different triangles, then P{X2 £ N},j^{Xi)nTl{Xi)) = 
0. So by the law of total probability 

A^andW := Pix^e n'pe{Xi) nriix^)) 
./™ 

= J2 ^(^2 e vVp£(^i) n riiXi) \ {Xi, X2} c t,) ^2} c t,) 

i=l 

Jm 

= 5I/^a„dWP({Xi,X2} C T,) (since P{X2 G iV^B(Xi) n r^^l) I {^1,^2} C T,) = MandW) 



^1=1 



= MandW E' 

where /iand('') is given by Equation ([8|). 

Likewise, we get /ior(f ) = fJ-oi-{r) wf^ where fiorir) is given by Equation 

Furthermore, the asymptotic variance is 

?and(r) = E [hlfir)hlfir)] - E [hlfir)] E (r)] = P ({X2, X3} C A^j;i;(Xi) n rU^i)) - (Aland(r))' . 
Then for Jm > 1, we have 



1=1 

= p({X2,X3}c7V^i;(Xi)nruxi)|{Xi,X2,X3}cre) 

Hence, 

J^andW = P({X2,X3} C N^pEiX,) H Tl{X,) \ {X,, X2, X3} C T,) f ^ ^ " (/landW)' 
= J^andW fxi^^^l +Aiand(r)' I E " (E ^' ^ 



i=l \j = l 



Likewise, we ge. . (sfc «,?) . ..(.)^ (E.t, »f - (St. 

So conditional on 3^„i, if T'andir) > then {Sf^^{r) — Atand(?')) — > A/" (0,i^and (»"))• A similar result holds 
for the OR-undcrlying version. 

Appendix 8: Proof of Theorem 16. 2t 

Recall that pfj'^^ir) is the version II of the relative edge density of the AND-underlying graph for the multiple 
triangle case. Then the expectation of pfj'^nir) is 



: ~ 1) 
2nt 



E 



since by ([J) we have 



ni{n.i - 1) 



k<l 



where /Ltand('") is given by Equation ([8]). Likewise, we get Por{r) = Poi{r) where Ho-cir) is given by Equation ([9]). 
Next, 



Var[pr/lM] =E 



i=l 



Var 



since and are independent for fc 7^ Z. Then by ^ we have 



Var 



^ -Var [/.-'(r)] + ^ Gov [/.-'i(r), /i-<^(r)] 



; - 1) 



So, 



Var [p- (.)] = E Var [.-(0] + £ ' ' 



Cov[;itf(r),/itf(r)]. 



1=1 



Here X;/=i ^"^a"^^ Var [/i?5'*(r)] = Var [/ifa'"*!^)] ■ Then for large n, and n, 



-Var [/i-'M] 



1^ E/=i 



Var 



since ^ = 2?!^ ^'^'^ '^i/?^ — > as , n — > cxd. Similarly, for large Ui and n. 



i=i 



Cov[/itf(r),/.-'i(r)] 



Hence, conditional on 3^„i, ~ Mand('")) — ^ A/" (0, 4 i>and('')) provided that i^and('') > where 

/iand(r) = A^and(?') and i>and(r) I'andl^') {j2i=i ^i) j {j2i=i ^i) ■ A similar result holds for the OR-underlying 



version. 




Figure 23: The vertices for Np^[xi,e) D r^(a;i,£) regions for xi G Ts in addition to the ones given in Figure 
[23] because of the restrictive nature of the alternatives. 




Figure 24: An illustration of the vertices for possible types of Np^{xi) n T\{xi) for xi £ Ts- 




Figure 25: Prototype regions Rt for various types of Np^{xi) D T'{{xi) and the corresponding points whose 
cc-coordinates are Sk values. 




Figure 26: Prototype regions Rt for various types of Np^{xi) D T'{{xi) and the corresponding points whose 
cc-coordinates are Sk values. 



